BREACH EROSION PROCESS OF 
HOMOGENEOUS EARTHFILL DAMS 
AND FUSE PLUGS DUE TO 
OVERTOPPING FLOW 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 
DOCTOR OF PHILOSOPHY 


by 

PRADIP KUMAR DAS 


to the 

DEPARTMENT OF CIVIL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

DECEMBER, 1997 


liiMM 


'2'1 

temn LtBRAif 

I. I. T ^AWfUi 


12g.596 

TH 

(^Ejjw/jp 

D9$b 



CERTIFICATE 



It is certified that the work contained in the thesis entitled “BREACH EROSION 
PROCESS OF HOMOGENEOUS EARTHFILL DAMS AND FUSE PLUGS DUE 
’I’O OVERTOPPING FLOW”, by Pradip Kumar Das (Roll No. 9410366), has been 
carried out under my supervision and that this work has not been submitted elsewhere for 
a degree. 


(T Gangadharaiah) 
Professor 

Department of Civil Engineering 
Indian Institute of 'fcchnology 
Kanpur, India 



7 December, 1997 


iii 


SYNOPSIS 


Dams play an important role in surface water management. The importance of a dam is 
greater in a counity like India where precipitation is not uniform in time and space. However, 
when a dam fails, it results in a sudden release of large quantity of stored water causing havoc 
in the downstream region. The damages in terms of life and property associated with a dam 
failure can be catastrophic. So the failure of a dam in any form is not desirable. Among the 
dams, earthen dams arc more vulnerable to failure. The main causes of the failure of an 
earthen dam are due to overtopping and piping. In this ‘International Decade for Natural 
Disaster Reduction’, (199() - 2()()() A. D.), a natural disaster like the dam breach problem, has 
its relevance to study. 

The present study is dealing with the process involved in the dam breach and the washout 
characteristics of fuse plug due to overtopping flow. The results of the present study on dam 
breach have a direct bearing on the field, particularly for flood warning and evacuation 
planning measure as remedial efforts towards the probable disaster mitigation. In more 
technical temis, knowing the initial conditions at the dam site, the routing of water and 
sediment discharge on the downstream can be carried out more effectively. Fuse plug works 
as a safely valve of an earthen dam. The height of fuse plug is equal to the difference between 
full reser\c)ir le\'el(FRI ) and maximum reservoir level(MRL). This is provided intentionally 
tt) be \\ ashed-oui during high flood and to be reconstructed after the passage of flood. The 
objective of providing fuse plug in the design of earthen dam is to make the spillway much 
smaller in height to handle floods of lower return period compared to the maximum probable 
flood at substantial savings to the project cost. In an eventuality of failure of crest gates to 
operate in time it also serves as an emergency outlet. 


project cost. In an eventuality of failure of crest gates to operate in time it also serves as an 
emergency outlet. 

The causes of overtopping of flood over the embankment may be due to improper design 
criteria of determination of probable ma.ximum flood (PMF) for the reservoir, sedimentation 
of reservoir, adopting the erroneous seasonal reservoir operational mles, or due to any 
unexpected flash flood. 

Embankment breach process is three dimensional in nature where erosion occurs freely in 
vertical, lateral and longitudinal directions of flow. As such, understanding of the dam breach 
process i.e. breach hydraulics and breach morphology and breach erosion is of paramount 
importance to predict outflow hydrograph and sediment discharge. During the last 25 years, 
though substantial numerical models on dam breach problem have been developed, yet 
relatively much less attention has been given to study the mechanism of dam breach with 
experimental observation support. In an experimental set up the measurement of variation of 
one or more controlling breach variables with time and space can be made to study their 
effects on dam breach problem. The dam breach process is also related to sediment transport 
under highly accelerated flow conditions. The aspect of accelerated flow has not been 
addressed in almost all the numerical models. On the contrary, these numerical models rely on 
sediment transport formulae meant to be used in streams with uniform flow. 

There is not much literature on the wash out of a fuse-plug. Dam breach process and fuse- 
plug wash out look outwardly to be same. Washout process in a fuse-plug is of two 
dimensional erosion where erosion occurs in vertical and longitudinal direction freely, erosion 
is restricted in lateral direction. This causes the flow characteristics to change. 

The washout process involved in fuse plug may be taken as a preliminary study towards 
the understanding of dam breach process behaviour. 

The present work, essentially a laboratory investigation, has been conducted in two parts, 

I and 11. 

Part I presents the mechanism associated with breach erosion process of homogeneous 
earthen embankment dam caused due to overtopping flow. From the results of the present 
investigation, a methodology has been suggested to analyse dam breach process behaviour. 



Part II comprises of the experimental investigations of wash out process of fuse-plug in a 
dam body. A procedure has been suggested from the analysis of present experimental results 
to estimate the major parameters in the washout process. 

The common variables involved in the part I and II studies can be grouped as; 
i. predictor variables and ii. criterion variables. The predictor variables controlling the dam 
breach problem are the dam geometry, capacity of the reservoir and the inflow discharge to 
the reservoir. The criterion variables are the variables for which a predicted value is necessary 
in design. In the present problem criterion variables can be listed as the breach discharge, 
sediment discharge, duration of breach, time of occurrence of peak sediment, and water 
discharges and the terminal breach morphological variables like final breach width, final 
depth of erosion, final breach cross section. To relate the measured criterion variables to 
predictor variables in functional relations, the empirical constants are involved. These 
constants are called either parameters or coefficients. In the present problem of dam breach, 
the parameters / coefficients comes out of the properties of the dam fill material. These are 
cohesion (c), angle of friction {f), bulk density (y^) and average grain size {d^u) of the dam fill 
material. 

With this in view experiments were carried out in Hydraulics laboratory of the Civil 
engineering Department of Indian Institute of Technology, Kanpur for both dam breach 
process and fuse plug wash out study. A 20m long tilting flume,0.60m wide and of 0.50m 
height was used to conduct these experiments for both part I and part II studies. 

In dam breach study, 17 number of experiments were carried out with varying model 
geometry and inflow to the reservoir. Out of these 17 number of experiments, 12 sets of 
models were constructed with fine cohesionless soil (type-A soil) and 5 sets of models were 
prepared with cohesive soil (type-B soil). The dam height were varied from 14.5 cm to 35cm 
keeping the upstream and downstream slope constant as 1V;2.5H for all the / /sets of models. 
The width of dam crest was varied from 15cm to 156 cm for the dam models with 
cohesionless soil, whereas it was kept constant as 15crn for all the five dam models with 
type-B soil. The angle of friction {(p ), cohesion (c), average biilk density of the soil (y J for 
the soil fill embankment dam model of before overtopping of flow and the average grain size 
diameter((ij^) for these soils are given in the Table 0.1. 



TABLE 0.1 


DESCRIPTION OF THE PROPERTIES OF SOIL USED IN DAM 
BREACH STUDIES AND FUSE PLUG WASHOUT 


Types of soil 

^50 

(mm) 

'!> 

(degree) 

C 

(KPA) 

Yb 

gm/cc 

w 

% 

Remark 

A 

0.16 

32.1 

0.00 

1.58 

12 

Cohesionless fine sand 
collected from River Ganga 
near Kanpur. India 

B 

0.28 

33.1 

1.60 

1.98 

14 

Cohesive soil collected from 

I.LT, Kanpur campus 

C 

0.18 

29.2 

0.85 

2.03 

15 

Mix of Type-A and Type-B soil 
in the proportion 3:2 

D 

0.60 

36.6 

0.00 

1.62 

12 

Coarse sand from River Yamuna 
at Kalpi 100km south of Kanpur 

E 

0.37 

35.2 

0.00 

1.65 

12 

Mix of Type-A and Type-D soil 
in the proportion 1 : 1 


Yb 

w = 

bulk density of soil (used in dam or fuse plug) before 
overtopping 

water content of soil in percentage 


The range of constant incoming flow(Q,„) to the reservoir for dam model with type-A soil 
and type-B soil was 169.7 - fdidxlO"* mVsec and 940-6088.6 y.lO'^vcvIztc respectively. For 
all these 1 7 sets of dam models, the upstream edge of the top dam crest was kept at 8m from 
the inlet section of the flume. The measurement of the water level, sediment bed level and 
breach width along the dam cross section with time were carried out at six and eight selected 
stations, respectively for dam model with cohesionless and cohesive soils. These 
measurements were carried out by different persons using point gauges with common 
reference of time shown on an electronic digital clock. 

In part II study, on fuse plug wash out, a model of earthen dam was made in wood with a 
rectangular breach section at the centre. The dam section had the dimensions as top width(B,0 
0.15m, bottom width (Bf^) 0.90m and height (Hq) 0.145m with upstream and down stream 
slope as 1 vertical to 2.5horizontal. The central rectangular breach section was 20cm width. In 
the present study the filled portion in the breach section of this model with compacted soil is 
called as fuse plug. Rajpal (1995) conducted sixteen number of experiments on the same fuse 
plug model with type-A , type-B and type-C soil. These data are also included in the present 



analysis. Here type-C soil is the mix of type-A and type-B in ratio of 3:2. In this part 11 
experimental study two phases of experiments were conducted namely, phase-A and phase-B. 
The purpose of experiments of phase-B are discussed subsequently. The upstream edge of the 
dam was kept for all the experiments in phase-A and phase-B at 15m from the inlet section of 
the flume. The inflow variation was of the range of 560.85-1267.8 xlO~^mVsec and 247.2- 
10838.3 X 10"" mVsec for phase-A and phase-B studies respectively. In phase-B study coarse 
cohesionless sand (Type-D) and mix of type-A andtype-D soil in 1:1 ratio (type-E), are used. 
In experiments on fuse plug the water levels were measured at 4 selected stations. The 
sediment bed profiles were marked on the two non-erodible sides of the fuse-plug at different 
intervals of time. This marking was carried out using coloured pencil. At an instant five 
persons made the observations at various locations. 


Part I : Mechanism of Breach Process of Homogeneous Earthfill Dam due 
to Overtopping Flow 

The mechanism associated with dam breach process can broadly be categorised into three viz. 
a) Breach moqjhology b) Breach hydraulics and c) Breach erosion. These three aspects of 
dam breach are implicitly related to each other. 

The breach morphology is analysed from the measured sediment level and the 
corresponding breach width at different selected sections. Depletion of the reservoir level 
during the breach process is continuously measured. From the measured longitudinal 
sediment bed profile and the eroded breach section, the volume of eroded sediment is 
computed. All criterion variables relating to breach flow and sediment discharge are computed 
from the three measured basic variables viz. water level, sediment bed level, and width of 
breach at different breach sections. Using some of the computed variables, along with known 
predictor variables like properties of dam fill and embankment geometry, various non- 
dimensional breach parameters are developed like Qpf^, Qap^, Qpsn^ CAPf^, t*c which are 
very often used in the analysis. The flow over the erodible sediment crest is considered to be 
broad crested weir type of flow. The coefficient of discharge for the flow over the crest is 
analysed in terms of ratio of head acting over the spillway to the crest height. The breach 



geometry, breach out flow hydrograph and breach sediment discharge graphs are analysed 
using measured flow information. 

The analyses of data are carried out in four phases: 

a. for the model data of cohesionless soil, 

b. for the model data of cohesive soil, 

c. for the combined data of cohesive and cohesionless dam models, and 

d. for the field and lab data analysis. 

The results of these analyses are compared with the observed values of two historical dam 
breach cases namely: 

i. Teton Dam, Idaho, U. S. A 

ii. Huaccoto Dam, Peru 

The validated results include the following: 

1 . Breach peak outflow discharge(^,yp). 

2. Time of occurrence of peak outflow discharge(r^^). 

3. Final breach width, depth of erosion, slope of the final breach cross 

section at crest, At respectively. 

4. Duration of breach(t 5 ) 5 ). 

5. Breach outflow discharge hydrograph (/ ~ 

The simulated results of above breach variables for Teton dam and Huaccoto dam failure 
obtained from numerical model, BEED (Breach Erosion of Earthfill Dams) developed by V. 
P. Singh and P. D. Scarlatos for U. S. army Corps of Engineers (1989) are also compared with 
the present estimated results. 

As there is no observed information on breach sediment discharge either for Teton Dam 
or Huaccoto dam, the simulated result of BEED model is taken for comparison with the 
present estimated values. This has been done only for Teton Dam case. These validated results 
are: 

1. peak sediment discharge (05^), 

2. time of occurrence of peak sediment discharge (t^^) and 

3. sediment discharge during at any time of breach(05). 



Finally a methodology has been suggested to study the dam breach process for any given 

situation. 

Observations on Dam Breach Process 

From the experimental study on dam breach the following critical observations are made. 

a. Planform of breach channel has contraction, throat and diverging sections. The position of 
throat moves upstream as the breach process continues. The cross section of breach during 
erosion will be almost trapezoidal in shape in the contraction section and rectangular at the 
throat and again trapezoidal in diverging section. 

b. The discharge over the breach crest can be simulated with the weir of credible crest or 
spillway type of flow. The coefficient of discharge in the use of weir type formula is found 
to vary with ratio of head acting over the weir to the crest height of sediment bed profile. 
The magnitude of discharge coefficient found to decrease initially and later on remains 
constant. 

c. The rate of lateral and vertical erosion are high during initial phase of breach process and 
decrease as final state of breach occurs. The rate of lateral erosion is higher than the rate of 
vertical erosion in cohesive soil embankment and is reverse in the case of cohesionless soil 
embankment breach. 

d. The depth of erosion, width of erosion and lateral slope of bank at the crest vary with drop 
in water level between the upstream and downstream of the embankment during breach 
process. The final breach width varies with magnitude of peak flow discharge. The final 
depth of erosion varies mainly with height of the dam. 

e. The mean velocity of flow in a breach is a function of average hydraulic radius and 

average energy slope( Sf ) and varies as: S°'^* . The flow velocity at downstream top 

edge of the embankment varies with critical velocity of flow and water surface slope. 

f Sediment discharge intensity q, is found to vary with water discharge intensity with a 
power as and sediment load per unit width varies as the square of average bed shear 
stress for dam model with cohesionless soil. But in a generalised analysis for any type of 



soil fill it varies in a quadratic form with the square of average shear velocity along with 
Qpsn and numbers. 

g. Relative energy loss in the breach is found to vary with average water surface slope, square 
root of Qpsisi number along with average shear velocity. 

h. The final breach geometry like breach width(5c^)y, depth of erosion(Dg^) and flow 
parameters like peak water discharge(2(f'/’) and sediment peak discharge(^ 5 p) their 
time of occurrence t^p and t^p respectively and time of breach % are analysed using lab 
data and field data. The functional relations are developed based on the geometry, soil 
properties of embankment, reservoir capacity and inflow characteristics into the reservoir. 
Two field dam broach failure data are included in the analysis. 

i. The variation of outflow hydrograph and sediment discharge graph based on the shape of 
hydrograph is presented. The validation of these functional relations are strengthened by 
including two field breach data of Teton dam and Huaccoto dam in the analysis. 

j. A method in the form of algorithm and a flow chart is presented for computation of 
variation of breach geometry, flow discharge, sediment discharge during the breach process 
and at the final stage of breach. 

The details of breach process evolved from the experimental study will be of helpful in 
getting a good insight into the mechanism of breaching and in the development of a better 
numerical model for dam breach. 


Part II : Study on Washout Process of Fuse Plug 

The experiments were conducted in two phases for Fuse-plug washout study. 

Phase-A This phase of the study refers to the experiments which were conducted to know 
the mechanism of wash out of the fuse-plug. The major objectives of this part of study were to 
obtain the following viz. a) water level variations in the reservoir, b) water level variation 
over the fuse-plug model and c) sediment bed variation with respect to time. With these three 
basic information, the involved mechanism in soil-water interaction under highly dynamic 
condition like washout process of fuse-plug are analysed. In this phase, six number of 
experiments (with cohesionless fuse-plug fill) were carried out. In the analysis, Rajpal 



(1995)’s data obtained from sixteen numbers of experiment on this fuse-plug model with 
cohesive and cohesionless soil fill are used. The inflow discharge were varied keeping the 
fuse-plug geometry constant. 

Phase-B The second phase of experiments were carried out on rigid bed profile analogous 
to the instant erodible bed profile of fuse-plug wash out. For this purpose six erodible profiles 
of run number 1404 were chosen and experiments were carried out with different 
condition by varying inflow discharge to the reservoir. Here, is head acting over the crest 
and is crest height of washout profile. In this category altogether 35 experiments were 
carried out. The main objective of these experiments were: 

a. to verify outflow discharge relation for a erodible bed of constant width 

b. energy loss characteristics during wash out of the fuse-plug. 

In an effort to know the effect of particle size on energy loss characteristics in wash out 
process, four number of experiments were conducted with coarse((i 5 ( 5 = 0 . 6 mm) cohesionless 
soil. 

Analysis of Fuse Plug Washout Process 

From the rate of depletion of reservoir water level, two scales viz. reservoir time scale and 
reservoir drop scale were obtained. These two scales and other non-dimensional parameters 
like A",,, Qpm 2 extensively used in the analysis to know the rate of sediment crest erosion 
and sediment discharge characteristics. Finally, a procedure is suggested for the prediction of 
wash put process of fuse-plug. The sediment discharge (^ 5 ) and water discharge (Qfy) can be 
computed at any time using the suggested methodology in addition to other erosional features 
of the fuse plug washout. 

Observations on Fuse Plug Washout 

1 . The flow hydraulics over the crest of erodible bed profile does not follow weir formula, 
rather the outflow discharge depends on head of water(;z„) over the crest and its height(Ac„)- 

2. Reservoir drop scale (/f.)and time scale (T.) are two most important parameters controlling 
the washout behaviour of the fuse plug. 

3. Rate of crest height erosion (dh^Jdt) with time depends on incoming flow velocity ^2g\ 
and its maximum value {dh^Jdt)^,. varies with , where q^p is peak outflow discharge 
intensity. 



4. Sediment concentration {C=q/qJ) is almost independent of average bed shear velocity V, 
However, it is found to be a function of drop in water level between upstream and 
downstream of dam AH. 

5. Sediment concentration during the washout of fuse plug and the ratio of energy in the 
downstream {Ei) to upstream {Ei) of fuse plug, (E/E,) are functions of 

6. From the experiment in phase-B the discharge relation for flow over nonerodible fuse-plug 
crest follows the weir type relation with discharge varying with head acting over the crest 

as . This observation confirms the variable nature of discharge coefficient q in 
washout process is due to the crest erosion. 

7. Explicitly energy loss due to erosion and transport of sediment is related with sediment 
concentration(C) for type-A soil. 

8. The duration of washout(%) is found to be varying between 1 to 2.5 times of T* depending 
on inflow discharge, soil resisting velocity and height of the fuse plug. 

9. The ratio of time of occurrence of sediment peak discharge (t^p) to time of occuirence of 
peak water discharge vary inversely to the square root of peak water discharge 
intensity. 

10. Peak water discharge intensity (q^^p) is proportional to nearly square root of inflow 
discharge intensity( ) 

The results of the present study have a direct bearing on field applications. An insight into 
the basic process involved in the dam breach and washout of fuse plug has helped in 
suggesting the methodology to predict the breach geometry, water discharge and sediment 
discharge during breach erosion. It is hoped that the insight knowledge obtained through this 
experimental study will help in better understanding of the flow mechanism involved in the 
breach process. The evolved functional relations among the breach variables will enable one 
to formulate the numerical model on dam breach and fuse plug wash out in better scientific 
and more logical way. 
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1 . Crest ; It refers to the position where erodible sediment bed profile height is maximum 
in breach section 

2 . Upstream and downstream top edge : It refers to the upstream and downstream end of 
the dam top length 

3 . Variables : Without any prefix, it refers to the controlling input and output variables in 
general involved in the process. 

Basic Variables : In the present study three basic measured variables refer to variation of 
water level, sediment bed level and width along the breach. 

Computed Variables : It refers to the computed from the measured basic variables like water 
and sediment discharge, average shear stress etc. 

Predictor Variables : Predictor variables are variables for which values are measurable 

Criterion Variables : It is the variables for which a predicted value is necessary for its 
estimation. 

4 . Parameter : A numerical measure of a property or characteristic that is constant for a 
system under specified conditions. 

5 . Coefficient : A variable or constant appearing in a mathematical relation, each value of 
which defines the specific form of the functional relation. 

6. Scales : It refers to some critical value of a particular variable, which are used to in the 
analysis like X50, 150, tjp etc. used in the present analysis. 
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7. Erosion Coefficient : In the present method it is defined as the ratio of rate of crest height 
variation to incoming velocity. 

8. Erosion Velocity : It refers to the average velocity of eroded soil from the fuse plug body. 

9. Drop : It refers to the drop in water level either at crest or between reservoir level and 
downstream of dam (AH) 

\Q. Analogous Rigid Bed: It refers to the rigid bed profile modelled in replica to the 
instantaneous mobile sediment bed during the washout process of fuse plug. 

1 1 .Soil resisting Velocity. This indirectly account for the resistance to erosion 

12. Embankment Dam : The homogeneous dam constructed without filter or core inside the 

dam body. 

\3. Erosion Process : It refers to the internal mechanism between flow parameter, breach 
morphology and soil properties of the dam fill materials. 

lA.Washout : Breaching process in connection to the fuse plug (through its width uniformly) 
this word is used as the fuse plug is meant to be washed away during high flood. 

\5. Normalised : It refers to the nondimensional form of variables by its corresponding 
maximum value. 



CHAPTER I 


INTRODUCTION 


1.1 An Overview of Dam Breach 

Among the natural resources, water is vital to life and economic growth. Dams play a major 
role in harnessing the surface water component of this precious resource by storing and 
distributing it in a proper planning and management for various economical use. Earthen dams 
are the most common among various types of dams constructed. It constitutes a vast majority 
of dams compared to other categories. Failure of a dam in any form is not desirable. When a 
dam fails, it results in the sudden release of stored water, which becomes a disaster in the 
downstream region. Among the dams, the earthen dams are more vulnerable to failure. The 
major causes of failure of earthen dams are overtopping and piping. Overtopping of flood on 
an earthen dam causes a breach in the body of the dam. The breaching phenomenon is a short 
time event due to which all parameters controlling the breach erosion are not possible to 
measure in field. However, under laboratory conditions, the phenomenon can be simulated to 
a good extent, which gives a better insight into the mechanism involved in breach process as 
shown in Figure 1.1. A better knowledge will help to develop more realistic numerical model 
for an earthen dam breach. 





Fig. 1.1: Simulation of dam breach phenomena due to overtopping flow in laboratory 
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1.2 Concept of a Fuse Plug 

Almost all dams should have spillway to permit excess flood discharge downstream to 
prevent the dam to be damaged. Gates on spillways permit the controlled operation of flood 
discharge. In the design of a dam, it is found that a weak granular embankment, which may or 
may not have a clay cover, is built at a suitable location of topography as indicated in Figure 

1.2 as a possible location of fuse plug to act as automatic excess flood releasing device. 
During high floods these auxiliary outlets which augment the service spillway by passing the 
larger flood flows are called fuse plugs. 

1.3 Types of Dam Failures 

The most common causes and modes of failure of earthen dams can be attributed to: 

a. Overtopping caused by extreme floods. 

b. Structural failure due to internal erosion (piping). 

c. Structural failure due to shear slide. 

d. Structural failure due to foundation problem. 

e. Failure due to natural or induced seismicity. 

It is found that more than 60 percent of earthen dam failures are caused by overtopping or 
piping (USCOLD, 1975; Middlebrooks, 1953). 

1.4 Causes of Overtopping of Flood 

The overtopping of flood over the embankment may be due to following reasons: 

a. Design criteria of probable maximum flood (PMF) for the reservoir 

i. Ignoring the future catchment yield due to hydro-geological changes. 

ii. Reduction of live storage capacity of reservoir due to sedimentation. 

b. Reservoir operational rules 

i. Spillway and sluice gate operation during flood. 

ii. Seasonal reservoir operational policy. 




FIGi 1.2 SCHEMATIC DIAGRAM FDR POSSIBLE LOCATION OF FUSE PLUG 
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1,5 Necessity for Dam Breach Study 

According to Gruner (1967a, b), the total number of dams that impose risk of serious damage 
in case of failure may well exceed 150,000. Jansen (1980) has noted that there have been 
approximately 2,000 dam failures around the world since the thirteenth century. Out of these, 
about 10 percent have occurred in 20th century itself Berga (1992) has pointed out that 41 
percent of historic dam failures were caused by the design flood having been exceeded. 
Loukola et. al. (1993) noted that 80% of nearly 87,000 dams in China are embankment dams 
and more than 90% are small dams. The dam failures in China associated with earthen dams 
are of the order of 93% and 96% with small earthen dams or small reservoirs. The main cause 
of failure was overtopping due to heavy rainfall and ensuing floods low design criteria for 
small dams and operational problem in appurtenant works. In India there are enough examples 
of earthen dam breach in recent past. A few of them can be cited, namely Machhu II 
(Irrigation scheme) dam (Gujarat), Panshet dam, Khadakwasla dam (Poona) Maharashtra, and 
Kodaganar dam (Tamilnadu) which breached either due to inadequate spillway capacity or 
due to inadequate provision of the outlet facility during emergency. 

It is easy to infer from the statistical figures of the probable dam failures, that it is 
essential to have contingency planning for evacuation during the dam breach with the help of 
reliable flood forecasting in the downstream flood plain. In some of the dam-break numerical 
modelling, the routing of flood and sediment downstream are mostly carried out by assuming 
the failure of the dam as an instantaneous phenomenon. But from historical earthen dam 
failure records it is evident that the duration of failure varies from about half an hour to 12 
hours depending upon the situations. That clearly indicates dam failure is a time-dependent 
and not a instantaneous process as assumed. Therefore, in order to accurately predict water 
and sediment discharge from the breach section of the dam, it is necessary to know as to how 
the breach process evolves with time. These quantities are controlled by the breach cross 
section or the morphology of the breach with time. The dam breach flood and its associated 
sediment discharge can be routed in downstream channel in a more accurate way to predict 
the stage and discharge at any reach, provided the breach fonnation in the dam can be 
predicted with time for a given condition. 
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1.6 Variables and Parameters Involved in the Problem 

In the dam breach study the two types of variables involved are: 

a. Predictor variables 

b. Criterion variables 

The predictor variables controlling the present problem are: 

i. Dam geometry 

ii. Capacity of the reservoir 

• Capacity and elevation relation, 

• Variation of discharge characteristics of spillway, sluice gate and other outlets, with 
head of water 

iii. Inflow discharge to the reservoir. 

The criterion variable is the variable for which a predicted value is necessary in design. In the 
present problem criterion variables can be listed as: 

a. Peak breach outflow discharge. 

b. Peak outflow sediment discharge. 

c. Times to peak water and sediment discharges and duration of breach. 

d. Terminal breach morphology. 

i. final breach width, 

ii. final depth of erosion, 

iii. final side slope of the breach cross-section. 

Some parameters or coefficients are needed to relate various criterion variables to predictor 
variables in a functional relation. Parameter, a numerical measure of property or 
characteristics that is constant for a system under specified conditions. In the present problem 
the parameters / coefficients are obtained firom the properties of dam fill material, namely, 
cohesion (c), angle of fiiction {(^, bulk density of the dam fill material and average grain 
size (Jjo). 
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1.7 Present Status of Knowledge on the Dam Breach 

Conceptually, the dam breach hydraulics and the breach morphology due to overtopping flood 
comprises the dam breach process in a broader sense. Enlargement of the breach affects the 
rate of discharge, which subsequently controls the rate of erosion. This process continues till 
the reservoir water is depleted or till the dam resists further erosion. 

Briefly the present state of the knowledge relevant to the dam breach problem may be 
described under the following subheads. 

a. Erodibility 

b. Mechanics of breach 

c. Hydraulics of breach flow 

d. Erosion of sediment and its transport 

a. Erodibility 

Hydraulic erosion entails a soil-water interaction. Erosion occurs when shearing stresses 
induced by fluid flow on soil surface are high enough to overcome the force holding the 
particles together and to cause the dislodgment of particles from the surface. A cohesive soil 
erodes in a quite different way as compared to a cohesionless soil. In a cohesionless soil, 
erodibility primarily depends on grain size distribution, density and shape of the grain while 
the main resistance to erosion is offered by the submerged weight of the sediment. In a 
cohesive soil, the amount and type of clay percentage, chemical composition of the pore fluid, 
degree of saturation, presence of organic material and other cementing agents play a 
significant role in its erodibility. 

b. Mechanics of Breach Formation 

The basic mechanism of breaching of earth fill dams is rooted in the erosion of the 
embankment material. Excessive shear stress on the surface induced by the flow of water over 
weak section of dam body initiates the process of erosion. The breach in the initial stages is 
predominantly ‘V’ shaped, three to four times as wider than its depth. As the erosion of 
surface progresses laterally the breach section grows. The extent of this lateral erosion 
depends upon the size of the reservoir and volume of water stored in it. 
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c. Hydraulics of Breach Flow 

Most of the earlier experimental and numerical studies assumed: 

i. Flow over the crest of the breach is analogous to the flow over the broad crested weir. 

ii. The unsteady flow is approximated by quasi-steady situation. 

iii. The tail water effect is considered suitably by reducing the discharge coefficient in some 
of the numerical models. 

iv. Manning’s or Chezy’s resistance equation used for flow depth calculation in the breach 
channel. 

Pugh and Gray (1984) professed that the flow over the entire top section of the breach 
can be assumed to be critical i.e. initially an S 2 water surface profile and supercritical flow 
thereafter along the breach. The governing equations for the flow over a dam are as: 

a. The continuity equation. 

b. Broad-crested weir formula. 

c. Breach shape geometry relation. 

d. The equation of resistance to motion. 

e. The equation of motion. 

d. Breach Morphology 

Initial cross-section of the breach is of a narrow valley of V-shape. As the erosion progresses, 
the shape modifies to trapezoidal cross-section. In the numerical models on dam breach, the 
shape of breach section has been assumed to be of various shapes like rectangular, triangular, 
trapezoidal, parabolic and semicircular as reported in earlier works on dam breaching. But 
from the experimental observation it is found, that the terminal cross-section of the breach 
may be taken to be trapezoidal. There are two cross-sectional shapes derived on theoretical 
considerations based on: 

i. the concept of stable section of greater efficiency as a simple cosine curve, and 

ii. the minim um cross-sectional area maintaining stability for a given discharge. Lane 


(1985). 



e. Erosion of Sediment and its Transport 

Due to lack of experimental evidence and not having a complete understanding of sediment 
transport relations under highly dynamic conditions like a dam breach, most of the studies 
reported so far have been forced to resort to sediment transport formulae applicable to alluvial 
streams, like Du-Boys formula(1879), Schoklitsch formula(1935), Meyer-Peter-Muller 
foirnula (1948) and Einstein-Brown formula (1950). These formulae are based on the tractive 
force theory i.e., sediment bed load is a function of average shear stress, and are applicable to 
homogeneous bed load in alluvial streams. 

The governing equations in breach erosion and sediment transport are: 

a. Continuity equation of sediment. 

b. Equation of motion of sediment. 

c. Shear stress relation with bed load. 

f. Stability of Breach Slope 

The side slope of the breach section becomes unstable during breaching when the gravity and 
hydrodynamic forces due to seepage become greater than the soil friction and cohesion and 
collapses, causing enlargement of breach section. Greater the breach flow, greater will be the 
breach erosion. 

1.8 Aims of the Present Study 

a. The present state of knowledge on earthen dam failures due to overtopping flood is 
postulated, mostly by assuming the principles of sediment transport in the rivers and 
hydraulics of flow over a weir. An insight into the mechanism involved in the breach 
erosion process is very essential for the better understanding of this type of earthen dam 
failures. With this aim, experiments are planned to study the breach process on 
homogeneous embankment dam models. Knowledge, thus obtained, is validated using field 
information. 

b. Fuse plug is a safety-valve for the dams. Little information is available on the washout 
process of the fuse-plug due to overtopping of flood. Therefore systematic experiments are 
planned to get a know-how on the washout mechanism involved in the erosion process of 
a fuse plug. 
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1.9 Organisation of Thesis 

Experimental work has been carried out in two phases and the thesis is being presented in two 
parts I and 11 


Chapter I 
Chapter II 
Part -I 

Chapter III 
Chapter IV 
Chapter V 
Part-II 
Chapter VI 
Chapter VII 
Chapter VIII 


Introduction 
Review of Literature 

Breach Erosion of Homogeneous Earthfill Dams Due to 
Overtopping Flow 

Scheme of Experimentation on Dam Breach 
Mechanism of Dam Breach Process 

Comparison of Dam Breach Results with Teton and Huaccoto Dams 

Washout Process of Fuse Plug 

Scheme of Experiments on Fuse Plug Washout 

Analysis and Results of Washout Process of Fuse Plug 

Conclusions and Recommendations 


This is followed in succession by References and an Appendix. 



CHAPTER II 


REVIEW OF LITERATURE 


2.1 Introduction 

A considerable body of literature has been accumulated over a span of last 150 years on the 
movement of flood wave resulting from instantaneous dam failure which in academic circle is 
coined as ‘Dam Break’ problem. But the historical earthen dam failure record derived from 
studies by Middlebrook:s(1953), Gruner(1969), Babb and Mennel(1968), American society of 
Civil Engineers(1975), Singh and Scarlatos(1989), Singh and Snorrason(1982), Ponce(1981), 
Central Water Commission, Govt, of India(1989), and Jansen(1988) clearly have shown one 
striking common denominator: the failure is anything but sudden. So the gradual failure of 
an earth dam is of particular interest to disaster relief planners, because the rate of growth of 
the breach strongly influences the peak and shape of the ensuing flood waves. As a result the 
significant interest on dam safety among public has come in the last two decades or so. The 
research on different facets of earth dam breach has been augmented through different 
empirical and numerical modelling without much experimental data support. Going through 
the relevant literature on dam breach it seems to be appropriate to present the different 
concepts and methodologies of prediction of breach characteristics in the following subheads: 
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a. Experimental Observations. 

b. Empirical Modelling. 

c. Mathematical Modelling. 

d. Ancillary literature relating to dam breach. 

2.2 Experimental Observations 

Powledge et al (1989a, b) carried out model and prototype research studies which have been 
conducted in the United State and Great Britain to evaluate how embankments for dams, 
levees, roadways perform when subjected to overtopping flows from probable maximum 
flood (PMF) or near PMF events. From their experimental observations, they concluded flow 
over an embankment passes through three zones, sub-critical flow over dam crest, 
supercritical flow over dam crest and supercritical flow on dam downstream slope. They have 
identified a potential negative pressure zone on the downstream slope immediately below the 
crest. The potential negative pressures can be quantified by superimposing the Nappe equation 
(theoretical ogee shape) on the embankment crest and slope profile from which sub- 
hydrostatic zone can be located. 

Powledge et al disclaim about the applicability of the tractive stress approach for erosion 
rate analysis during the breaching process. The initiation of breach starts near the slope 
discontinuity, such as the toe or base of the dam. They have observed that the overflow 
discharge may be computed by using weir flow equation for both granular and cohesive 
embankments. In granular embankment, once the overtopping occurs, the seepage on the 
downstream slope accelerates slope erosion. The erosion phenomenon starts at or near the toe, 
where seepage is exiting and the turbulence of the impinging flow is greatest. Then quickly 
the erosion process propagate rapidly up the slope and cause the lowering of the crest, 
allowing increase in flow rates which accelerate the erosion process. For cohesive 
embankment, generally there is more resistant to erosion than granular fill due to existence of 
low seepage on the slopes because of low permeability. Factors which influence the erosion 
process are; 
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i. Embankment configuration 

ii. Maximum velocity attained by flow 

iii. Discontinuities: such as bemis, outlets, creates, local changes in momentum producing 
turbulence and erosion 

iv. Presence of tail water 

Their main aim of the experimental study was to evaluate performance potential of 
different protective measures to prevent breach erosion. Focus on the mechanism of breach 
process was not given. 

In August’ 1995 the sea defense at a MAFF / English Natural Experimental site at 
Tollesbury in Essex was deliberately breached. This prototype breaching was a part of a series 
of experiments into an approach to flood defence known as ‘managed set-back’. The 
breaching operation, carried out as part of the MAFF flood and Coastal Defence Division 
strategic research program was carefully controlled and monitored by HR Wallingford to 
study the behaviour of clay embankments under simulated flood conditions. Information’s 
about breach details are not available to the author. 

The report from Streamline {December’ 1995) can be cited here in connection to the dam 
breaching: 

[There are number of coastal lagoons in urban zones in New South Wales, Australia. 

These are usually blocked by a sand barrier built up by wave activity and reach 2- 
3m above sea level. Heavy rainfall causes water level in the closed lagoons to rise 
rapidly putting urban area at risk of flooding unless an entrance to the ocean is 
scoured through the barrier. This process can be initiated naturally by overtopping, 
or digging a small “starter channel” across the berm.] 

HR Wallingford has studied above problem at the Wamberal Lagoon in New South Wales. 
Following are observed from the prototype non-cohesive embankment breach model test. 

i. initial vertical cutback across the breach 

ii. fomration of a curved weir upstream in the sand plug behind the beach berm 

iii. rapid widening and deepening of the gap in the berm and final sub-critical phase 
These two experimental studies gave only qualitative information on embankment breach. 
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2.3 Field Observations 

From the 52 historical dam failure cases Singh and Scarlatos (1988) studied inter-relationship 
among terminal depth of erosion, top and bottom width of the breach section and side slope 
of breach section shown below: 


SI No 

Relation 

Qualitative values 

1 

B/b 

= 1.06 -- 1.74, with mean value of 1 .29 and standard deviation 0.180 

2 

B/d 

= 0.84 - 10.93, with mean value of 4.18 and standard deviation of 2.62 

3 

tan 9 

= 1.0 

4 

tf 

1/2 hr to 12 hr 


Where B = final top breach width, 

b = final bottom breach width, 
d = final depth of erosion, 

0 = angle of side slope of the breach section with vertical, and 
t)- = failure time. 

The qualitative relations bring an insight into the problem of dam breach. 

2.4 Empirical Model 

The models express formalized concepts of real-life problem. They consists of two elements : 
input variables and empirical constants. Again the input variables can be categorised into 
two, namely: 

a. Criterion variable 

b. Predictor variables 

The criterion variable is the variable for which a predicted variable is necessary in design. 
Predictor variables are used to make estimates of the criterion variable. Like in the present 
topic of interest, capacity of the reservoir, inflow hydrograph are examples of predictor 
variables, whereas the peak breach discharge, final breach width, time to peak discharge, 
duration of breach, final depth of erosion are examples of criterion variables. The empirical 
constants are derived by fitting the model with measured values of the criterion and predictor 
variables. These constants are either called coefficients or parameters. So parameter can be 
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defined as a numerical measure of a property or characteristic that is constant for a system 
under specified conditions, whereas coefficients is a variable or constant appearing in a 
model, each value of which defines the specific form of the model. 

Major contributions to the modelling of breach process of an earthen embankment dam is 
from Singh and his associates. They looked in to the breach processes in depth and developed 
both empirical and mathematical models. Their contribution can be considered as a milestone 
in the modelling of embankment breach process. Singh and Scarlatos (1988) have developed 
an empirical model by relating two basic governing equations in dam breach problem i.e., 
reservoir water balance equation and rate of erosion. The water balance equation is: 

(2.1) 

at 

Where = water surface area within the reservoir, /= inflow to the reservoir, Q= flow 

due to crest overtopping, <2* = breach flow, (dH/dt) = rate of change of water level in the 
reservoir. They have brought down the equation 2.1 to a simple form by equating rate of 
change of volume in the reservoir to mean breach flow velocity (u) times the wetted breach 
area (Ai,) assuming inflow to the reservoir (/) and dam crest overtopping flow (Q) to be 
negligible. The average breach flow velocity (u) is calculated by adopting hydraulics of broad 
crested weir flow (Chow 1959, Pugh and Gray 1984) as: 

u = a,{H-Zf' (2.2) 

Where Z= the bottom elevation from reference datum, and H= the water level surface 
from a reference datum. For critical flow conditions, the empirical coefficients: 

2 

? 

A =0.5. 



SO 


u = 1.7 - Z) ml s . 
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With the assumption erosion is directly proportional to shear and subsequently proportional to 
the water velocity. Then equation (2.1) reduces to: 

^=-a,u{H.zf (2.3) 

Value of P 2 is related with the concept of Laursen, the rate of sediment transport is a power 
function of mean velocity. The values of P 2 vary in the range 4 to 6. If the value of P 2 is an 
integer, the solution to equation (2.3) becomes a closed form one. 

They adopted three different breach shapes namely, rectangular breach cross-section with 
constant width, which enlarges in vertical direction, triangular breach with constant side slope 
and enlarges similar to rectangular breach, trapezoidal breach cross-section with constant 
bottom width and constant side slopes. With known initial conditions, H=Ho, and Z=Zo at t=to 
and stated assumptions, closed form solutions are developed for each breach cross-section. 
Depending on the value obtained for [32, two different cases, i.e., linear erosion (P 2 = TO) and 
non-linear erosion (|32 ^ 1.0) are studied. In all these cases it has been tried to determine 
reservoir water level and breach bottom elevation with time i.e., H(t) and Z(t). The erosivity 
coefficient a 2 was estimated through calibration, based on maximum outflow discharge and 
the failure time. For 16 historical dam failure case 02 is computed and is reported, it varies 
within a range 0.0008 ~ 0.008 for linear case and 0.0001 5 ~ 0.0021 for non-linear situation. It 
is found that the overall performance of the linear rate of erosion is better than the non-linear 
erosion rate. However the nonlinear erosion models gave better results for the rising limb of 
the breach hydrograph, while the linear models simulated better the recession limb. 
Calibration of erosivity coefficient in their model can be considered as major drawback for the 
application of their model. This could be overcome only when experimental data is available 
for prediction of erosivity coefficient. Validation of erosion coefficient with experimental data 
is needed for the application of their model to the field problems. 

Macchione(1989) derived an analytical solution to the work of Singh and Quiroga(1988). 
He expressed the reservoir storage volume as a function of the depth of water (H) in the 
reservoir. He proposed the exponent to H, to vary from 1 to 4. Regarding breach geometry he 
stressed that the final breach shape would be trapezoidal. Initially a triangular breach is 
formed, when it reaches the natural ground, the breach expands laterally. He stressed that the 
breach flow can be taken as critical, without undue loss-of accuracy. He has used du Boy’s 
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equation for movement of sediment. His analytical procedures describe the evolution of 
breach hydrograph at the beginning for triangular breach and then for trapezoidal breach. He 
has defined explicitly the expression for the reservoir water level variation with time for both 
the breach sections triangular and trapezoidal shape. There is no mention to differentiate 
triangular breach from trapezoidal breach with respect to time. .His visualisation of breach 
formation from triangular to trapezoidal shape keeping the angles of sides of breach constant 
does not satisfy the laboratory observation on the dam breach. 

A close form dimensionless analytical solution to Equations 2.1 and 2.3 is tried by Singh 
and Quiroga(1988) for rectangular breach section. In this approach they have tried to relate 
the water surface elevation(Z) in a non-dimensional fomi to non-dimensionalised H, they have 
used dam bottom elevation (Hb) and Hq, the reservoir water level at the starting of the breach, 
t = 0. Similarly the breach bottom elevation Z is non-dimensionlised with Hp and Zq, where 
at t = 0, Z(0) = Zq. From the given dam geometry Zq and H), can be found out. Then they 
have related H and Z with respect to time. The time (t) is also non-dimensionlised with 
breach formation time (tf). From this solution they have come out with following view: 

“...For a breach to develop completely, the rate of erosion at the breach must be 
larger than the rate of depletion of the reservoir. This means that peak discharge 
occur when the difference between non-dimensional reservoir level elevation and 
breach bottom elevation is maximum. This normally happens at time is equal to the 
breach formation time” 

With this approach, only the rising part of the breach hydrograph can be obtained, not the 
complete picture of hydrograph. A single overall erodibility factor is used for the rate of 
sediment erosion. 

2.5 Mathematical Modelling 

A good number of mathematical models for simulation of earth dam failure have been 
developed in the past twenty five years or so, though the models vary considerably in their 
physical content. 

These models are: a) Cristofano model (1965), b) Harris-Wagner model (1967), c) BRDAM 
model (1977, 1981), d) Ponce-Tsivoglou model (1981), e) Lou model (1981) . 



Review of Literature 1 8 


f) DAMBRK model (1984a and:1984b), g) SMPDBK model (1984), h) Breach model (1984b) 
and i) BEED model (1986a, 1987a, 1987b, 1988a). 

The various components of these models are 

i. Reservoir hydraulics 

ii. Hydraulics of flow over dams (breach hydraulics) 

iii. Breach morphology 

iv. Flow routing 

From these four angles these models can be reviewed, giving special importance to the breach 
hydraulics and breach morphology, as these two aspects are addressed in the present study. 

a. Reservoir Hydraulics 

The inflow routing and reservoir water elevation study comes under the component of dam 
breach modelling. Each of these models route the inflow using storage routing method and 
computes changes in the reservoir water level by volume balance. Outflows through 
spillways, outlet works and pdwer plants during the breach constitute a part of the storage 
routing. 


b. Breach Hydraulics 

Except PT(1981), Lou(1981) and Nogueira(1984) models, all other models use the broad- 
crested weir formula. For a rectangular broad-crested weir with b defining the width of the 
breach channel and hy the head acting over the crest, the discharge of flow over the weir Q is 
given by 


e=*Vi' 


3 ^ 






= \.lbH^ 


(2.4) 


where b is the width of breach, *// the flow depth and Q the discharge. 

The PT model employs a relation between breach top width and flow rate. This relation 
was applied from inception to peak flow, after which the breach width was kept constant and 
equal to the peak flow width. Lou and Nogueira models utilize respectively St. Venants 
equation and critical flow equation. 
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A majority of the models do not separately consider the flow through the breach channel. 
The PT, Lou and Nogueira models use the St. Venant equations for modelling this channel 
flow, whereas the BREACH and BEED models use a steady uniform formula like Manning’s 
equation. The majority of the models never treat the flow through the crest and breach 
charmel, separately rather combinedly as hydraulics of flow. 


c. Breach Geometry 

Each model characterizes the breach morphology differently. Except the analytical 
expressions given in Lou and PT model for breach geometry, other models assume a breach 
shape. The DAMBRK, SMPDAM and BEED models use a trapezoidal shape but .the 
program has the facilities to use rectangular and triangular shapes also. The Cristofano model 
uses trapezoidal breach assuming a constant bottom width. The same assumption was also 
made by the PT model. The HW model uses a parabolic shape, whereas the BRDAM model 
represents the breach by a triangular shape with 45° side slopes and parabolic bottom. The 
BREACH model expresses breach width as a function of the depth of flow in the breach. Lou 
model considered that the breach formation due to overtopping grows in space and time until 
an equilibrium section is attained. Based on this mechanism and hydraulic considerations, a 
most efficient stable section is derived for describing breach shape as: 

y = DCos[f^f^ (2.5) 


where y = depth at a distance x, D = maximum depth at the center, x - distance from the 
center, and B = top width. This shape results in the smallest width having greatest hydraulic 
radius in order to keep channel stable. Similarly, Nogueira model(1984), following the work 
of Lane(1955) and effective shear stress analysis, an equation was derived to represent the 
breach profile as : 


/ 

T = To Cos 


xXzx\(f) 


( 2 . 6 ) 


k To 2 

where, y = depth at a distance x, yQ = maximum depth at the center, x = distance from the 
center, and (p = angle of repose of the bed material. The profile is a cosine curve and has a 
minimum cross-sectional area for a given discharge while keeping the stability. 


i 
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d. Sediment Transport 

The sediment transport equations that have been used are Schoklitsch formula, the Meyer- 
Peter and Muller equation, the DuBoys equation, and the Einstein-Brown equations. Three 
models namely PT, Lou and Nogueira, coupled the transport equation with Exner sediment 
continuity equation; 


dt 




dt 


= 0 


(2.7) 


Where Q= Sediment transport rate by weight; p - porosity of the material forming the bed; 
B = average charmel width; Z = bed elevation, density of bed material. The change in bed 
elevation can thus be computed as a function of sediment transport rate and bed material 
properties. The Cristofano model derives its own erosion equation by equating the rate of 


change of vertical and lateral erosion within the notch. He assumed to the rate of change of 
water flowing through the notch for a unit area as: 


= exp 


lXdX\(f) 

h 


( 2 . 8 ) 


where = Volume of sediment discharge; Q), = Volume of water discharge through the 

breach, = Constant of proportionality, / = the base length of the flow at any time, and h = 
the hydraulic head at any given time. The Lou model has provision whereby sediment 
transport is expressed as a function of the duration of failure, erodibility index, and velocity of 

C 

water as: E = (2.9) 


where Cg = constant, = duration of failure, = velocity of water and E = mass of soil lost. 


e. Breach Slope Collapsing 

Among the mathematical models discussed only tliree models, Nogueira, BREACH and 
BEED model, account for side slope instability and collapsing in a explicit way. The 
Nogueira model considers six cases to analyse lateral erosion by sliding, based on the 
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assumption that the sliding will occur into the breach along a line with the slope equal to the 
tangent of the angle of internal friction of the material. The SEED model uses the contour 
method to analyse the slope collapsing assuming saturated soil conditions. With similar line, 
BREACH model analyses the forces to compute collapsing of the breach side assuming dry 
soil conditions. 

f. Downstream Flood Routing 

The flood routing facility in the downstream channel and flood plain is in-built in the BEED 
and DAMBRK models. Wliereas BEED model used Muskingum method, the DAMBRK 
model employs the dynamic flood routing based on the St. Venant equations. The BRDAM 
and SMPDAM models route the flow only in downstream channel. The former uses dynamic, 
wave routing, and the latter uses dimensionless curves produced by the DAMBRK model. 
Other models have not incorporated flow routing in their model. 

In the Table 2.1a and 2.1b, the different aspects like the components in various 
mathematical models and their treatment are summarised. The information about the different 
mathematical models listed in table 2.1(a), and 2.1(b) are collected from the book on ‘Dam 
Breach Modeling Technology’ by Singh(1996) and are presented in two broad headings viz. 
breach hydrulics and breach morphology. 


TABLE 2.1a 

COMPONENTS OF VARIOUS MATHEMATICAL 
MODELS FOR DAM BREACH EROSION 


Bre ach Hydrauli cs BreaiA Morphology 


Mode! 

Author (year) 

Flow over 
the Crest 

Flow through 
the Breach 

Flow through the 
Breach on the 
Downstream face 

Sediment 

Transport 

Breach Shape 

Side Slopes 
Collapsing 

Cristofano Model 

Cristofano (1965) 

No 

Yes 

No 

Yes 

Yes 

No 

HW Model 

Harris and 

Wagner (1967) 

No 

Yes 

No 

Yes 

Yes 

No 

BRDAM Model 

Brown and 

Rogers (1977) 

No 

Yes 

No 

Yes 

Yes 

No 

PT Model 

Ponce and 
Tsivoglu (1981) 

No 

Yes 

Yes 

Yes 

Yes 

No 

Lou Model 

Lou (1981) 

No 

Yes 

Yes 

Yes 

Yes 

No 

Nogueira Model 

Nogueira (1984) 

No 

Yes 

Yes 

Yes 

Yes 

Yes 

DAMBRK Model 

Fread (1984a) 

Yes 

Yes 

No 

No 

Yes 

No 

SMPDBK Model 

Wetmore and 

Fread (1984) 

Yes 

Yes 

No 

No 

Yes 

No 

BREACH Model 

Fread (1984b) 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

BEED Model 

Singh and 
Scarlatos (1989) 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 


In Table 2.1(b) the various aspects pertaining to breach hydraulics and breach morphology 
adopted in different mathematical models on dam breach are presented. 
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TABLE 2.1b 

TREATMENT OF COMPONENTS OF VARIOUS MATHEMATICAL 
MODELS FOR DAM BREACH EROSION 



Model 

Flow over 
the Crest 

Flow through 
the Breach 

Flow through the 
Breach on the 
Downstream face 

Sediment 

Transport 

Breach Shape 

Side Slopes 
Collapsing 

Cristofano Model 

N.A 

Broad Crested 
weir formula 

N.A 

Empirical 

Trapezoidal 

N A 

HW Model 

N.A 

Broad Crested 
weir formula 

N.A 

Modified 

Schoklitsch 

Parabolic shape 
with 45*’ 
side slopes 

N.A 

BRDAM Model 

N A 

Broad Crested . 
weir formula 

N.A 

Schoklitsch 
bed load formula 

Parabolic shape 
with 45*’ 
side slopes 

N.A 

PT Model 

N.A 

Relation between 
top width and 
flow rate 

Manning’s 

equation 

Mayei -Peter and 
Mull Cl equation + 

Ex net equation 

Constant Width 

N A 

Lou Model 

N A 

Dynamic 
wave modelling 

Dynamic 
wave modelling 

Du Boys’ bed load 
+ Einstein’s 
suspended 
load formula 

Cosine Shape 

N A 

Nogucira Mcxlel 

N.A 

Dynamic 
wave modelling 

Dynamic 
wave modelling 

May Cl -Peter and 
Mullei equation + 
Exnci continuity 
equation 

Cosine Shape 

Morphological 

mechanics 

DAMBRK Model 

Broad Crested 
weir formula 

Broad Crested 
weir formula 

N.A 

N A 

Rectangular 
trapezoidal or 
triangular 

N.A 

SMPDBK Model 

Broad Crested 
weir formula 

Broad Crested 
weir formula 

N A 

N A 

Rectangular 

shape 

N.A 

BREACH Model 

Broad Crested 
weir formula 

Manning’s 

equation 

Manning’s 

equation 

Mayei-Peter and 
Muller equation 

Rectangular 
Trapezoidal shape 

Balance of forces 

BEED Model 

Broad Crested 
weir formula 

Broad Crested 
w'cii formula 

Quasi -steady state 
approximation 

Einstein Brown 
formula 

Rectangular, 
triangular or 
trapezoidal shape 

Contour method 


2.6 Ancillary Literature Relating to Dam Breach 

The following literature are not having direct bearing on Dam breaching as such. But the 
information helps to understand some of the aspects of dam breaching. 

Annandale (1995) has enlightened the idea of erodibility with a rational approach for 
earth materials, ranging on a continuum from cohesionless granular soil through massive hard 
rock. His graphic relationship between the erodibility index of earth material and the rate of 
energy dissipation for the complete data set of 150 field observations and other published data 
on initiation of sediment motion exhibits an erodibility threshold that applies continuously 
across a broad range of materials. The erodibilty index is used to characterize the ability of 
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earth material to resist erosion. The rate of energy dissipation is used to characterize the 
relative magnitude of the erosive power of water. The erodibility index is dimensionless.The 
erodibility index is used to characterize several geological parameters. The rate of energy 
dissipation is calculated for a variety of flow conditions that can lead to erosion. 

Graf and Song(1995) have studied the behavior of bed shear stress in non-uniform and 
unsteady open channel flows. They have concluded the friction velocity V*, increases in 
longitudinal direction of flow in accelerating, but decreases in decelerating flow. In unsteady 
flow it reaches its maximum in the rising branch of the hydrograph. They also mention 
friction factor is essentially constant in a non-unifonn flow, but needs further investigation. 
In an unsteady flow the friction factor reaches its maximum in the rising branch, even before 
V», of the hydrograph and subsequently decreases. 

Aragon(1996), have developed hydraulic shear stress model for a uniform, steady flow on 
a slope of a maximum of water and granular material at high concentration draws attention, as 
the concept is new and the approach is elegant. His model includes stresses generated by 
particle-particle interactions along with these generated by fluid-bed interaction. Bagnold’s 
model of dispersive stress is modified by introducing Kinetic theory for granular flow. 
According to him an interesting finding from his experiment is; there is not a continuous 
increase in velocity with depth of flow as stated by Manning’s formula for whole range of 
concentrations. His hydraulic shear-stress model gives an insight about the mechanical 
behaviour of the flow along with velocity-flow depth relation. 

Lai and Shen(1996), in their research paper on ‘Flushing Sediment Through Reservoir’ 
have come out with some findings which it is felt is useful for the present study of dam 
breaching in particular to the morphological aspect of breach. From both laboratory and field 
data, they have found that the outflow sediment discharge in the flushing operation of 
sediment in the reservoir, can be well related to a hydraulic parameter which is a function of 
outlet discharge, water-surface gradient and flushing channel width. The equilibrium width of 
flushing channel is found to be proportional to Q °/ , where is the dominant outflow 
discharge. 
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2.7 Fuse Plug 

Tiney and Hsu(1962) conducted model study on fuse-plug in the laboratory adopting 1:20 and 
1 :40 scale. From their experimental analysis of the lab study data on fuse-plug washout, they 
found: 

a. The washout is a function of grain size of the materials the larger the size, the slower the 
rate of washout 

b. The washout rate for a fuse plug is between 1.05 and 1.2 times faster than the 
corresponding half-scale model. 

c. Increasing the volume of rockfill decreases the washout rate slightly. From this 
experimental observation and theoretical analysis, they have brought out these conclusions: 

• The mechanics of washout can be explained on the basis of the laws of sediment 
transport. 

• The rate of recession of the eroding face varies approximately with the 1/3 power of the 
scale ratio of fuse plug constructed with identical materials (non-cohesive). 

In addition to the above laboratory model test, they carried out large scale field model tests to 
verify the design assumptions like depth of water required to effect a successful breach, 
mechanism of the breaching action at the pilot channel, completion of breach before lateral 
erosion and washout process for lateral erosion. They also observed the effect of core on rate 
of washout, completeness of washout. Based on direct observations and the motion picture 
record, the rate of washout results are summarized as follows: 

i. Completion of breach: 1 min and 27sec. 

ii. Rate of washout (uniform material): S.bft/min or 0.0284m/sec. 

iii. Rate of washout (well graded material): 1.4 ft/min or 0.0071m/sec. 

iv. Hold-up time of the core was insignificant. 

V. Fuse-plug was completely washed out to the concrete foundation. 

Chee (1988) conducted a series of experiments in the laboratory to study the scale effects 
associated with modelling the erosion of granules embankments. In his investigation he 
carried out experiments with the smallest one measuring 150mm wide x 120 mm high to give 
a range with upstream slope 2:18 to 6:1 downstream slope 2:1 to 3:1 of model scales of up to 
1:12. A wide range of material sizes and embankment shapes were utilized in the test. The 
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model dams were breached by overtopping and washed away at constant discharges. In the 
washout of the dams, special care was taken to ensure that the bed material was not subjected 
to backwater effect. Photographs were continuously taken during the scour process to capture 
all the significant phases of erosion phase through the Plexiglas walls of the flume. The 
washout rate was computed from this pictorial time record by enlarging the photographs. In 
studying the washout time of the dams, he has used two separate methods namely, one from 
experimental dimensionless equation with numerical values of the coefficient obtained from 
the laboratory data by statistical analysis and secondly from well known sediment bed load 
equations. The time averaged values of the measured variables are used to formulate non- 
dimensional equation of the wash out rate. He has related C = sediment concentration 
with height of dam {Hq) critical flow depth for q,^, mean grain size (d 5 o), specific gravity of 
grain ( 5 ^), angle of repose of material in radians {0) and acceleration due to gravity (g), and 
mean width of dam has come out with a relation 



where E = initial distance between flume bed and dam crest, different to as dam is 
placed on a deck 

It has been shown, the experimental time scale for washout using Equation (2.10) depends on 


Op/6^ = 65 ., model scale {Lp/L„), dr= dp/d„ and S = (iS^-l) as: 



Conclusions from his analysis of Experimental data are: 

Washout time scale using Einstein bed load transport formula is independent of the model 
scale and depends on the properties of the granular material. Regarding the study of scale 
effects, he has also drawn the conclusion: for grain diameter scale, the model sediment grain, 
specific gravity decreases as the model material size increases and as the model sediment size 
decreases, the scale effect is reduced. 
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2.8 Scope for the Present Study 

From the present review of the literature it is found that the following aspects of dam breach 
problem should be addressed for the better understanding of the mechanism associated with 
its process. 

Dam breach problem due to overtopping flood lacks necessary experimental data support, 
particularly in the process behaviour. In the earlier experimental observations both lateral and 
vertical erosion process with time and space are not considered in the analysis. Reservoir 
capacity, inflow discharge to the reservoir and dam geometry and dam fill properties are not 
taken together for analysis in experimental data. 

Regarding numerical modelling on dam breach, quite a good number of models are 
available. Some of the models have predicted well the historical dam failure cases pertaining 
to breach outflow hydrograph. But the dam breach problem essentially associated with the 
erosion of sediment and its movement under highly accelerated flow condition. So the bed 
shear stress and sediment discharge relation needs a special attention to have a functional 
relation among them along with other breach parameters. In the literature the investigators 
have mentioned their reservations about the use of sediment transport formulae in the 
dynamic situations like dam breach which is meant to be used for alluvial stream for uniform 
flow situation. 

With these shortcomings in the present state of knowledge on dam breach due to 
overtopping of flood, the basic model studies are planned in the laboratory environment to 
study the dam breach process and to verify the obtained result with some field data on 
historical dam breach. 

In fuse plug washout study there is not much reported literature either on experimental 
study or numerical modelling except a very few to name. Although both the dam breach and 
fuse-plug washout problems more or less come under the same category of study, the former 
is three dimensional and latter is two dimensional. In case of fuse plug there is restriction of 
lateral erosion. Keeping this state of knowledge on fuse plug in view, basic model study (for 
the better understanding of the washout process in the fuse plug) is planned in the laboratory . 
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BREACH EROSION OF HOMOGENEOUS 
EARTHFILL DAMS DUE TO OVERTOPPING 
FLOW 

Part I deals with the erosion aspects of breach caused by overtopping flow on 
homogeneous earthfill dams. The details of experimentation carried out on the breach 
formation in homogeneous cohesionless earthen embankment and cohesive earthen 
embankment models are presented in Chapter III. 

The mechanism of dam breach process includes the morphological aspects, flow 
hydraulics and erosion characteristics. These are studied, first for cohesionless soil 
embankment and for cohesive soil embankment separately and later on the combined 
study has been made. Analysis of the laboratory data combined with field data is carried 
out for final breach geometry and flow parameters. These aspects are presented in 
Chapter fV. 

Chapter V includes further analysis of the lab data including two field data namely 
Teton dam breach and Huaccoto dam breach to predict the breach hydrogaph, sediment 
discharge graph along with final breach geometry and flow parameters. A methodology is 
suggested for the computation of breach process based on the analysis presented in 
Chapters IV and V. 



CHAPTER III 


SCHEME OF EXPERIMENTATION ON DAM 
BREACH 


3.1 Scheme of Experiments 

The breaching of an earthen embankment dam due to overtopping of flow is simulated by 
modelling earthen embankments with two different soils and varying embankment geometry 
in the Hydraulics laboratory of Indian Institute of Technology, Kanpur. The major objective of 
this experimental study was to understand the dam breach process due to overtopping flow. 

In the present chapter, details of flume, dam model, soil used for model preparation, 
measurement of criteria variables, the procedure to compute other criteria variables involved 
in the problem have been presented. Besides these the photographs taken during the breach for 
a dam model with cohesionless soil are presented to familiarise oneself with the different 
phases of breach process. Statistical analysis of the measured data is carried out to assess the 
magnitude of the deviations from the mean values. 

a. Details of Flume 

The flume used for these experiments is a tilting flume 20m long, 0.60m wide and 0.50m in 
height. The details of flume dimensions, position of embankment model in the flume and 
dimensional details of embankment are shown in Figures 3.1(a) and (b). The water was 
supplied to the flume from a constant overhead tank which got its supply from an 
underground sump, through a pump. Flow straightener was kept at the entry point of reservoir 
(storage area between the dam model and the water inlet section )to check the unsteady wave 
formation in the reservoir during the time of experiments . 
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Fig 3 . 1 bi SECTIONAL ELEVATION OF TILTING FLUME WITH DAM MODEL 
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Water supply into the reservoir was regulated by operating a valve in the inlet pipe. The test 
section for all the experiments was maintained at 8m from the inlet. A calibrated V-notch was 
used for measuring discharge passing through the flume after the breach flow comes to steady 
state. 

b. Model Details 

Two types of soils are used in the preparation of the embankment dam models, namely 
cohesionless (type-A) soil and cohesive soil (type-B). Henceforth the cohesionless soil and 
cohesive soil used for the model preparation will be referred as type-A and type-B soil 
respectively. 

i. Model with Cohesionless Soil 

The type-A soil, cohesionless soil was collected from the river bed of Ganga near Kanpur for 
model preparation. With this cohesionless sand, twelve experiments were carried out with the 
variation of model geometry and incoming discharge to the reservoir. The volume enclosed by 
the flume’s two side walls and the upstream of the dam model is called here as the reservoir. 
In all these experiments the inflow is kept constant throughout the experiment. The details of 
the model geometry are given in Table 3.1(a). The height of the dam (Hq) was varied from 
14.5 cm to 35 cm as shown in Table 3.1(a). Similarly the length of the dam crest is also varied 
from 15 cm to 156 cm, keeping the upstream and downstream slope of the dam constant to 1 
vertical (V) to 2.5 horizontal (H). 


TABLE 3.1a 

DIMENSIONS OF DAM MODEL WITH COHESIONLESS SOIL 


S. No 

Run Nos. 

Hd 

(cm) 

Bl 

(cm) 

L-TD 

(cm) 

V:H 

** 

Ldb 

(cm) 

Aq 

(era") 

Reservoir 
Capacity (cm^) 

(cm .'see) 

1 

1005 

35.0 

60 

150 

1:2.5 

325 

8312.50 

2363025.0 

420.00 

2 

wmm 

30.0 

60 

150 

MKBI 

300 

6750.00 

2025450.0 

1349.60 

3 

mimikJ 

25.0 

60 

mem 

MWHI 

279.0 

5412.50 

1687875.0 

1539.70 

4 

0505 

19.5 

60 

msm 

1:2.5 

247.5 

3875.63 

1316542.0 

569.70 

5 

hh 

14.5 

60 

145 

1:2.5 



978967.5 

291.00 

6 


15.0 

60 

156 

1:2.5 



1012725.0 

1586.60 

7 

MiHiM 

■09 

HI 

156 

1:2.5 


2787.63 

978967.5 

1616.40 

8 


■09 

HI 

154 

1:2,5 

MHll 



169.70 

9 


14.5 

HI 

142 

1:2.5 

214.5 

2584.63 


1586.60 

10 

MEW 

14.5 

HI 

60 




978967.5 

1036.13 

11 

2604* 

14.5 

HI 

15 

■09 

87.5 

743.13 

978967.5 

336.4 


2504* 

14.5 

60 

45 


117.5 

1178.13 

978967.5 

1866.60 


* These data are not used in the analysis due to insufficient informations.' 
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TABLE 3.1b 

DIMENSIONS OF DAM MODEL WITH COHESIVE SOIL 


S.No 

Run Nos 

Hd 

(cm) 

SCL 

(cm) 

Ltd 

(cm) 

V:H 

Ldb 

(cm) 

Aq 

(cn?) 

Reservoir 

Capacity(cm^) 

9'" 

(cm'Vsec) 

1 

1203 

12.0 

60 

15 

1:2.5 

75.0 

540.00 

7.97x10^ 

940.00 

2 

2703 

16.9 

60 

15 

1:2.5 

99.6 

968.50 

1.12x10* 

1129.85 

3 

0404 

20.0 

60 

15 

1:2.5 

115.2 

1305.80 

1.32x10* 

1082.20 

4 

1104 

23.6 

60 

15 

1:2.5 

133.1 

1750.39 

1.55x10* 

1330.00 

5 

2204 

26.4 

60 

15 

1:2.5 

137.2 

1859.20 

1.60x10* 

6088.60 


Upstream & downstream slope of dam cross-section V;H 


Notations used in Table 3.1(a) and (b) are explained as follows; 

//£, = height of the dam, 

Bci = crest length of dam model, 

= bottom width of the embankment along the centre line of the flume. 

Ltq = top width of embankment, 

V:H = bank slope of vertical : Horizontal 

y4£, = cross - sectional area of the embankment along the centre line of the flume 
Reservoir capacity = volume of water stored in the flume upstream of embankment 
= Inflow discharge into the reservoir. 

This slope of the dam was chosen keeping in view, that the ratio of Vertical to Horizantal 
(V;H) varies from 1:2 to 1:3 in the field. In some cases the top width of the embankment dam 
model is not proportionate to the earthen dam geometry. Larger top width is provided in some 
cases, because this gives enough time for the measurement of breach process. 

ii. Embankment Dam Model with Cohesive Soil 

The type-B soil, cohesive soil used for the preparation of embankment dam model was 
collected from the campus of Indian Institute of Technology, Kanpur. In total five experiments 
were carried out in this category. In case of cohesive soil embankment model the height of the 
dam varies from 12.0 cm to 26.4 cm, keeping the top width constant as 15 cm. The upstream 
and downstream slopes were maintained as V:H = 1:2.5. In this case top width was 
maintained proportionate to dam geometry. This is because the breaching process in cohesive 
soil model takes more time compared to cohesionless soil, giving sufficient time for 
measurement of breaching process. The details of dam models with cohesive soil is given in 
Table 3.1(b). The inflow was kept constant in these experiments also. 
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c. Properties of Soil used in Embankment Model 

Soil samples collected from each embankment of the dam were used for measurement of the 
angle of friction ( (j) ) and cohesion ( c ) by conducting undrained triaxial test in the Geotech 
laboratory of the Institute. These samples used for triaxial test are used for measurement of the 
dry density and water content of the soil after putting it inside the oven for 24 hours. The 
particle size distribution analysis was carried out for the soil used in the embankment fill. The 
type-A soil and type-B soil were subjected to sieve analysis for the size distribution. The 
details of soil properties are shown in Table 3.2. Moisture content was kept at 14 percentage 
during the construction of embankment dam model with type-A and type-B soils. 


TABLE 3.2 

SOIL PROPERTIES OF COHESIONLESS 
AND COHESIVE SOIL USED IN MODEL PREPARATION 


Type of soil 

<^50 

(mm) 

<1) 

(degree 

c 

KPA 

(Yb)av 

9810-^ 

N/m-' 

Cohesionless Soil 
(Type-A soil) 

0.16 

32.1 

0 

1.58 

Cohesive soil 
(Type-B soil) 

0.28 

33.1 

1.6 

1.95 


The particle size distribution curv'es for type-A and type-B soils are shown in the Figure3.1(c). 

d. Preparation of Embankment Dam Model 

To prepare the embankment dam model the following procedures were adopted for both 
cohesive and cohesionless soil. The soil was mixed properly after adding water to it. The 
amount of water added to the soil was judged purely based on the experience for better 
compaction of the soil. There after this soil was laid in the layers of 5 cm thick on the dry bed 
of the flume. Soil was compacted well with a hand hammer uniformly over the top surface. 
Scratches were made on the top surface before laying another layer. The next layer of the 
thickness 5 cm was spread over the surface and was again compacted. This process was 
continued till the model dimensions were achieved. Then the top surface, upstream and 
downstream face of the dam were compacted uniformly. Slight depression in the centre of the 
top of the model along the centre line of the flume was made in order to control the position 
of the breach occurrence. Then the soil samples were collected from the model body for the 
measurements of soil properties. The positions where soil samples were removed were again 
filled with soil with proper compaction.. The water was allowed upstream of the reservoir 
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keeping constant flow. As soon as the overtopping of flow occurs, the breach process sets in. 
Measurements of water surface levels, sediment bed levels and width of breach were carried 
out at prefixed positions along the breach longitudinal section. The methods adopted for the 
preparation of embankment dam, both for cohesive and cohesionless soil were same for all the 
runs conducted. 

3.2 Some Observations of the Breach Process 

Some critical observations made during the breaching process of cohesionless and cohesive 
soil embankment dam models are; 

a. Observations of Erosion Process due to Overtopping of Flow in Cohesionless Soil 

Embankment Dam 

i. As soon as the incoming flow or reservoir water level exceeds the embankment level, 
water starts over flowing first on the dam crest and then on the downstream slope as 
indicated in the photographs, shown in Figures 3.2(a), and (b) respectively. These 
photographs are taken from downstream end sequentially, away from the 
embankment, looking in the upstream direction. The height of this model was 0.25m 
and top width of the crest of the embankment was 0.15m with both upstream and 
downstream slope equal to V:H=1:2.5. 

ii. As the flow traverses the downstream slope completely and touches the bed, the 
erosion near the toe of the embankment is initiated. Lateral widening of channel at 
the toe of the embankment and deepening of channel on the downstream face are 
visible as in Figures 3.2(c) and (d). Erosion on the flow path particularly in the 
downstream is initiated because of steep gradient associated with high shear stress. 

iii. Toe erosion starts receding upstream and by widening its path as well as having 
steeper vertical fall. These can be seen in Figures 3.2 (e) and (f). During this process, 
the flow channel formed at the dam crest and downstream face of the embankment 
deepens having almost v - shaped valley or narrow deep valley. 
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Fig. 3.2(d) : Lateral enlargement of erosion at toe of the dam and formation of 
gully on downstream slope of the embankment 
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Fig. 3.2(e) : Receding of toe erosion in upstream direction 



Fig. 3.2(f) : Deep vertical cut in the downstream face of the dam along with 
receding of toe erosion in upstream direction 
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iv. Formation of bank collapse and washout of collapsed bank soil, formation 
contracting section upstream and diverging section downstream with contracted 
throat are seen in Figures 3.2(g) and (h). 

V. During the receding crest of the flow, the flow channel gets widened in the 
downstream face with possible bank collapse. This phenomenon can be seen from 
Figures 3.2(i), (j), (k), and (1). 

vi. During this receding process of crest, one can observe that the flow over the crest 
behaves like flow over a erodible weir with gradual change in water surface in the 
upstream and steep drop in water surface level over the crest. On the downstream end 
one can notice pool formation and sometimes flow associated hydraulic jump. This is 
shown in Figures 3.2(m), (n), and (o). 

vii. Once the receding of the crest moves in the top section of the embankment dam, the 
erosion process is fast, causing bank collapse and bed erosion. This can be seen from 
Figure 3.2(p). 

viii. As soon as the flow crest moves upstream, beyond the upstream edge of the 
embankment ,the crest width enlarges. The alignment of the crest is curved. This 
causes the breach section having a large width at upstream edge, contraction at the 
centre and enlargement of width at down stream end. This can be seen in Figures 
3.2G)to3.2(p). 

ix. By this time, water level in the reservoir has lowered considerably. The rate of 
erosion process slightly reduces. The flow at upstream face will be similar to flow 
over spillway with erodible crest. During this process, the crest level reduces and 
bank collapses causing widening of the breach, as shown in Figures 3.2 (q) to 3.2 (r). 

X. The above mentioned process continues till the crest is almost washed out and 
channelisation of flow through breach section forms. This is shown in Figures 3.2(q) 
and 3.2(r).This process continues till emptying of reservoir or occurrence of gradual 
varied flow having a very minimal erosion. This erosion is similar to the erosion in 
channel with gradually varied flow. 
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Fig. 3.2(g) : Widening of the path of erosion on downstream slope of the 
dam with steeper vertical fall and bank collapse 



Fig. 3.2(h) : Under scouring on the downstream slope of the dam with breach channel 

having contraction section upstream and diverging section downstream with 
a throat are visible 
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Fig. 3.2(i) : Widening of the breach channel by bank collapse and spillway 
type crest formation are visible 



Fig. 3.2(j) : Washout of soil from bank collapse with contracted throat 
and spillway type of crest flow at upstream are visible 
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Fig. 3.2(k) : Widening of flow path of water, contracted throat and 
spillway type of crest flow are visible 



Fig. 3.2(1) : Spillway type of flow over the erodible crest with large width at upstream 
edge, contraction at the centre and enlargement of breach width at 
downstream end are visible 
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Fig. 3.2(m): Flow over the crest behaves like flow over erodible spillway 



Fig. 3.2(n): Pool formation at downstream of erodible crest and flow 
associated hydraulic jump 
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Fig. 3.2(o): Widening of breach with pool formation at downstream of the crest, and 
hydraulic jump formation at the erodible spillway at the upstream 



Fig. 3.2(p): Curved crest spillway type flow, hydraulic jump form at the toe 
of the spillway flow, more lateral erosion on side bank are visible 
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Fig. 3.2(q): Reduction in reservoir level, lateral erosion touching one side of 
the flume, flow concentration on one side are visible 



Fig. 3.2(r); State of breach process nearing to the final state of erosion is observed 
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The whole process of erosion of breach occurs within around 5 to 10 minutes in laboratory 
for dam model with cohesionless soil. 

b. Erosion Process due to Overtopping of Flow on the Cohesive Soil Embankment 
All the steps of erosion process indicated in cohesionless soil embankment are occurring in 
the erosion process of cohesive soil embankment also. The main difference in the erosion 
process of the embankment with cohesive soil is slower to cohesionless soil embankment 
erosion. Duration of breach observed for the cohesive soil embankment in the laboratory is 
about 20-45 minutes depending upon the intensity of inflow and channel geometry of the 
embankment model. The other main observation is the occurrence of the cascading of the flow 
near the downstream edge of the embankment during part of erosion process. 

3.3 Data Collection 

a. Measurement of Basic Data During Breach Erosion 

To study the breaching process of the dam due to overtopping of the flow, it is necessary to 
measure the water level(hvi,), bed level(h(;) and width of the breach at different longitudinal 
positions with respect to time. Vernier point gauge with an accuracy of 0.10mm has been used 
for measurements of h^^,, and width of the top surface of water has been measured by an 
ordinary scale with an accuracy of 1.0mm. The vernier point gauge was modified for the 
measurement of h,. by providing a flat bottom at tip of the pointer by fixing a thin cork. From 
this basic data one can compute water and sediment discharge and breach geometry' for given 
flow condition. During the measurement of these breach variables for cohesionless 
embankment it was observed that the erosion process w'as very fast. In order to increase the 
erosion time, the top width of the model was increased from 15 cm to 156 cm. 

The position of the measurement of water surface, bed level, and breach width are 
indicated in Figure 3.3(a). Table 3.3(a) contains details of measured quantities at each section. 
Analysis of erosion process was carried out after completing the water surface, bed surface 
and breach width variations with time. A few shortcomings were observed during the analysis 
of dam breach model data with cohesionless soil fill, particularly the absence of bed level 
measurements at dowmstream end, and water surface level variations for upstream of the 
embankment. These shortcomings were taken care of by installing more measuring stations 
for the embankment models prepared out of cohesive soil. 
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Fig 3,3a i LONGITUDINAL CROSS-SECTION OF EMBANKMENT DAM MODEL WITH CDHESIONLESS 
SOIL INDICATING LOCATION OF SELECTED MEASURING SECTIONS 
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Figure 3.3(b) indicates the positions at which measurements were carried out for cohesive soil 
embankment models. Table 3.3(b) gives the details of measured quantities at different sections 
for cohesive soil embankment model. 

b. Final Measurements after the Attainment of Complete Breach Process 

The bed level, water surface levels, breach width and cross-sections of the breach at different 
locations in addition to the specified measuring sections were carried out. By end of the 
breach, flow in the breach section has reached steady state. The measurement of discharge 
through breach was carried out by measuring the discharge over the v-notch attached to the 
tank located at downstream end of the outlet flume. 

TABLE 3.3a 


VARIABLES MEASURED DURING DAM BREACH AT DIFFERENT 
SECTIONS OF DAM MODELS WITH COHESIONLESS SOIL 


Section 

Hi 

2 

(B, 10 cm 

3 

® 60 


mgm 

6 

Downstream 

Remarks 

1005 

WEM 

s w 

SB 

SB 

SB 

w 


0805 

HSfl 

s WB 

SB 

SB 

S^ ' 

w 


0705 



SB 

SB 

SB 

w 



w 

S WB 

SB 

SB 

SB 

w 


■lEbfcM 

w 

S WB 

SB 

SB 

SB 

w 


■tltliM 

w 

WB ' 

- 

■SOI 

SB 

w 


WEBm 

w 



"SB 

- 

w 




sw 

w 

S 

- 

w 



— 

■SBi 

SB 

SB 

- 

w 



w 

WSB 


- 

- 

WB 



TABLE 3.3b 

VARIABLES MEASURED AT DIFFERENT 
SECTIONS OF DAM MODELS WITH COHESIVE SOIL 



1 

2 

(aa) 


4 

5 

6 

7 

8 

HDSBi 

w 

w 

w 

WSB 



WSB 

WSB 

HBHHi 

w 

w 

w 

Kaasi 



WSB 

WSB 

jH^HH 

w 

w 

w 

B15i51 



WSB 

WSB 

HDEEHI 

w 

w i 

w 






2204 

w 

w 

w 




WSB 

WSB 


Where W Water level 

S ==> Sediment Bed level 

B => Breach width at corresponding water surface. 


Hence forth the section of measurement of water level in the reservoir, section- 1 and section- 
2 for cohesipnless and cohesive soil embankment rspectively will be called as section-(aa) in 
the course of discussion. 
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Fig 3.3bi LONGITUDINAL CROSS-SECTION SHOVING LOCATION OF SELECTED 
MEASURING SECTIONS FOR COHESIVE SOIL EMBANKMENT MODEL 
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c, Basic Data Presentation for the Embankment Dam Model of Cohesionless Soil 
Basic data collected for the run no.-1005 is shown as a typical set. The run no.- 1005 means, 
the experiments were carried out on the 10th day of the 5th month of the year 1995. The 
variation of water level in the reservoir for section- 1 and sediment bed level for the 
section-2 is shown in Figure 3.4(a). Variations of water level and sediment bed level 

) at section - 2 are shown in Figure 3.4(b). The variation of water and sediment bed level 
are shown in Figure 3.4(c) for section-5. Figure 3.4 (d) shows the water level variation for 
section-6. The sediment bed level variations are shown in Figures 3.4(e) and (f) for sections 3 
and section 4 respectively. The variation of breach width are shown in Figures 3.4(g), (h), (i) 
and (j) for sections-2, 3, 4 and 5 respectively. 

d. Basic Data Presentation for the Embankment Dam Model of Cohesive Soil 

The basic data collected for the run no. -0404 is presented as a typical set. Here the run number 
0404 means the experiment was conducted on the 4th day of the 4th month of the year 1996. 
Figure 3.5(a), contains the variation of reservoir levels with time for sections 1 and 2. Figure 
3.5(b) contains the variation of the water surface level for section-3. The variations of water 
surface level and bed level are shown in Figure 3.5(c) for section-4. The Figures 3.5(d), (e), 
(f), and (g), contain variations of water surface level with time for sections 5 , 6, 7 and 8 
respectively. The breach wddths at water surface level for sections 4, 5, 6, 7 and 8 are 
presented in the Figures 3.5(h), (i), (j), (k), and (1) respectively. 

3.4 Computation of Parameters Controlling Breach Process 

The parameters controlling the breach process can be broadly classified into two broad 
categories, namely flow parameters and sediment parameters. Definitions of some of these 
parameters are shown in Figure 3.6. 
a. Flow Parameters 

The following are the flow parameters viz., , Q,,., Sj-, F, ^ AE. 

: the difference in the water level at section-(aa) to the sediment crest at breach. The 
reservoir level at section-a is chosen based on the water surface variation from the inlet to 
the breach section. It was observed in most of the cases that the maximum water level was 
occurring at section 1 in the case of cohesionless embankment and section 2 in the case of 



Scheme of Experimentation on Dam Breach 51 


cohesive embankment. This level was considered for discharge computation. 

Q^.: the water discharge flowed through the breach at any instant of breach process. This is 
computed using flow continuity equation as shown below. 

[(& A0]= Change in volume of water in the reservoir during time At 

where , Q^. = Water discharge through breach section 

Qi„ = Inflow discharge in to the reservoir, which is kept constant 
throughout the experiment and is measured after complete 
breaching of the embankment. 

The change in the reservoir volume is computed by knowing change in water surface level 
from the variation of the reserv'oir level and multiplying it with water surface plan area. 


AE : total energy difference between energy level at section-(aa) and downstream toe section 
of breach. Velocity at toe section is computed by knowing the cross-sectional flow area 
and its corresponding breach discharge. 


: Energy slope. This is the difference in total energy between the upstream section-aa to the 

down stream (toe) of the dam divided by the inclined sediment bed length between 
upstream edge of the dam to the toe of the dam. where , 

I,, = (Ah,^ + 


Ahp. Difference in sediment bed level between the upstream edge to the toe of the 
embankment. 

L : Horizontal distance between the above mentioned sections. 

In this analysis it is assumed that the energy at upstream edge of the embankment is equal 
to the energy at the section-(aa). 

V, : Average bed shear velocity and is computed as , 


_ I- _ _ -. 1/2 

Where, g = acceleration due to gravity 




L' 

i' I. 

V'- 




121596 
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20 


RUN NO. - 0404 

Section — 3 
Water Surface Level (WL) 


0 1000 

Time (sec) 

Fig. 3. 5a: Variation of Reservoir 

Water Levels with Time at 
Sections-I and 2 



0 1000 ; 

Time (sec) 

Fig. 3. 5c: Variation of Water Surface 
and Sediment Bed Levels 
at Section— 4 with Time 


000 2000 3000 

Time (sec) 

Fig. 3. 5b: Variation of Water Surface 

Levels with Time at Section~3 


Time (sec) 

Fig-3. 4e: Variation of Water Surface 
.and Sediment Bed Levels 
at Section— 5 with time 
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Fig.3.5e: Variation of Water Surface 
and Sediment Bed Levels 
with Time at Section-6 
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Fig.3.5f: Variation of Water Surface 
and Sediment Bed Levels 
with Time ot Section— 7 
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Fig.3-5g: Variation of water surface Flg.3.5h: Variation of breach 
and sediment bed levels width at section— 4 

with time at section-8 with time 




Fig.3.5i: Variation of breach 
width at section— 5 
with time 


Fig.3.5j; Variation of breach 
width at section— 6 
with time 




Fig. 3, 5k: Variation of breach 
width at section— 7 
with time 


Fig. 3. 51: Variation of breach 
width at section— 8 
with time 
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Fig. 3.6: Definition sketch for the dam breach analysis 
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: Average value of the hydraulic radius at three sections, namely at the upstream 
top edge, downstream edge of the embankment and at the toe of the breach. 

AH; Difference of water surface level between upstream of embankment (section 1, for 
cohesionless embankment section 2 for cohesive embankment) to the water level at 
downstream toe of the embankment. 
hc.„: Height of crest level measured from the bed of embankment. 

: Height of the sediment profile at a section x measured from upstream edge of the dam. 
h[ : Height of sediment profile at the downstream top edge of the embankment. 

: Reservoir level measured at section-(aa) 

; Height of embankment dam. 
b. Sediment Parameters 
The sediment parameters are A^, and C. 

A^'. This is defined as the eroded breach cross sectional area. This quantity is computed as the 
product of depth of erosion and the average width of erosion at a section considered. This 
quantity is computed at three sections namely upstream edge, downstream edge of the 
embankment and toe of the dam. 

\//. Volume of sediment eroded. This quantity is computed as the product of average eroded 
area between the two sections and the length between them. In the case of cohesive soil 
embankment models, more sections are used for the computation of . 

: Sediment discharge at the end of breach section. The cumulative volume of sediment 
eroded is plotted with time and an average curve is fitted to this data. The ratio of 
difference in cumulative eroded sediment volume divided by the corresponding time 
interval is considered as the total sediment discharge out of the breach section. 

C : Sediment concentration. This quantity is computed by dividing sediment discharge by 
breach outflow water discharge at a given time. 
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3.5 Frequency Analysis of Deviation of Data from the Mean Value 

Measurement of water surface level, sediment bed level and water surface width for a section 
are plotted with time. Typical variations of these measurements may be seen in Figures 3.4 
and 3.5. Fluctuations in data measured are clearly visible. A mean curve was drawn through 
these measured data using regression analysis. The data measured from this mean curve are 
taken for further analysis, not the actual fluctuating values. In order to see how much error 
one makes by using mean value instead of fluctuating values, a frequency analysis for the 
deviations from actual value to mean value is carried out. Typical variations of these 
frequency plots are shown in Figures 3.7 (a), (b), 3.8(a), (b), and 3.9(a), (b) for measurement 
of sediment bed level, {h^, water surface levels {h^) and width of breach at water surface 
level. The standard deviation indicates the average order of fluctuations associated with the 
measured data. The left hand side and right hand side crossed points of the dotted line in 
Figures 3.7(b), 3.8(b), and 3.9(b) shows /i-cr and fi+cr respectively. In x-axis of these figures, 
range means difference between the observed and fitted value of the variable. The symetrical 
frequency plot with high peak and narrow base around mean indicates there is least difference 
between observed and fitted value of the variable. The unit of /u, and a in these figures are 
expressed in centimeters. The magnitude of standard deviation in the measurement of 
sediment bed is high in comparison to other measurements. This is because erosion and 
deposition process associated with bank collapse and also with some uncertainty in locating 
sediment bed level. Since the mean curve obtained through regression analysis is used for 
extracting data for further analysis, error involved in the data will be minimum. 

3.6 Concluding Remarks 

In the present chapter besides describing the experimental setup and procedure in details, a 
typical set of measured variables like water level, sediment bed level, breach width at selected 
sections are presented for dam model of both cohesive and cohesionless soil fill during breach 
process. The approach adopted for the computation of other breach criteria variables are 
demonstrated. A typical frequency analysis of the deviations of experimental data from the 
mean curve are presented to access the order of the fluctuations associated with data. 

Since the data extracted are from the mean curve, the error in the analysis of breach 
process is hoped to be minimum. 
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CHAPTER IV 


MECHANISM OF DAM BREACH PROCESS 


4.1 General 

In this chapter effort has been made to know the details about the mechanism involved in the 
dam breach process, with the help of experimental data obtained from the model studies. The 
data on water level, sediment bed levels, and width of breach are measured both in time and 
space during the dam breach process due to overtopping flood for the homogeneous dam 
models of cohesive and cohesionless soil with varying flow parameters and dam geometry. 
The effect of dam breach causes variation in the discharge of water and sediment at the 
downstream end of the dam. These two variables are very important so far their routing in the 
downstream of dam is concerned. To estimate these two quantities with time, it is essential to 
know the evolution of breach process, particularly the breach morphology and breach 
hydraulics with time. The study on dam breach process involves to know the dynamic 
variations of the following criterion variables, namely: 

a. Variation of crest and the sediment bed level in the longitudinal direction of breach. 

b. Variation of breach width at crest and at other sections along the breach. 

c. The rate of variation of crest height and its width. 

d. Water level variations at different sections of the breach. 

These are the basic criterion variables to which other related criterion variables like energy 
loss, sediment discharge, water discharge are related. 
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To know these stated criterion variables, the available predictor variables in dam breach 
process are; 

a. Inflow hydrograph to the reservoir 

b. Capacity of the reservoir (reservoir-water elevation relation) 

c. Dam geometry 

These two types of variables, in a more common way the input and output can be related 
through parametric study with the properties of the dam fill material like angle of friction, 
cohesion, median particle size, and bulk density of dam fill before overtopping. 

In the present simulation of physical scaled model study both input and output are 
measurable. To study the governing mechanism involved, the cause and effect can be 
functionally related through suitable controlling parameters. 

4.2 Parametric Formulation 

In an experimental study one has to resort to nondimensional parametric study to take care of 
the scale effects. The results so obtained from the analysis of the experimental data can be 
applied to any prototype situation. Considering this aspect in view and knowing the possible 
variables controlling the dam breach process, the following nondimensional parameters are 
used, which can be grouped under the following head. The S3mibols are defined in the list of 
symbols. 

a. Breach Hydraulics 
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c. Time Components of Breach 
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d. Energy Loss Parameters 
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e. Sediment Transport Parameters 
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Some specially coined nondimensional numbers are used frequently in the present study. 

These nondimensional quantities are presented here with an intention to acquaint them before 
their appearance in the analysis. The individual descriptions of these parameters are made 
wherever they are used. 
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4.3 Analysis of Breach Process in Cohesionless Soil Embankment 

Cohesionless embankment breach data are analysed first. From this analysis one gets 
information on the variation between the predicted and criterion variables. With this 
knowledge, identification of major controlling variables for the analysis of cohesive soil 
embankment breach data is carried out . This analysis reveals the importance of cohesion in 
the breach process. Finally an attempt has been made to analyse the breach process by 
combining both the model data on cohesionless and cohesive soil embankments. Field data 
available in literature are used along with laboratory data to develop the relation to predict the 
important variables of breach process. A method is developed to predict the breach processes 
for a given inflow, dam geometry reservoir capacity and embankment material properties. 

Analysis of breach processes involves study of breach morphology and breach hydraulics. 
Breach morphology involves variation of water surface, sediment bed surface and width of 
the channel with time during breach process. Water discharge relations, sediment discharge 
computation, energy loss, bed shear stress and their internal relationships are studied in the 
breach hydraulics. Breach erosion involves the variation of crest height, crest width and their 
rate of erosion, volume of sediment eroded and transported out of breach section. 

Cohesionless soil embankment model are built using river sand collected from river 
Ganga near Kanpur, India as detailed in Chapter III. As soon as the flow overtops the 
embankment, erosion of soil sets in. Erosion of the soil on the embankment body leads to 
development of channelisation. Erosion of the soil at the toe of the dam leads to the formation 
of the crest, where flow drops suddenly. With increase in the magnitude of overtopping flow, 
the crest starts receding upstream and finally reaching the upstream face of the embankment 
dam. During the receding of the crest towards upstream direction, the widening of the channel 
bed and bank occurs by erosion. This process of channel formation is defined here as the 
breach process. Development of breach is analysed by studying the morphology of the breach 
first and later the breach hydraulics. The formation of the eroding crest during breach 
development may be looked upon as a erosive crest similar to spillway or weir. Analysis of 
flow relation with crest geometry and flow depth are carried out. During the process of crest 
receding, widening of the crest and erosion of the crest occurs simultaneously. Hence, the rate 
at which crest erosion or widening takes place is studied further. Growdh of the crest width 
and reduction in crest level with variation of the upstream and downstream flow conditions 
are analysed. The volume of water flown out of breach and the amount of the sediment 
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eroded in the breach are analysed next. The erosion of sediment depends mainly on the 
excessive erosive force caused by the flow. This aspect is analysed by computing the average 
shear stress (r/) of flow and relating with sediment concentration {C=Qs/Q^ ) at a given 
instant of time. Work done by the flow by eroding the sediment causes loss in energy. This 
loss in energy is related to the flow condition upstream and downstream of the embankment. 
Finally, the time taken for breach process and final breach dimensions are analysed. 

a. Breach Morphology 

Morphology of breach contains variation of planform of the breach, longitudinal variation of 
reach bed level, w'ater level, and cross-sectional variation of breach sections. 

i. Morphological Variations of Breach 

Breach morphology is related to the variation of breach geometry with time. Variation of 
breach geometry depends on outflow hydrograph. For run 0405, morphology is presented for 
a typical outflow hydrograph as shown in Figure. 4.1 . 
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Morphological variables like the water surface, sediment bed surface and breach width before, 
at and after the occurrence of the breach peak discharge are presented in Figures 4.2 (a), (b), 
(c) and (d) respectively as a typical set. These Figures indicate a gradual variation of the water 
surface upstream of the crest and rapid variation after the crest. The sediment surface level 
remains almost parallel to the water surface level. The gradual widening of the breach width at 
the upstream section, followed by a narrow channel width where water surface drops rapidly 
is observed. Widening of the breach width in the downstream section 6, 7 where rapid 
variation in the water level occurs is observed in every run. After the occurrence of peak 
flow, as the breach process progresses the channel contraction occurs right at or near the 
upstream face of the embankment. The gradual widening of the breach occurs in the 
downstream direction from the upstream face of embankment upto the downstream face of 
the embankment. These features are observed in all the cohesionless embankment models. 
The bed level and water surface level in the breach section at the end of breach remains 
similar to the gradually varied flow in the channel. 

ii. Breach Cross Sectional Variation with Time 

The cross sectional area eroded due to breach flow at upstream edge and down stream edge of 
embankment dam for four different instants of time during the breach are plotted in Figure 
4.3. It may be observed that the variation of cross-sectional shape during most of the breach 
erosion process for cohesionless soil remains rectangular. Due to the instability of the slope, 
banks may collapse. The bank soil collected in the breach channel due to bank collapse will 
get eroded quickly, retaining back the unstable rectangular shape. This process continues till 
the final stage of breach occurs. At this stage, the breach section remains more of trapezoidal 
shape. The area of erosion at the top upstream edge of the embankment is presented as a 
function of time for run 0405 as a typical case as shown in Figure 4.4. The volume of 
sediment eroded is plotted against time of erosion in Figure 4.5 for the run no. 0405. It reveals 
that the slope of the volume of erosion curve is steeper in the beginning and becomes fairly 
constant thereafter. 
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Fig. 4.3: Variation of breach cross-sectional area before, at and after peak discharge 
at sections-2 and 5 for dam model with cohesionless soil, for run no. 0405 
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Fig. 4.4: Variation of eroded area at the upstream edge 
of the embankment for run no. 0405 with time. 



Fig 4 5: Cumulative volume of eroded sediment from breach 
section of the dam body for run no. 1005 
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b. Breach Hydraulics 

Breach hydraulics contain the discharge computation through the breach, analysis of 
coefficient of discharge {cf) of the assumed broad-crested weir type flow situation in dam 
breach, the analysis of head (h,) acting over the crest and maximum height of the erodible 
sediment bed. Variation of average velocity, breach hydrograph, energy loss and bed shear 
velocity are also analysed. 

i. Discharge Characteristics 

The photos on breach process presented in Chapter III reveals that the flow through the breach 
behaves like flow over a spillway. From this consideration, the flow through the breach is 
analysed as the flow over a spillway in which crest is getting eroded. In order to study what 
part of time during breach process the flow behaves like flow over spillway and what part of 
time flow behaves like a channel flow, it is essential to analyse the variations of reservoir 
water level, water level at downstream (near toe of the dam), and the crest variation during the 
breach process .These variations are shown in Figure 4.6. From this figure one can observe, 
when sediment crest level equals to the downstream water level by getting eroded, till then 
flow can be treated as firee overflow similar to spillway flow. Till the downstream water level 
fairly coincides with upstream water level, it can be treated as weir flow under submerged 
condition. After the downstream water level fairly coincides with upstream water level, the 
state of flow in the breach is similar to channel flow. The discharge characteristics of the 
overtopping flow is analysed firom these considerations. Discharge intensity , discharge 
per unit width of breach at the upstream top edge of the dam) is plotted against , the 
difference between upstream reservoir level and the sediment crest level as shown in Figure 
4.7(a). By using the regression analysis, the power law relation is developed for this data and 
it indicates, 

(Ala) 

However the scatter of the data is quite large. The coefficient of discharge {c^ ) is computed 
based on the weir type flow as: 




(4.1b) 
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is plotted against the ratio of ( where is the height of the crest above the bed of 
the flume. The Figure 4.7(b) indicate that q^/ flghl is inversely proportional to . 

However, its relation is given below: 



= 0.06x 



- 0.87 


(4.1c) 


The frequency analysis for the data about the mean line is analysed and plotted in the Figure 
4.8. The mean ( // ) and standard deviation (cr) for this frequency analysis are: 

// =0.0218 
(j =0.0165 

Based on the classification of the flow as carried out in Figure 4.6, and relation shown in 
Figure 4.7(b), it may be postulated that the free weir flow exists for < 0.4. From the 

Figure 4.7(b) the order of magnitude of variation is 0.7-0.22 for free overflow. 


ii. Computation of Head of Water Acting Over the Weir Crest (/i„) 

As a representative case, the variation of /i^ with time for run-1005 is plotted in the Figure 4.9. 
It may be observed that h^, starts increasing with time and attains a maximum value. With 
further increase in time starts decreasing asymptotically and reaches certain value 
depending on the breach geometry and amount of incoming flow to the reservoir. The 
variation of with time is mainly during free overflow. From this observation it was planned 
to analyse with time, as h^^ / and t / t},up, where h^p and are corresponding 
maximum value of and its time of occurrence. Using the data of all the runs for 
cohesionless soil embankment, /i„ / h^p is plotted against t in Figure 4.10(a). The relation 
between them is given as: 

(4.2a) 

K 

where = 2.5 (rising limb) and /7y= 1.0, 1.2, 1.25 for falling limb as shown in Figure 4.10(a). 
As already indicated the variation / h^p follows a steep rising limb and then gradual 
receding limb. The position at which the final value in the receding limb reaches depend on 
the magnitude of the incoming flow to the breach and breach geometry. The h^p is found to be 
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function of height of the dam inflow to the reservoir ( ) snd cross sectional area 

( ) of the embankment. Using these variables, a dimensional analysis is carried out and it is 

found that h,p /Hi^ is related to the Q.^ / ^2gH^ . This relation is plotted in the Figure 

4.10(b). It may be observed that the curve gradually increases with increase in X-axis value 
and reaches maximum value around two. It decreases very steeply with further increase in its 
X-axis value. Using regression analysis an equation relating these is given as: 

(K / )= 0-2068 + 288.7X - 323 1 03X^ (4.2b) 

where, X = 

The time scale ( ) is nondimensionalised with and is plotted against nondimensional 

as shown in Figure 4.10(c).The functional relation among them is written as: 

^ = X 0.1 85 (4.2c) 

^95 

where X = — ^ 

With computed value of h^p from Equation(4.2b), ti,„p can be found out using Equation(4.2c). 
With reference to Figure 4.5, t 95 can be determined. 

Hi. Breach Average Velocity 

The average breach velocity ( ) is defined as the breach outflow discharge {Qw) at any time 
divided by its corresponding average flow cross sectional area, {Aj-). The average cross 

sectional area is the area average of the cross sectional areas computed at upstream edge, 
downstream edge of the embankment and at the centre of these two. This average velocity is 
plotted against average hydraulic radius ( ) computed at three different sections is stated. It 
may be noted that the average velocity decreases very slowly with increasing as indicated 
in the Figure 4.1 1(a). Same average velocity is plotted against energy gradient, , computed 
between the upstream edge and the downstream edge of the embankment. It may be observed 
that average velocity increases with Sj in Figure 4.11(b). The relations of 5^ and with 

average velocity does not match either Chezy’s or Manning’s equation. 
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(Qw/(A,)ov)=((R.)J-''“““*37.56 


Fig. 4.11(a): Variation of overage breach flow velocity(y„) 
with average breach hydraulic radius (Rhjav 


(Qw/(A,)ov) = (Sf)a,"*99.1 


• • 7 • • 


4.11(b); Variation of average breach flow _velocity(V,,) 
with average breach energy gradient(Sf)av 
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The average velocity in the dimensionless form 3s[QJA^)/fgR^Sj. is plotted against 

non-dimensional water level variation at crest as shown in Figure 4.11(c). From this figure it 
may be observed that the dimensionless average velocity increases with increase in the 
magnitude of drop in the water level at the crest. The average velocity is further analysed in 

relation to Manning’s equation. The magnitude of Q,y / Aj. is plotted against ^ and 

also is plotted against x ” obtained from the regression analysis. Both the relations 


found to coincide each other, indicating mutual balance between Sj- and . as indicated in 

the Figure 4.11(d). The data as plotted in Figure 4.11(d), when reduced appropriately 
simplifies similar to Manning’s relation: 

S^ = ^;’“x5;'^^^x43.8 (4.3a) 

A^ 

From regression analysis. 


Q. 


= (5 


0.09 / 0^4^x45.35 


(4.3b) 


Overall this relation indicates nonconfirmity to Manning’s equation in breach flow. This may 

be due to flow accelaration associated during the breach. 

iv. Flow Velocity at the Downstream Top Edge of the Embankment 

The VI is the velocity of flow at the downstream top edge of the embankment. Analysis of 

F; is carried out by introducing critical velocity, (gy’^) for nondimensionalising FJ and is 

plotted gainst average water surface slope, with same Q^rp - number as shown in Figure 

4.12(a). Here is the critical depth of flow at the downstream top edge of the embankment. 

Its magnitude is computed assuming wide rectangular flow as: 


= 


g 


1 

2^3 


(4.4a) 


where ql = water discharge intensity at downstream top edge of the dam 




ig. 4.11(d): Comparison of Manning's Equation in breach flow 
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Relation between them is obtained from regression analysis as given below. This relation 
involves the knowledge of the water intensity (^J. 

f ^ 

c *2, 

S„x — ^ (4.4b) 

On similar lines the non-dimensional flow velocity at the downstream top edge (FJ ) of the 
embankment with corresponding shear velocity at that section V,' is plotted with as shown 
in Figure 4.12(b). The relation obtained from regression analysis enables to predict the shear 
velocity at the downstream end for a known flow condition. The functional relation among 
them is given as follows: 

( Q V'”r 

(4.4c) 

From Figure 4.12(b) it can be observed that for low values of the scattered data deviates 
more , whereas for higher the scatter in the data is less. This indicates that the relation is 
good during initial stages of breach process, 
v. Breach Hydrograph 

Breach hydrograph is the relation between outflow through the breach plotted against time, as 
shown in the Figure 4.13(a). This hydrograph has a steep rising limb with a peak followed by 
a gradual receding limb. The receding limb discharge after a long duration reaches equal to 
the magnitude of the inflow discharge into the reservoir. Considering the variation in inflow 
discharge it was thought to analyse rising limb and receding limb separately. The rising limb 
is non-dimensionalised by peak breach outflow discharge (Q^p), whereas receding limb is 

non-dimensionalised hy (j2w ~ Qin)^ (Q/tp ~ Qm)- magnitudes are plotted in Figure 

4.13(b) with nondimensional time where V' is the time under consideration and t^p is the 
time of occurrence of peak outflow discharge. 


V 

= 161.64 X 

K 




Mechanism of Dam Breach Process 83 



Fin 4 12fb'): Variation of the ratio of downstreann top edge breach velocity to 

^ downnstream shear velocity with average water surface gradient 
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Fig. 4.13(a): A typical breach outflow hydrograph for run 


no. 1005 
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Fig. 4, 13(b); Nondimensional breach outflow with nondimensional time 


Receding Limb 
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The relations of rising limb and receding limb are obtained by regression analysis as indicated 
below. 



(for rising limb) 


(4.5a) 


Ql^ Qin 
QwP ~ Qin ^ 



(for receding limb) 


(4.5b) 


The magnitude of peak discharge {Qwi^ and its corresponding time (r„p) are related with 
embankment geometry and inflow characteristics in the section where the combined analysis 
of cohesionless fine sand and cohesive soil embankment is carried out. 


vi. Energy Loss During the Breach Process 

Energy is computed at upstream of the embankment level and at downstream end of the 
embankment as the sum of datum energy head, potential head (depth of water) and the 
velocity head. The difference in energy between these two sections is considered to be the 
energy loss during the breach process. This energy loss nondimensionalised with upstream 
energy (E|) at a given instant of drop of water {AH), is plotted as showm in Figure 4.14(a). 


Here the drop {AH) is defined as the difference in water level in the reservoir to the 
downstream end of the embankment. This drop is non-dimensionalised with dam base length 
{Ids). The relation between them is obtained from regression analysis as: 



0.66 


(4.6a) 


A good correlation exists between them. To predict drop {AH ) during the breach, AH is 
nondimensionalised with H^ and is plotted against //% as shown in Figure 4.14(b). It may be 
noted that drop AH decreases with increase in t/% and the relation between them is given as: 
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c. Breach Erosion 

Variation of sediment bed profile, crest height and its width during breach process are 
analysed. Eroded sediment volume and its transport relation are analysed further. 

i. Variation of Sediment Bed Profile with Time 

As the water flows over the embankment, sediment gets eroded. Magnitude of erosion 
depends upon the shear stress acting on the sediment bed. It is observed that the sediment 
erosion starts from the downstream end of the embankment, recedes upstream within a short 
time. Once it reaches upstream edge of the embankment, the sediment gets eroded over the 
whole surfiice along the breach. These surfaces are measured along the centerline of the 
breach with time and they are plotted in Figures 4.15(a) and (b). It is observed that in the case 
of a dam of larger height the downstream bed profile gradually flattens with increase in time. 
In the case of lower height dam sediment erosion profile is slightly flatter compared to dams 
of larger height. These profiles are studied by nondimensionalising sediment bed level (/?<-) by 
its maximum level at the crest, which is usually located at the upstream top edge of the 
dam. This distance(X) along the embankment and along the flow is nondimensionalised by 
the value where }f.,„ occurs. The nondimensional plots of against X/X^ are plotted 
in Figures 4.16(a) and (b). It may be observed from nondimensional variation that for higher 
embankment dam model the variation for Hq 19.5 cm to 35 cm is in one pattern. For lower 
dam models of //;; equal to 14.5 cm the variation is in a slightly different way. However in 
both the cases, there is a steep rise to the maximum value, then gradual decrease with 
increasing X values. The relation between them obtained from regression analysis is a 
polynomial curve and is presented in Figures 4.16(a) and (b). 

ii. Variation of Crest Height ( hc „, ) 

The sediment bed profile variation with time is of similar pattern as shown in Figures 4.16(a) 
and (b) for all the cohesionless embankment models. With this in view, the height of 
sediment crest is non-dimensionalised with final depth of erosion (£)e,.)and plotted against the 
nondimensional time t/t^Q in Figure 4.17(a) . Here is the time at which 50% of the 
maximum erosion depth (Dg,) has taken place. It was observed, the relation between 
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Fig. 4.16(o); Nondimensional variation of sediment bed 
height (h^) at any location of breach with 
its nondimensional distance, X from heel 
ofthe dam (for dams of height=1 9.5cm, 
30 cm and 35cm) 



Fig. 


4 16fb)' Nondimensional variation of sediment 

height (h^) at any location of breach wi 
ts nondimensional distance, X I'rorn nee 
the dam (for dams of height-1 4.5cm 






Fig. 4.17(a): Nondimensional Variation of crest of erodible sediment bed 
(hem) at any time of breach with its nondimensional time 
of occurrence. 
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nondimensional depth of erosion with nondimensional time follows simple exponential decay 
curve as: 


Where, Y = 

^ “ ^^^ 50 . 2 nd Dg^= final depth of erosion 

The time scale Uq is studied by nondimensionalising it as . This value is plotted 

against nondimensional Q.„ as shown in Figure 4.17(b). The functional relation among them is 
as follow: 


hJg'Ho =-3194.91 + 160087Y-9.289x 10’Y- (4.7b) 

( This equation is not valid for X=0) 

Where X = 

From the Figure 4.1 7(a) it will be observed that the crest height can be predicted with time of 
erosion which is used in predicting sediment bed profiles along the breach from Figures 
4.16(a) and (b). The final breach depth of erosion at crest (D^^) is related with peak outflow 
discharge {Q\yp) and height of the dam {Hf) in a non-dimensional form as shown in Figure 
4.17(c). The functional relation developed for Der with these variables may be written as: 

- 0.043 

(4.7c) 

Depth of erosion is also correlated with inflow discharge parameter, defined as 
Qj['hXgfT) as shown in Figure 4.17(d). Using regression analysis, relation between 
them is written as: 


D.. 


Hr 


052 


a 


WP 




D, 


Hr 


0.71 


Qi. 


0.015 


(4.7d) 


For prediction of either Figure 4.17(c) or Figure 4.17(d) can be used depending upon the 
availability of the information. 



Uo^Cg/Ho) 


t5o(g/Ho)“-®=-31 94.91 4-1600! 
Here X = (Q,X/(Ao(gHo)°-®) 
Not Valid for X=0 


(Qln)ov/(Ao*(g 


g. 4.17(b): Variation of time see 
form with Qapn numb 


(D.r/Ho)=(Qwp/(gHD=)°-") 


5 \ 0 ^ S - 0.043 


^^.000 0.002 0.004 0.006 0.008 0.010 0.012 

Qwp/CqHd ) 

Fig. 4.17(c): Variation of final depth of erosion 
at crest, (Der) in nondimensional 
form with peak breach discharge(Qwp) 
and height of the dam(HD) 
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The relation between the depth of erosion and drop of water level at crest is also 
presented as an alternative approach to predict crest height variation. Here depth of erosion 
(Il^-u is non-dimensionalised with height of the dam (H^) and plotted against 
nondimensional water level drop at the crest as shown in Figure 4.17(e). Here is the crest 
height. It may be obserwcd from this figure that there is a good agreement between sediment 
crest level and water level at the crest with respect to the water level in the reservoir. The 
relation between them is obtained from regression analysis as: 


^hm) _ ^ Drop of water level at cre^tV 

~l in ; 


(4.7e) 


From the Equation (4.7e), by knowing the water surface at the crest one can predict, crest 
level of the dam at any time of breach. 



Fig, 4.17(d): Variation of final depth of erosion at crest, 
(Der) in nondimensional form with overage 
inflow discahrge to the reservoir, height (Hq) 
and cross-sectional area (AD)of the dam 
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iii. Rate of Crest Height Variation 


The rate of crest variation idh^„,/dt) is plotted against crest drop flow velocity This 

drop-flow velocity is computed based on the difference in water level upstream to the crest. 
The {dh,,ydt) is plotted against i^2gh, ) as shown in Figure 4.18(a). There are two distinct 

limbs, rising and receding limb. During the rising period, the rate of crest erosion is less 
compare to the receding period. During the transition state from rising to receding limb, rate 
of crest erosion has multiple values. 

Crest sediment erosion rate is also plotted with crest flow velocity (V^,^ ) in the Figure 
4.18(b). Crest erosion is a single valued function of The crest flow velocity is also 
associated with peak discharge Qnp and height of the dam . The magnitude of is 
computed by knowing discharge over the crest 2f,-and wetted flow area at the crest as: 

^ — (4-8a) 

Wetted flow area at the crest 


where Wetted flow area at the crest = B(,,. 

This was studied particularly to avoid the multiple values as indicated in the Figure 4.18(a). 
The relation between the crest erosion rate with flow velocity obtained from regression 
analysis is written as: 


dt 


(a, 


1.2 


1 . 5 * 10 *^ 


(4.8b) 


It may be noted in Figure 4.18(b) that the crest flow velocity is in the same unit of 
measurement as that of the rate of crest erosion. These two figures involves flow velocity 

(Fy, crest flow' velocity ), where one is obtained from breach discharge and breach 

flow area at the crest and the other one from the difference in water level upstream and at the 
crest.The prediction of is easier compared to prediction of V^,^. However the sediment 
crest erosion rate related to hii is a multiple valued function. Whereas the sediment crest 
erosion rate is a single valued function of . It may be advocated to use single valued 
functional relationship. The rate of sediment crest erosion is subjected to further analysis 
because of its importance in the prediction of the discharge through the breach and also to 
study in the dimensionless form. 




Fig. 4.18(b): Variation of rate of crest height, (dhcm/dt) 
with crest flow velocity, (V«,c) 
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The rate of decrease of sediment crest is nondimensionalised with average shear velocity and 

Qtyp ! yjs^D plotted against nondimensionalised drop in water level at the crest with 

respect to height of the dam Hq. These values are plotted for all the embankment models as 
shown in the Figure 4. 1 8(c). It may be noted that the rate of erosion decreases with increase 
in drop water level at the crest. The relation between them is as follows: 


1 


V, dt 




- 3.14 S- 




(4.8c) 


iv. Variation of Sediment Bed at the Downstream Top Edge of the Embankment 

Top width in some of the embankment model varies 3 to 10 times the height of the dam. Due 
to longer top width, the sediment bed level variation at this section also becomes important in 
the sediment erosion analysis. The bed level of the eroded sediment at the downstream edge 
of the embankment is studied with respect to drop at the downstream edge and both the 
quantities are nondimensionalised with Hp as shown in Figure 4.18(d). Initial erosion is not 
accounted in this figure. There appears to be a good agreement between the height of 
sediment bed, and the drop in water level at that section. The relation is obtained through 
regression analysis as follows: 

-^ = 1.47xe (4.9) 

Where Drop' = H ^ +K,) 

Where is the sediment bed level at the downstream top edge of the dam and Drop' is the 
drop in water level with respect to top level of the dam. 

v. Lateral Erosion at the Crest 

Lateral erosion of the crest is considered at the upstream top edge of the embankment. 
Erosion of the breach section at the crest results due to the reduction in the sediment crest 
level and widening of breach bank simultaneously. It is observed from the experiments that 
the breach widens faster than the crest height erosion. This can be seen clearly from Figure 
4.19(a) that the ratio of the rate of crest erosion to the rate of crest width erosion is less than 
unity from the starting of the breach process. This ratio remains fairly constant for longer part 
of the duration. Later in the final stage the rate of lateral erosion is much faster. 



((dh^/dt)/V.)/(Qwp/(gHo^)° 
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Drop in water level at crest/Ho 


Fig. 4.18(c): Variation of rate of nondimensional rate of 
crest height, (dh^,^/dt) with nondim ensional 
drop of water at that section 



(HD-(h-c+h',c))/Ho 


Fig. 4.18(d): Variation of sediment bed height at 
downstream edge of the dam with 
drop of water level at that section in 
a nondimensional form 
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In order to study the breach widening at crest at any time of breach, the variation of 
breach width with drop of water level at the crest from the reservoir level are plotted by non- 
diinensionalising both the quantities with It is observed that the breach crest width erodes 
with increase in drop as shown in Figure 4.19(b) the relation between them is given as: 

0.57 

(4.10a) 




Hr 


Hd~K-K 


Hr 


X 2.5 


ro ^ 

Where , {H^ -h„, - is called the drop of water level at the crest. This equation is 
obtained from regression analysis. The final magnitude of the breach width is also plotted 
against the maximum outflow discharge of the breach hydrograph in dimensionless form as 
indicated in the Figure 4.1 9(c). The relation between them is given as: 


(K)f 

Hr. 


' 0, 


- wp 


X 4.02 


(4.10b) 


In order to include the effect of the downstream water level on the erosion process of the 
crest width the plot of B^.,. /Hi^ against AH/Hd is carried out as shown in the Figure 4.20, 
where Alf. is water level difference between the crest and downstream end of the embankment 
i.e., toe of the embankment. It may be observed that /H^ decreases as AH/Hj:) increases. 
The relation between them is written as per regression analysis: 


B. 


Hr 


-=712.02e 


-i,9.ir 


(4.10c) 


Here 



f \ 

Qwp 




1.77 


Crest width (B^.,. ) being a major breach morphological variable it is related with both drop in 
water level at crest in Equation (4.10a) and drop AH^ .These two equations indicate B^.^ 
increases with time. It increases exponentially with AH^ whereas it varies in a power law 
form with the drop in water level at crest. But Equation (4.10c) is preferable. 



Bcr/Ho ■ (dh^/dt)/(dB„/dt) 
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4.19(a): Comparison of rate of crest height erosion to its 
lateral erosion with time in a nondimensional form 



Fig. 4.19(b): Variation of nondimensional breach crest 
width at any time of breach, (Bcr) with 
nondimensional drop of water at crest. 
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((Bcr),)/Ho)=X°-^%4.02 
Where. X=(Qwp/(gHDy=^ 


0 0.01 0.02 


Qwe/(gHV)°-^ 


Fig. 4.19(c); Variation of nondinnensional final breach 
crest width, (Ber)f with Qp^ number 


Y=71 2 02’^‘e" 

Where ■x=(AHe/Ho)/(Qwp/(gHo")°-")’-^ 


(AHc/HD)/((Qwp/(gHD j) 


Fig. 4.20: Variation of nondimensional breach crest 
width at any time of breach with drop in 
water level between crest and toe of the 
dam, AHc and Qpn 
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vi. Rate of Lateral Erosion at Crest 


The rate of width erosion with time is equally important in the study of erosion proeess along 
with the rate of vertical erosion. Widening of breach width at the crest causes the change in 
magnitude of How intensity at the crest. Keeping this in view, the rate of width erosion with 
time is plotted in the Figure 4.21(a). From this figure it maybe observed that the rate of width 
erosion decreases with time. To represent this variation in non-dimensional form, dB„ / dt is 

nondimensionalised with respect to ^2^ and time (t) with % The t,, is defined as the 

time at which 95% of the sediment erosion occurs in the breach. This is plotted in the Figure 
4.21(b). Using regression analysis the relation is established as: 


Idt) 





xL.HxlO"' 


(4.11a) 


The rale of erosion of the width at the crest is mainly a function of the magnitude of the 


crest flow velocity (U,,,). To test this postulation, rate of erosion is again plotted against the 
crest flow velocity in dimensionless form. The relationship exists among them is indicated in 


the Figure 4.21(c). The relation obtained among them is as follows: 


bik, 

\ dt 


/ \ I 62 

U / > 




UT 


SR 




X 0.0128 


(4.11b) 



(dB„/dt)/(2gHof* ^ ^ dB„/dt^ 





of Dam Breach Process 1 05 



nondimensional crest flow velocity 
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vii. Variation of Width at Downstream Top Edge of the Embankment 


The breach width at downstream top edge of the embankment is denoted as 5;. Its 
qualitative variation is presented in the section containing the description of breach 
morphology. Its variation is studied using similar parameters that are used for the crest width 
study, Bcr ■ The non-dimensional variation of B,/ with is plotted against AH/H^ in the 
Figure 4.22(a). The relation is given as follow: 

^= 164 .. ( 4 . 11 ,) 

It may be observed that there is a sudden increase in the width at initial stages of the erosion; 
as erosion proceeds, width remains almost constant. 

Rate of erosion of downstream width is plotted against the flow velocity along 
withy^"j, as shown in Figure 4.22(b). The results indicate two distinct variations, one 


following the rising limb and another one with a falling limb. The rate of erosion of width is 
almost linear in both the states of flow, but however, there is a wide gap between both of them 
indicating the c.xistence of transitional state of flow. 

The relative variation of crest widths at the downstream top edge and the upstream top 


B' . t 

edge denoted as -- is plotted against — as shown in Figure 4.22(c). The functional relation 


B.. 


‘95 


among them is: 


B_^ 

~B... 


t 

r r V 

f 0 

0.067+ 3.65 3.95 X 

— +1.4x 

^95 




, Valid upto—<\.0 (4.1 Id) 


‘95 


The regression curve is drawn for the data. It may be observed that this ratio becomes greater 


than 1 .0, when -- is 0.4 and more. This indicates that in the beginning the crest width 

V)ii 

erosion at the upstream edge is more compared to that at the downstream edge of the 
embankment. As the breach erosion continues the crest width at downstream end enlarges 
gradually and becomes equal to that at the upstream edge at t = 0.4 t 95 . Later on the crest 
width at downstream edge widens and remains fairly constant after r/% =_0.8. Further 
onwards its magnitude remains constant at an average value of 1.2. 



d(B'„)/dt 


KR erosion a 
nondimensio 
nstream of t 


° Rising Limb , Y=X ' ^0.00016 

A.„A,A Falling Limb. Y=x’-°^^^0.00094 

★ * ♦ ★ ★ Transition 


V'^/(Qwp/(gHD'’)°’'^“) (cm/sec) 


Fig. 4.22(b): Variation of rate of breach width at 

downstream top edge of the dam with 
the flow velocity at that section 






nondimensional time 
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viii. Variation of Sediment Volume 

The volume of sediment eroded is computed, knowing the cross - sectional area of sediment 
eroded at different sections along the breach sections where measurements were taken. The 
typical variation of cumulative volume of sediment eroded with time is shown in Figure 4 5 
The rate of cumulative volume of sediment at a given time is considered as the sediment 
discharge (&). It is observed that the sediment erosion is a function of the magnitude of the 
flow discharge {Qn). Based on this consideration, the cumulative volume of sediment is 
non-dimensionalised with cumulative outflow volume of water (IVJ at any time and is 
plotted against t/tg; as shown in Figure 4.23(a). The functional relation among them may be 
written from regression analysis as: 


V T/ 

y v 

Mm* r 11/ 


= 0.03 X 






(4.12a) 


Here % is the time at which 95% of the cumulative sediment eroded. The non-dimensional 
value of varies almost inversely with non-dimensional time variation as indicated in 


Figure 4.23(a). Instead of taking Jl/y, the discharge ratio of sediment and water is considered 
for further analysis. 

ix. Variation of Sediment Discharge Intensity 

In order to study sediment discharge intensity ) with flow discharge intensity (^,„ ), plot 
of qs versus is carried out in Figure 4.23(b). Here, sediment discharge intensity (q^) is 
obtained from ^^as q^= Q^/ average breach width . The magnitude of q^ is plotted against the 
water discharge intensity r/,,, as shown in Figure 4.23(b). Here the water discharge intensity is 
obtained as q,^ ~ Qifu'Avcrage breach width. The average breach width is the average of the 
breach width at the top edge, middle of the top section and downstream top edge of the 
embankment. Using the regression analysis, the relation between them is as follows: 

=9.86x1 (4.12b) 
It may be noted that varies as qj'^^, where exponent 1 .43 is much below what is observed 
in the sand bed stream. 
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Fig. 4.23(a); Variation of nondinnensionalised cumulative eroded 
sediment volume by cumulative flow of water in 
breach section with nondimensionalised time. 



Fig. 4.23(b): Variation of sediment discharge intensity, (qs) 
with breach water discharge intensity, (q*) 
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Bed shear stress is considered to be important parameter for the sediment discharge. 
Sediment discharge intensity (q,) is plotted against square of shear velocity (V.^) in Figure 
4.23(c), where V, ,T(, - bed shear stress, /?„, = density of water. The relation 

between q, and fv’ indicates a well defined relation for the value of F,= 7.0 cm/sec and above 
in a power law form, q, oc( vJ) ' as given below, 

(4.12c) 

On the same plot, polynomial variation between q^ and V/ is plotted. This relation is shown 
in the following equation. 

r/, = - 0.0456 + 0.001 - 1.25 x 1 0'' x (4.1 2d) 

where X - (F; ) 

In Figure 4.23(c), it is clearly shown that there is a deviation from the power law relation for 
F. < l.Qem / sec . This situation occurs at the final stage of breach process. 

X. Variation of Sediment Concentration 


The ratio Qs'^Qw can also be termed as the sediment concentration. By plotting with 

non-dimensional variation in water surface level between upstream reservoir level and 
downstream end of the embankment AH, is shown in Figure 4.23(d). Sediment concentration 
is high during the initial phase of the erosion and gradually reduces with decrease in AH/Hq 
values. The maximum concentration reaches around ten percent (10%), which is considered 
to be hyper concentrated level. The magnitude of hyper-concentration is 8 percent and above. 
The concentration study involves for the whole breach section. The relation obtained through 
regression analysis is as follows; 

£l. = _ 0.005 1 -t- 0.14526 - 0.3826' -K 0.34526' 

a 


Here 26 = 


AH 

Ho 


(4.12e) 
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xi. Sediment Load Computation 

The product of qs with bulk density of sediment (/*) of the embankment soil is considered as 
the sediment load. It is non-dimcnsionalised with jq , d^Q and where unit weight of 
sediment, d^o = median size of particle, = sediment resisting velocity, defined as: 

^^*^50 *tan^+ — 1 (4.12f) 

y \ Piy ^PivJJ 

where, 

Pir V Piy / 

^ = angle of repose 

c = cohesive stress, obtained from triaxial undrained shear test 
p ^ - density of sediment particle 

The dimensionless sediment load is plotted against non dimensional shear velocity, V,' / 
in Figure 4.24. The relation between them is obtained from regression analysis and is given as: 

/ , \ 1.938 

g;= 0,0082.^ (4.12g) 

y^sRJ 

It can be observed that for smaller values of Vp IV^r, there is a considerable scatter in the 


sediment load. 
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4.4 Analysis of Cohesive Soil Embankment 

In the previous section, the analysis of cohesionless soil embankment has been reported. It 
may be noted that the variation of flow parameters, geometric parameters and sediment 
parameters are interrelated to bring an insight into the breach processes of cohesionless soil 
embankment. In the present section the analysis of cohesive soil embankment is carried out in 
the light ofcohesionlcss soil approach. 

The basic data comprises of variation of water level, variation of sediment bed level and 
the corresponding breach flow width at different time along the longitudinal section of the 
breach which are presented in Chapter III, Figures 3.4(a) to 3.4(m). These figures contain the 
data for the run number 0404 for the sections 1 to 8. Only water levels are presented for the 
sections 1, 2 and 3 which are in the reservoir, upstream of the embankment. For other sections 
all the three data arc presented. Using these basic information further analysis on breach 
morphology, breach hydraulics and breach erosion are carried out. 

a. Breach Morphology 

The variations of breach planfomt, bed and flow profiles along the flow in the breach and 
breach cross-sections arc analysed. 

i. Breach Planform Variations 

During the breach erosion, it was observ'ed that the crest of the erodible bed profile is moving 
from the downstream towards the upstream as the erosion progress. Also widening of breach 
section follows a particular pattern during the breach process. In order to study the variation 
of depth of flow at crest and its corresponding width variation during the breach process, four 
situations namely two situations before the occurrence of breach peak outflow ( ), one at 

the time of and other after peak outflow are shown in the Figures 4.25(a), (b), (c), and (d) 
respectively. It may be observed that the position of crest is very near the upstream edge of 
the embankment. The width is contracting slowly from upstream to the downstream edge. As 
the discharge increases the crest moves upstream and reaches upstream face of the 
embankment before the attainment of peak flow. By the time peak flow is attained, the 
exposure of the hard bed occurred due to erosion. After the peak flow occurrence, width 
reduces upto the upstream top edge of the dam and remains fairly constant further onward. 



Width (cm) Wl . BL Width (cnn) WL . BL {cm) With fcm) W1 . 01 (cm) Width (cm) Wi . Bl (cm) 
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Water level 
Bed level 


Run No.: 0404 (t=500sec) 
Before Breach ^Peak Discharge 
Q,« 1 1 48.5 cm /sec 




Run Nc : 0404(t=15503ec) 
At Peak Disch^arge 
0^*2653 cmVsec 


f ig 4,25c 



301 


fS »1D0 -75 -50 -25 

lofMgitwdjnol Drstorn:® from wp*lr«im edge of CmborWtment (cm) 


Run No : 0404(t« lOOOsec) 
After Breach Pcok Discharge 
CU- 1209 cm /sec 


f'tg 4.254 

Fiqs.4.25o. b, c and d: Voriations of Water levels, sediment Bed Levels 

ond Breoch Width along the Longitudinal direction of Dam breach 
at different time for Cohesive Soil embankment model 
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The constriction of the width causes accelerated flow as indicated in Figures 4.25(a), (b), and 
(c). 

ii. Variations of Bed and Flow Profiles Along the Flow in the Breach 

The variation of water level and width variation at sections, the downstream edge (section - 6) 

and section in between downstream edge and downstream end, denoted as section - 7 follows 
slightly different variation compared to other. Water level and bed level fall very rapidly in 
these two sections compared to other sections like sections 4, 5 and 6. This is due to steep 
erosion in the bed causing a valley like flow. This type of channel fomiation starts from 
downstream end as can be seen from the Figures 4.25(a), (b), (c), and (d) which indicates 
around 500 seconds at section - 7 and at 1000 seconds at section - 6 for the given mn. The 
width variation at sections 6 and 7 is more like a rectangular section. At the section - 5 width 
variation increases rapidly and later decreases gradually after the peak outflow discharge. 
The shape of the cross-section remains trapezoidal during the erosion process. Figures 4.25(a), 
(b) and (c) indicates free weir flow situation in which gradual drop in water surface is upto the 
crest and from the crest downstream the steep reduction of water surface profile occurs. 
Almost same water surface level between upstream and downstream indicate submerged weir 
type flow for F'igure 4.25 (d). Slope of the downstream eroded sediment surface is steeper 
initially and as the breach process continues its magnitude decreases. 

iii. Breach Cross-Sectional Variation 

The flow cross sectional profiles for sections 5, 6 and 7 are plotted for run no.0405 in the 
Figure 4.26. Here section 5 indicates section at upstream top edge of the embankment, section 
6 is at downstream top edge of the embankment and section 7 is nearer to the downstream end 
of the embankment. Water surface level is shown with dotted line and solid line for the bed 
level. When time (t) is equal to 500 seconds from start of overtopping of flow, cross sections 
at 5 and 6 arc more or less trapezoidal in nature, but in section 7 it is almost like a rectangular 
narrow valley. At time (/) equal to 1000 seconds, the cross sectional shape at section 5 
continues as trapezoidal but at section 6 it becomes almost rectangular. This trend in variation 
in cross sectional shapes at sections 5, 6 and 7 continues at the time of peak flow and after 
peak flow as indicated for t=- 1550 seconds and t = 1 900 seconds respectively. It may be 
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observed in the beginning of erosion that at the downstream end a narrow steep rectangular 
shape of cross-sectional erosion takes place but at the upstream end it is more like a 
trapezoidal shape of erosion . Later on as the erosion process continues, the widening of 
downstream section occurs and it becomes more like trapezoidal shape. However, in the 
upstream edge, section 5 the cross sectional shape remains almost same with a rnagnimde of 
side slope remaining constant. At the section 6 after initiation of breach, the cross sectional 
shape becomes more rectangular in shape. 

b. Breach Hydraulics 

Computation of discharge over the breach section is carried out on the line of analysis of 
cohesionless embankment by plotting discharge intensity (gj against head acting over the 
crest (/ij as shown in Figure 4.27(a). Using regression analysis, the relation between them is 
written as; 

r/^,=5.45x;i'-^‘-^ (4.13a) 

The exponent of /i„ is nearly equal to the discharge over the spillway. There is a wide scatter 
observed in the relation. With this scatter in view, is nondimensionalised with respect to 

2 

head acting h,, as - c. , which is a function of discharge coefficient for weir flow, 

(q). This is plotted against the ratio of depth of flow acting on the crest to the height of the 
crest from the bed level as indicated in Figure 4.27(b). It may be observed that / {Igh^Y' 
decreases with increase in the value of {h^ I ). This decrease is rapid upto ) < 

0.4. For higher values it remains fairly constant. As observed in the discharge analysis of 
cohesionless soil embankment, free weir flow occurs at the crest for ^ 0.4 ; for 

higher values the flow becomes submerged weir flow and later on open channel flow. The 
steep drop in the values of / {2gh^y^ fox {h^! ) < 0.4 indicates the effect of crest 

erosion, even though discharge relationship behaves similar to weir flow. For the comparison 
of the variation of with I ) the relationship obtained for cohesionless soil 
embankment is incorporated in the form of dotted line. Observing the scatter in the discharge 
relationship data, it was thought, a statistical representation of the deviation is needed. With 
this in view the statistical deviation from mean value is plotted against the frequency of 
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ccurrence as shown in the Figure 4.27(c). The high frequency are observed around the 
deviation of the magnitude of ± 0.04, having two peaks. 

c. Breach Erosion 

The vanation of rate of crest height, width at crest , rate of width erosion and rate of cross 
sectional area at crest are presented for cohesive soil embankment breach, 
i. Rate of Crest Variation 

Rate of crest variation is represented as dlfjj/dt. It is nondimensionalised with upstream flow 

velocity, ^2gh,, along with and is plotted against r/% as shown in Figure 4.28(a). The 
relation between them is given by polynomial function as: 


=9.8x10-' -6.4xl0-'f— 1+2.3x10-'^ 

\ ^95 / V 

where , = [{“’'-/i) / 

It can be seen from Figure 4.28(a) that there is a steep reduction in the rate of crest erosion 
towards the end of breach process. 

Rate of crest variation is also related with drop in water level between upstream and 
downstream of embankment (AH). The nondimensionalised rate of crest erosion 

(dKJ‘it)if-gK /(&»)' is plotted against nondimensional drop in water levels AH /Hq 

as shown in Figure 4.28(b) .The functional relation among them can be written as: 




Two curves, one from power law relation and another from polynomial fit are obtained . The 
polynomial functional relation appears to be better representation of the data. From Figure 
4.28(a) and (b) the rate of variation of crest height of sediment bed can be obtained in any 


time during the breach for the corresponding drop in water level {AH) . 
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ii. Variation of Crest Width 

Qualitative description on the variation of crest width is described under the breach 
morphology section. It is necessary to know quantitative variation of width at the crest 
position also. Breach width at crest (B,,) is nondimensionalised with height of the 
embankment (74) and is plotted against the drop of water level, at the crest which is also 
nondimensionalised with 77^. It may be observed that the breach width attains its maximum, 
almost equal to the 1 .9 times TTq and decreases linearly with the drop in water level crest This 
observation can be seen in the Figure 4.29(a). The functional relation among them is, 

= 1-896 -0.058Z (4.15a) 


Where X = [Drop of water level with respect to H^)! 

The relation between the breach width {B,,) at the crest (section 5) and the breach width (5.) 
at section 7 located near the downstream end of the embankment are studied. The breach 
width ratio (Bj / 74) is plotted with in Figure 4.29(b). The functional relation among 
them is obtained from regression analysis as: 

\ 0.44 

=O.U2-26.02.ir+17992.9X’ - 1979090 /r'+ 6408 xlO’A" (4.15b) 

I gill J 



where X 



During initial phase of erosion, {By /^r) is small, as erosion continues this ratio of width 
attains high value (0.3) at X = 0.07 and there after remain constant. This can also be observed 
in Figures 4.25(a), (b), (c) and (d). 



(B7/Bj((Q«./(gH 
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iii. Rate of Width Erosion of Cohesive Soil Embankment Model 
The rate of width erosion at section 5, the upstream edge of the embankment top is analysed 
by plotting dB,,/dt against time (t) as shown in Figure 4.30(a). It may be observed that the 
rate of width erosion increases very steeply within a short period and remains fairly constant 
with magnitude zero throughout the breach process. This indicates that the dimension of the 
breach width attains maximum in the initial process of breach erosion and remains althrough 
constant during the remaining part of the breach erosion. Nondimensionalised rate of width of 
erosion is plotted against non-dimensionalised time, t/tg^ as shown in Figure 4.30(b) The 
dimensionless rate of width erosion as shown in Figure 4.30(b) attains peak very soon after 
initiation of the breach process and remains zero althrough. It is found in the breach process 
that there are few cases of sudden fluctuation in rate of breach width erosion .This may be due 
to collapse of the banks. 



Fig. 4.30(a): Variation of rate of width erosion at crest 
for dam model of cohesive soil with time 
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iv. Rate of Variation of Cross-Sectional Area at Crest (A^,) 

In the last subsection the rate of depth and width erosion at crest are analysed. Taking 
together, the depth and width of erosion into account, the rate of variation of eroded crest 
cross-sectional area is analysed with drop in water level between the upstream and 
downstream of the dam. This is carried out in nondimensional way as shown in Figure 4.31. 
The functional relation among them is as follows. 

^„ = 0.00196-0.01366.Y + 0.0419Z--0.03657Z^+0.0066Z^ (4.16) 

where = dt) ! ^2ghJ ! and Z = zlFf/Ho 

It may be obsen^ed that the rate of erosion of crest cross-sectional area increases and attains a 
peak near dU/Hp = 0.7 and then decreases upto AH/Hq s 0.3. Before discussing the 
combined analysis of breach process for the embankments consisting of cohesionless and 
cohesive soil separately, the variation of the variables like head acting over the crest (/zj, 
height of crest (//„„), width of crest {B,,,), the rate of erosion of crest and rate of width erosion 
at crest arc presented for both cohesionless and cohesive soil dam model together for 
comparison. 

d. Comparative Study Between Cohesionless and Cohesive Soil Embankment Breach 
Observations 

The presence of cohesion in the soil adds more resistance to the erosion. This increase in 
resistance causes the decrease in the rate of crest erosion and rate of width erosion compared 
to cohesionless soil. With this observation, the head acting over the crest (hj, height of crest 
(h^J and width of crest fB^^) are nondimensionalised with height of the dam (H^) and plotted 
against as shown in Figures 4.32(a), (b) and (c), respectively where t is time at any 
instant from the beginning of the erosion and tpj is the time by which 95% of the breach 
erosion has taken place. These data are chosen for nearly same height of the embankment 
dam, however, there is variation in the magnitude of inflow to the reservoir. Still these plots 
give a qualitative comparison in behaviour of breach process with cohesive and cohesionless 
soil embankment dam models. 

The Figure 4.32(a) indicates the variation of with t/ tg^ for the embankment 

models built using cohesionless soil (type-A soil) and cohesive soil (type-B soil). The h^H[) 
for cohesionless soil decreases from start of the overtopping of the flow, whereas in cohesive 



(dBer/dt)/(2gHD) 
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Fig. 4.30(b); Variation of nondimensionai rate of breach 
width with nondimensionai time for cohesive 
soil embankment dam 



Fig. 4.31: Variation of nondimensionai eroded crest area (Agr) 
with nondimensionai drop in water level (AH/Hd) 
betwen upstream and downstream of the dam. 
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Fig. 4.32(a): Comparison of the variations of crest level of 
erodible sediment bed for the breach model of 
cohesive and cohesionless soil embankment dams 



Fig. 4.32(b): Comparison of the variations of head of water level 
over the crest of erodible sediment bed for breach 
model of cohesionless and cohesive soill 
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Fig. 4.32(c): Comparison of the variations of crest width of 

erodible sediment bed for the dam breach model 
of cohesionless and cohesive soil embankment dams 
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soil it starts slightly later. However, their decreasing relation as observed is fairly parallel to 
each other indicating the crest erosion is similar in nature. 

The variations of the head acting over crest {hf) for both the soil are shown in Figure 
4.32(b). The maximum magnitude of occurs in early phase of the breaching for 

cohesionless soil (type-A soil) compared to cohesive soil (type-B soil) and also the magnitude 
of the peak values are different. 

The comparison of width variation for both the soil is shown in the Figure 4.32(c). It may be 
observed that width enlarges faster initially, attains peak and starts decreasing then slowly 
increases and reaches constant as the end of the breach process approaches for cohesive soil 
(type-B soil). The increase crest width for cohesionless soil (type-A soil) occurs during ?/% < 
0.15, reaches a peak with a magnitude of BJHq = 0.8, then decreases with increase in t/tg^ 
upto 0.55 and remains fairly constant thereafter. The magnitude of peak of the breach width 
for cohesive soil is of the order of three times whereas in the case of cohesionless soil it 
is of the order of for this particular chosen run. The sudden rise in the width in case of 
cohesive soil is an indication of collapse of the bank at that time. In the computation of the 
sediment volume eroded at any instant of time the average of width at that instant is 
considered, not the abrupt deviation in this value because it was very difficult to take into 
account the bank collapse and its erosion at any instant of time. 

The rate of erosion of the crest height and the width are non-dimensionalised with Hjg and 
are plotted in the Figures 4.33(a), (b) respectively. The magnitude of the non-dimensional 
crest erosion for cohesionless soil is much higher compared to its corresponding value for 
cohesive soil. The values of t/ tg^ for peak rate of crest erosion are almost equal for both 
cohesive and cohcsionlcss soil. The magnitude of non-dimensional rate of crest erosion for 
cohesionless soil is of the order of 0.012 whereas for cohesive soil is of the order of 0.00008. 

The variation of the rate of crest width indicates that most of the erosion occurs in the 
initial period of the breach in case of cohesionless soil (type-A soil), whereas in cohesive soil 
(type-B soil) it occurs in the later stages of the breach. The magnitude of its peak for 
cohesionless soil is of the order of = 0.00013, whereas for cohesive soil it is = 0.0044. 

It may be summarised that the magnitude of the rate of erosion of crest height and crest 
width is higher in cohesionless soil compared to cohesive soil embankment breach. This is 
due to extra resistance offered by cohesion in cohesive soil. 
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Fiq. 4.55{a): Comparison of the rate of variations of crest 
height of erodible sediment bed for the dam 
breach model of cohesionless(type-A) and 
cohesive (type-B) soil 



Fiq. 4.3j(b): Comparison of the rate of variations of crest 
vv'idth of erodible sediment bed for the dam 
breach model of cohesionless (type— A) and 
cohesive soil(type-B) 
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4.5 Analysis of Breach Process for Cohesive and Cohesionless Soil 
Embankments 

In the section 4.4d, comparison of the breach process was carried out for cohesionless and 
cohesive soil embankments. In the present section it is planned to analyse the combined data 
taking soil properties into consideration. In this section the analysis contain variation of water 
surface (//^„) in the reservoir, head acting on the crest (hj and crest height energy 

variation (AE), sediment concentration (Qs/Qi^) and sediment loadCg^ ) along with time 
scales. After this, the final values of breach width, final depth of breach erosion, duration of 
breach are analysed. The occun-ence of peak outflow discharge (Q^p) and its corresponding 
time are also analysed. From this combined analysis it is hoped to evolve a general 
method of prediction and evaluation of the breach process. 

a. Variation of Water Surface in (he Reservoir 

Water surface level measured at a section 75 centimeters upstream of the embankment is 
plotted against time taken in the form (i-t,)/(i^yt,), where /» is the time at which maximum 
water level the reservoir occurred during the breach process. The ratio of h^./hy.^ is 

plotted against as shown in the Figure 4.34(a). It may be seen that the gradual 

reduction in water level as the duration of the breach process increases. In the final stage of 
erosion the water level in the reservoir attains certain magnitude after which the reduction in 
the reservoir level is almost negligible. This variation is mainly due to the intensity of inflow 
discharge (Q,,,) to the reservoir. If Qp, is large the erosion in the breach is higher and final 
stage of water in the reserv'oir is lower. This variation can be seen in Figure 4.34(a) for both 
the embankment breach processes. As is obvious from the Figure 4.34(a) there is wide 
distribution in with (t - 1.)/ - 1,). In order to fit a single curve to these 

scattered points, the modification of X-axis is considered. The value on X-axis ((3) for 
/i„ = 0.7 is read and subjected to regression analysis. From regression analysis the 
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relation between X intercept for =0.7 is correlated with oc- 

where, X = 

This relation is shown in Figure 4.34(b). Time scale t„ time at which maximum value ofh^, 

occurs from the start ofthe flow overtopping, is nondimensionalised with soil resisting 

velocity (Fsa) and plotted against number as shown in Figure 4.34(c). Polynomial 

functional relation is obtained from regression analysis and is written as : 

y = H 9503 + 3.04 x 1 O' jlT - 1.02 x 1 0® (4. 1 7b) 


Where Y 


and X 


Vsn 


a 


WP 


<^APN 


Using results of I'igurcs 4.34(b) and 4.34(c), the results presented in Figure 4.34(a) are 

Qpsn 


replottcd as (/z„.//j„.,,,„ ) against 


^ h)<, “ t, Jl, Fc0u ) 


in the Figure 4.34(d). The relation shown 


in Figure 4.34(d) shows a fairly good agreement and is represented in polynomial functional 
form as: 

h 


Where , a 


■ 1.0- 0.0289a -0.084«' + 0.01a' 

('-'■) JQ> 


(4.17c) 


fn 

iPSN 




b. Variation of Reservoir Water Level in relation to Downstream Water Level 

Ml is the difference in water level upstream of the embankment, in the reservoir and the 
downstream end of the breach channel. The variation of this quantity is non-dimensionalised 
with its maximum s aluc, denoted as AH- and is plotted against {t - 1—) /(tg^ - 1—) as shown in 
Figure 4.35(a) for all the runs, where t** corresponds to the time of occurrence of AH- . It may 
be observed tliat as the breach process continues the magnitude of AH decreases and reaches 
almost 0.1 at final stage. The values of X-axis (fl) in Figure 4.35(a) are obtained with intercept 
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X 


Ficj, Vnricition of /?. the X-intercept at y=0.7 in the 

fi(p 4.34a for different runs of experinnents on 
dom breach model of cohesive and cohesionless 
soil with X, a function of peak outflow 
discharge, soil resisting velocity and the 
height of dam 



Qwp/(Ao*(9Hd) ) 


Fiq. 4.34(c): Variation of reservoir time scale (t.) with 
peak breach outflow discharge, height of 
the dam and its cross sectional area. 
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Fig. 4.34(d): Variation of reservoir water level(hv,r) with time in a 

nondimensional form along with Vsrn and Qpsn numbers 
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of y-axis = 0.7, in Figure 4.35(a). This intercept is analysed further. It may be observed that 
the intercept varies with and number as indicated in the Figure 4.35(b) The relation 
of f2with <2/>.v,v and Fy/j/vf’btained from regression analysis is as follows: 

n = 2.02/:'’"''’ 


(4.18a) 


'1 


K = X V 

The reser\'oir drop time scale, f.., corresponding to the maximum difference of water level, 
(AH*), between upstream and downstream of the dam is plotted in the form ac^ainst 

^SR 

Qapn number us shown in I ii^inc 4.35(c). Xhc functiontil relation between them in the form of 
quadratic equation obtained from regression analysis may be written as: 

F = 84610.2 f 3.83 x lO' J 1.192 x 10' 

/ V . y. QWP 

whem Y = ~ and X - — 

K« AjgH^ 


(4.18b) 


Using correlated parameters for ,X-axis from Figures 4.35(b) and (c) the data shown in 


Figure 4.35(a) is replaced as against 

All. 


( t-t ^ 

f V ^ 

SRN 

V^95 

/-)0.67 

C idpSN / 


in Figure 4.35(d). The 


agreement between the data and cur\^c obtained from regression can be considered to be good. 
The relation between them is represented in polynomial functional form as: 

AH 


AH, 


where , i' 


0.968 + 0.056f - Oils- + 0.026f^ 


(4.18b) 


[ - t,. 

^ 6,5 - 1, 


* V ) 

\^PSN ' SRN ) 


0.469 


c. Variation of Crest Height 

The magnitude of water Hewing in the breach depends on the head acting over the crest QiJ. 
Magnitude of is the tlifference in water level in the reservoir and sediment bed level at 
crest. Variation ol' sediment bed level at the crest is important for further analysis. With this in 
view, crest level (/?,.„,), is measured from bed of the flume to the crest level of the sediment 
bed, which is nondimensionalised with height of dam (Ho) and plotted with (1-(AH/AHJ) 
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along with Qpsn ^ shown in Figure 4.36. Crest level decreases with increase in {1-(AH / 
AHJ). From the regression analysis, polynomial functional relation is obtained as: 

h.. 


H, 


where, x 


■■ 0.91 - 16 . 652 ' + 1902 ' - 868 . 42 ^ 


(4.19a) 


1 


A/-/ 

A//J 


♦ ;/« 0 
^ SRN 


The scatter in the data about the mean curve obtained from regression analysis is high. It 
may be noted that the variation of for water level rising from overtopping level to the 

maximum level on the embankment is not included in the analysis. By the time AH attains 
ail the erosion [u-ocess is already started and nearly 5% of embankment height must have 
eroded by that time. I'his observation can be seen in the Figure 4.36. 


d. Rate of Crest Height Erosion (<//«, 

The rate of crest erosion along with time is essential for the breach process study. Rate of 
crest erosion is nondimensionalised with velocity of flow approaching the crest, denoted as 


^2gh~ . Dimensionless rate of crest erosion is plotted against difference in water level 
between the upstream and ilownstream of the embankment and nondimensionalised with Lps, 
the base length of the embankment. These values are plotted individually for each run. Power 
law relation is fitted with regression analysis. The exponent and multiplying constant of the 
power law n,, imd K 22 ‘ire analysed with the ^SRN as indicated in Figures 4.37(a) and (b). 
Using related values of and K 22 , the rate of crest erosion is plotted again as shown in 
Figure 4.37(c). It may be observed that the scatter in data is high , but there is a general 

variation of increasing trend oi'dh^^dt with AH/ Lpg as: 


dt y 

75T 


'*22 


V LpgJ 


(4.19b) 


Incorporating the values of A:,, and 032 in Equation (4.19 b), the relation may be simplified as. 



= 0.1<^'^ 


(4.19c) 
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0^01 




0.5 


n22-0.0003*(VsRN)~ 


Where , n22 is the exponent 
of the power law relation 
between individual run of 
((dhem/cit)/(2gh,)^-^) 

Versus dH/Loa 


l-iq. 4,.')/((i); Vciriation of exponent njj with soil 
resisting velocity nunnber(VsRN) 



Fiq. 4..5 7(b): Variation of (< 22 . sediment erosion _ 
coefficient with soil resisting velocity 
number(VsRN) 
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Fig. 4.37(c): Variation of nondimensional rate of erosion of crest height(h 
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where , ^ ■■ 


A// 


N (0.0003 /(/”=) 


X X 1.8x10-" 




It may be observed here that the data of both types of the embankment merged together, 
can be seen from Figure 4.37(c). 


as 


e. Variation of Head of Water Over the Crest (/i„) 

The head of water over crest of sediment bed during dam breach is noted as /?„. This is the 
difference between reservoir lcvel(/;, J and crest of sediment bed profile(/!^J as: 

K = Kr 

In the analysis of cohcsionless and cohesive embankment separately, is non- 
dimensionalised with and plotted against / / 0 ,,,^, in Figure 4.38(a). The value of for 
breach discharge calculation can be obtained at any time of breach, if the scales h^pmd 
can be estimated. Attempt to correlate t,,,,/, for combined data was not successful. Instead, the 
analysis oUi„p based on reservoir level variation is attempted here. 

It is realised for the computation of/;,, using the breach variables like reservoir level(/ 2 „,,.) 
and crest of sediment bed( /;,.„,) will be more appropriate compared to time as shown in Figure 
4.38(a). Prior to this in sections 4.5a and 4.5b the method to predict h^,. and /;„„ is discussed. 
One can directly use liquations (4. 1 7c) and (4.19a) for /;„ computation. 

In another way of analysis /;„ is related to reservoir level (/;„,,. ) in a normalised form, for 
both the cohesive and cohesionless embankments data as: 

= X (4.20a) 


w'here, .v = 


ilKr„,-Kr) 

(iKm ~Kr-) 


, where is the reservoir water level 


corresponding to . 

In this equation to predict at any time of breach reservoir level to be known for the 
estimated exponent and scales like: ?], and h,^p . The variation of h^/h^ip is plotted 

for combined analysis of both the cohesive and cohesionless embankments with normalised 
reservoir water level (/;„,,.) as shown in Figure 4.38(b). Even though the trend in variation is 
same but large scatter can be seen in both rising and falling limbs. Falling limb data is further 
analysed using the relation given as equation (4.20a). The value of exponent {rf) of the 
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Equation {4.20a) varied from 0.06 to 0.68 as indicated in Figure 4.38(b). All the data fall in 
the range of exponent rj variation between 0.06 to 0.68. The is the reservoir level 
corresponding to top of the dam, (h^J is the reservoir level acting at any time under 
consideration, and is the reservoir level correspond to h^p.. 

1. The variation of i] is plotted against as indicated in Figure 4.38(c) using regression 
analysis, a power law relation is given as; 

7 ~ ^s«A/ ^ ^4.1 (4.20b) 

2. It was obsei-vcd in the experiments that the magnitude of is slightly higher than the 
height of the dam. Since the difference in their magnitude is not appreciable, an assumption of 

= Hp is considered as very close to the reality. 

3. The /i,,,,.. plotted against nondimensionalising with maximum water level in 

the reservoir. There is a good agreement between these two variables as can be seen from 
Figure 4.38(d). A polynomial relation is obtained based on regression analysis. 

= 1 0,49K„, - 40.5V, - 4039.4F,;,, (4.20c) 

Krn. 

4. The magnitude of h„p is the maximum value of //„ occurring during the erosion process. 
The value is needed to predict the variation of /i„ as indicated in Figure 4.38 (b). The 
magnitude of is computed as; 

fhp = °f h^„,,coiresponding to h^.^. ) (4.20d) 

The value of /j„„. is non-dimensionalised with reservoir level at that time is plotted 
against as shown in Figure 4.38(e). The functional relation among them is obtained by 
regression analysis as: 

ttsaiA = 0.443 (4.20e) 

\Kr^) 

The magnitude of head acting over the sediment crest of breach can be computed by 
knowing and from Figure. 4.38(b). The magnitude of h^p and rj can be 

found from the relations developed in Figures 4.38(c), (d) and (e). 
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Fig. 4.38(1)): Variation of normalised head of water over erodible crest 
of sediment bed with normalised reservoir water level 
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f. Variation of Energy Loss During the Breach Process 

The modes of energy dissipation during the breach process can be categorised in the 
following headings: 

a. Energy spent in eroding and transporting the sediment. 

b. Due to fonnation of hydraulic jump intcnnittently. 

c. Due to boundary friction. 

d. Water fall occuixing particularly in the cohesive embankment. 

It was difficult to compute energy loss under the above sub items. Combined energy loss 
is computed as the difference in energy between the upstream of the embankment (Ej) to the 
downstream of embankment (E^) is denoted as AE. This energy loss AE is non- 
dimensionaliscd witli upstream energy (E/) and is plotted against the variation in water level 
in the upstream and downstream of the embankment {AH ) as shown in Figure 4.39. This AH 
value is non-dimensionaliscd with the base length of the embankment Using the ratio of 
average shear velocity! E.) with sediment resisting velocity along with Qp^;^, the 

regression analysis is carried out. The polynomial relation between AE/E; and AH/Lp)s is 
written as: 


AE 

E, 


where , S 


O.On + 0 .S 7 S- 0 . 23 S" 


X Q 


(4.21) 


A//] 

0.41 

s/ 

rnM 

0.24 

V 

^ ^DB ^ 

A 




0.51 

PSN 


This functional relation for 6 was obtained from regression analysis of the data. Energy loss 
will be higher during the initial stage of erosion and it decreases as the breach erosion reaches 
final staue. 


g. Sediment I.oad Variation 

Volume of sediment eroded is computed from breach geometry at any given instant during 
breach. The bulk density of compacted soil in the embankment is known by measurement 
from the sample collected from the embankment. The product of volume rate of sediment 
discharge, Q^, with bulk density (/{,) is considered to be sediment load (gj = Qs Yb)- Sediment 
Ev is non-dimensionaliscd with ( d^oV^ji ) and is plotted against non-dimensionalised 




Fig. 4.39: Variation of nondimensional energy loss, AE with 
nondimensional AH 
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average shear velocity( V, ) as Ff / F^ along with Qps^ and as shown in Figure 4.40. 
Using the regression analysis, relation between them is written as: 

^,v “ ~ 2 . 12 + 4 . 3^/ - 0 . 17^* (4 22) 

'' p- ) 

yvhere,i//= ~ x x F^^f 

This equation is v^alici lor ^ 0.48. Wide scatter in the data indicates a poor agreement The 
relation shows that g, varies almost square of the shear stress. 

h. Variation of Sediment Concentration 


Sediment discharge {(7,v) is nondimcnsionalised with water discharge (Q,^), denoted as 
sediment concentration (O at a given instant of breach time. Using regression analysis 
between the variables /i,/7//, AH/L^n, Ff IFp^ and 2m a relation is obtained as 

indicated in Figure 4.41. The exponent of the variables are obtained from the regression 
analysis. Separate curves arc fitted for cohesive soil and cohesionless soil and also for 
combined data. Power law relation of the type given below in Equations (4.23 a), (4.23b), and 
(4.23c) respectively for cohesionless soil, cohesive soil and combined one are obtained. 

C. = g*'”* X 0.071 (Fme Cohesionless Sand - Type /Cj (4.23a) 

C. = c" X 0.073 {Cohesive Soil - Type B) (4.23b) 

C. = X 0.067 (Combined Data) (4.23c) 


where, <^ = 



- 0.95 

X 




0.15 


X 2 


0.16 

PSN 


/■ \ ().404 

1 iis r/ 0.222 

and C. = X F^.^^ 

\QwJ 

It may be observed that the power of the variable ^ and multiplying constant for both soils 
data and for combined data remains fairly constant. The scatter of data is very high. Hence 
these relation must be used with caution. 




Fig. 4.40; Variation of nondimensional sediment load gs*, with 
nondimensional bed shear stress 
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Fig. 4.41; Variation of nondinnensional sediment concentration 
with hu, AH, LnR, V, and Qpsn 
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4.6 Analysis of Laboratory Data Along with Field Data 

It was necessary to verify the trend and order of magnitude of some of the parameters 
evolved from the laboratory data, with the field data. These evolved parameters are: peak 
outfiow discharge ((?„,.), erosion of breach width depth of erosion (£>,,). These 

parameters are non-dinicnsionalised using flow and geometric parameters. Field data are 
collected from the book on ‘Dam Breach Modelling Technology’ by Vijay P. Singh. This 
book contains the data of 85 dam breaches all over the world. Data contains year of 
construction and failure along with the information on dam geometry and storage capacity of 
the reservoir. For these dams when breach occurs, the final dimensions of the breach geometry 
like depth of erosion, top width of erosion , and duration of erosion along with maximum 
outflow discharge arc recorded. These field data on dam breach contain part of the 
infomiation and only a few dam breaches contain all the required information. These data are 
used along with experimental data for analysis of breach parameters. 


a. Peak Outflow Breach Discharge (Owe) 

The maximum discliargc of outflow' is related to capacity of the reservoir and height of the 

dam as indicated in the I'igure 4.42. The relation between them is as follows: 

2 



0.00215 X 


^ Capacity^ ^ 


(4.24a) 


The magnitude of the inflow into the reservoir is not considered in this analysis because of its 
non-availability. Also the dam fill material properties are not known in most of the cases. The 
analysis of the model data indicates that the soil properties are also important variables in the 
prediction of Q^p.. 

b. Depth of Ero.sion {D^r) 

The final depth of erosion is plotted against in a dimensional form and with capacity of 
reservoir in a nondimcnsional form shown in Figures 4.43(a) and 4.43(b). The figure shows 
that height of dam (///)) and capacity of the reservoir are very important variables in the 
prediction of final depth of erosion. The relation between them is written as: 

= 0.737 or = 0.74 

and (MDj/Qwih = 33719.7 x {{H d /C apacityf^f^^^ 


(4.24b) 
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c. Top Width of Erosion (BJ^ 

Width of erosion is plotted against peak discharge (Q^p), keeping regime flow analysis in the 
rivers in view. Regime flow analysis indicates that width of the river is a function of the 
maximum dischaige passing thiough. The final breach width is plotted against peak discharge 
for both laboratory and field data. The relation between them indicates that is 

0 4 * ^ 

proportional to the • Lai and Shen (1996) showed from their experiments that the 

equilibrium width (5^,) of the channel is related to the dominant discharge as: 

5. = 10.91 (4.24c) 

Wliere 5,. is measured in meters and Q,j is in cubic meters per second. These experiments 
were conducted for the study on the flushing of sediment through reservoir. The magnitude of 
the exponent observed for the breach of a dam is very near to the exponent observed for the 
channel fomiation. In tlte present mixed data analysis of both field and laboratory, final width 
is non-dimensionalised with height of the dam and is plotted against dimensionless Qt^p in 
terms of ///) and capacity of the reservoir as shown in the Figure 4.44. The relation between 
them is as follow: 



28.4 


Qwp 


0.57 


/g Hi"' Capacity°-^ 


(4.24d) 


There is a wide scatter in the data due to the exclusion of inflow and soil properties. 




f^fcchanism of Dam Breach Process 1 58 


d. Relation Between Breach Fronde Niimhi>r rir % ^ 

TO «n<l Breach Shape Faclor (i).) 

The breach Fronde number is defined as: 


F = — ^ __ 
and breach Shape factor is defined as: 


(4.25a) 


These parameters are coined by Vijay P. Singh and P. D. Scarlatos 
Breach Froude number is related to .9^as: 


(4.25b) 

(1989). They showed the 


F, = 0.2 S/ ” 

This analysis is carried out 'using laboratory and field data 
shown in figure 4.45 and in a functional fonn may be written 


(4.25c) 

and the relation among them is 
using regression analysis as: 


/•; = 0.048 


(4.25d) 



Fia. 4.45: Relation between breach F/oude number (Ft) 
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4.7 Discussion of the Results 

In this chapter efforts has been made to know the details about the mechanism involved in the 
dam breach process, with both cohesive and cohesionless soil models by varying both flow 
parameteis and cmbankinent geometry. Required variables are measured with respect to time 
at suitable location along the dam body and in the reservoir. The dynamic nature of the 
breach problem is analysed using these data. 

The mechanism involved in this problem can be broadly categorised into three domain 
namely, 

a. Breach Morphology. 

b. Flow Hydraulics. 

c. Breach Fmosion. 

a. Breacli Moi jrhology 

i. From the experiments, it is observed that the breach cross-sectional area does not remain 
same along the breach at any lime of breach, except towards the end of the breach. This is 
shown in Figure, v 4.1(a) and (b) for cohesionless soil and in Figures 4.25(a), (b), (c), and 
(d) and 4.26 for cohesive soil embankments. This needs to be taken appropriately in using 
the routing procedure for flow through breach channel. This variation of breach geometry 
implies that there will be no uniform erosion along the length. 

ii. The cross-sectional area of the breach changes from trapezoidal at the upstream edge (or 
near crest) of the dam, to rectangular (or steep side slope) at the downstream edge, where 
discontinuity of surface slope occurs. Again the shape modifies to trapezoidal shape upto 
downstream end of the dam. The computation of sediment load eroded at any time is done 
through taking the average of cross-sectional area between any two measuring sections. 

iii. At the initial phase of breach there is cascading of flow particularly in cohesive soil. 

iv. The limitation or the shortcomings of the present study mostly goes to lack of automation 
in recording of the bed profiles both in longitudinal and lateral direction and water surface 
profile along breach with time. During the initial period of breach the water will be turbid, 
particularly near the sections 6 and 7 (as shown in photographs attached in Chapter III), 
near the downstream top edge and at the section between downstream edge and toe of the 
dam. So the bed level measurement with time by the help of point gauge or even with 
suitable modification to point gauge as described in Chapter III may not be error free. To 



f^schanism of Dam Breach Process 160 


overcome this limitation, actual values of water level, bed level and water surface widths 
are plotted with time. An average curve passing the points is fitted. The readings 
conesponding to this average curve are used as data points in the analysis. 

V. Collapsing of bread) banks are not analysed explicitly. Corresponding width enlargement 

due to side collapse is considered in computation of sediment volume. Conditions under 

which banks are collapsing is not considered in the present analysis. 

b. Breach Hydraulics 

i. From the measured and Q^,, it is confinned that the broad-crested weir formula can 
suitably be used in the dam breach hydraulics. The magnitude of discharge coefficient 
varies with the hciglit of sediment crest level and head acting over it. 

ii. The effect of tail-water is not considered in the analysis. 

As observed from the experiments and also from the historical dam failure cases, dam body 
never gets eroded \'crtically upto its bottom. It can be seen from Equation (4.24b), Der = 
0.74 //;, for both lab and field data. So the chances of submergence of flow for very low 
value is rare and even if it happens, that will be towards the end of the breach. By 
that time the flow' in the breach channel will be of channel type flow without much further 
erosion cither in vertical or lateral direction. So one may fairly assume free weir flow 
situation without undue loss of accuracy. 

iii. An interesting observation about the breach crest is that it moves in upstream direction, 
particularly the flow' crest moves very fast near the zone of downstream top edge of the 
dam, where there is discontinuity of slope surface. Though the crest moves in upstream 
direction, it is observed during most of the time of the breach process, the crest lies near 
the upstream top edge of the dam. Though it depends on the top width of the dam section, 
still, it will be appropriate to locate the crest near the upstream top edge of the dam without 


substantial en'or. 
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C. Breach Erosion 

i. Erosion in the breach section occurs both in vertical and lateral directions. Vertical erosion 
is mainly due to erosion caused by flow velocity. Lateral erosion or widening of breach is 
due to banK collapsing. Oank collapsing occurs due to caving of the banks at the comer of 
the base. Caving causes hanging of banks, with water flowing below resulting in their 
collapse. Permeability of soil fill plays a dominating part in the erosion process. Soil fill 
having less pemieability is more resistant to vertical erosion. In the present analysis, 
explicitly the effect of permeability of the soil fill is not considered. 

ii. The dominating factor in breach erosion is the acceleration of flow. This is controlled 
particularly by the magnitude of water level drop between the upstream and downstream of 
the dam. 

iii. Cascading and water fall types of flow occur along the breach in cohesive soil 
embankment. These are absent in cohesionless soil embankment breach. Cause for these 
types of flow is mainly due to pre.sencc of cohesion. 

iv. The flow during initial phase of erosion appears like dense muddy flow in cohesionless soil 
embankment. 

4.8 Concluding Remarks 

a. In the breach hydraulics, reservoir water level (A,,,,.), head over crest {hj, water level 
difference {AH) are suitably non-dimensionalised and analysed for their variations. 

b. The corresponding scales for the above variables are obtained for the use at any other given 
situations. The scales are related with known parameters of dam geometry and properties 
of its fill materials; reservoir size and breach peak outflow discharge. 

c. The related time component associated in these variables are also analysed. In the breach 
morphology, width and depth of erosion are studied extensively along with their rate of 
variation. 

d. The photographs arc taken at different stages of the cohesionless embankment breach 
which indicates the process involved in breaching of the embankment. 

e. Though photographs give qualitative information about the changing nature of morphology 
with time during the breach process, the rate of variation of depth of erosion of crest(/3<,„) 
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of sediment bed and width of erosion at crest (S„) reveal quantitative infomation about 
their mutual relation. 

f. Erosion of sediment and its transport occurs mostly with accelerated flow conditions in the 
dam breach process. Sediment transport relations available in literature are developed for 
the uniform flow condition. Due to this basic change in flow condition, transport formula 
used in most of the numerical modelling of dam breach is not used in present analysis. 

g. The average shear velocity ( F. ) in breach channel along with water drop {AH) and the soil 

resisting velocity {V^f) are used in developing functional relationship between sediment 
transport and flow properties. 

h. There is a strong relationship between the peak water discharge and peak sediment 
discharge and their corresponding times of occurrence during the breach process. 

i. The maximum depth of erosion and maximum width of erosion at crest are interrelated to 
each other. 

Now with this background of dam breach analysis, the functional relations developed can be 

used to validate the field data. 



CHAPTER V 


COMPARISON OF DAM BREACH RESULTS 
WITH TETON AND HUACCOTO DAMS 


5.1 Introduction 

The developed functional relations among the system parameters from the analysis of 
experimental data need to be tested with prototype data for their validations. As any numerical 
model heavily relies on the experimental observation to test the validity of the methodology 
adopted in computational algorithm, the experimental results similarly needs a well accepted 
numerical model for its comparison . Sometimes the system is so complex and dynamic in its 
behaviour, like the present problem of dam breaching, it is not easy to get the process 
information easily. This compels the investigators to bank upon the systematic, well planned 
experimentation on scaled models. In open literature very few experimental results are 
documented in this topic. Those which are available do not deal so much on the mechanism of 
dam breach process behaviour to bring much light on the subject. Immense effort has been 
made in the development of numerical solution to dam breach problem. Very few approaches 
have attempted with a proper perspective insight into the basic mechanism. Due to this, many 
numerical approaches become invalid for a given situation. A few approaches, to some extent, 
like BEED model are successful to a fair degree of acceptance. The knowledge on mechanism 
of breach process from the experimental studies will be helpful immensely in further 
improving the numerical approach to solve real life problem on dam breaching. However, it is 
necessary to validate how far the functional relationships developed from the model 
experiments will predict prototype behaviour for their authentication and further use . Almost 
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none of prototype data is available during the breach process. This gap may be filled by both 
qualitative and quantitative observations of the breach process by the model study in a 
laboratory set-up. The final verdict on the laboratory data is only how close they predict final 
breach geometry of the prototype dam failure. For the validation of the present study, two 
prototype dam fiulures information are used as reported by Singh (1996) namely; 

a. Teton Dam 

b. Huaccoto Dam. 

Ponce and Tsivoglou (1980), Singh and Scarlatos (1989) have verified their result with these 
two observed dam failure data . The present laboratory analysis result is compared with the 
result of BEED (Breach Erosion of Earth-Fill Dams), a numerical model for the breach 
associated with Teton dam and Huaccoto dam. Before going into the validation of the result, a 
brief note on Teton dam, Huaccoto dam and BEED model will enable one to know their 
background. 

a. Teton Dam 

Teton dam, Idaho, U.S.A. an carthfill structure of height 93m, experienced failure on 5th 
June’ 1976. The details of the dam geometry and the properties of its fill material are given in 
Table 5.1 . The dam had a crest width of 10.5m. Clay, silt, sand and rock fragments were used 
as dam fill. The full reservoir capacity was 3.1x10'^ m’\ The dam cross-section is shoA^m in 
Figure 5.1(a). 



Elevation (m) Elevation 
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Max Water Surface 

1623.0m 1,625.2m 



Fig. 5.1(a): Cross-sectional details of Teton dam 



Fiq 5 1(b)- Cross-sectional view of Huaccoto natural dam 
built by landslide material. 
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table 5.1 


CHARACTERISTICS OF DAM GEOMETRY AND SOU Ftl i 
PROPERTIES OF TETON DAM AND HUACCOTO D^ 


.No 

Description 


Teton Dam 

Huaccoto Dam 

1 

R e s erv oi r Capa city 

“ ~(mT 


8.87x10* 

2 

Dam Height 

(m) 

93 

= 170 

3 

Lateral Length of Embankment L| 

DB (m) 

550 

3800 

4 

Crest Elevation above mean Sea level 


2630 

5 

Top wid th of Dam 

(^dt) 

29.6 


6 

Aq (cross-sectional area) 

(m^) 

24,986.71 

= 43,6500 

7 

d^o (median size of soil) 

(mm) 

3 

11 

8 

(j) (fiictional angle) 

(degree) 

40 

40 

9 

('(cohesive strength) 

(K.PA) 

49 

40 

10 

Yj^(unit weight of dam fill) 

kgW 

2300 

2300 

1 1 

12 

average inflow,,. 

'Fype of Soil 

(cumec) 

106.2 

Silty sand and some clay 


b. Huaccoto Dam 

The Huaccoto natural dam was created due to landslide that occurred in Cochacay Creek, a 
tributary of the Montaro River in Peru on April 25, 1974. The landslide had a volume of 15.8 
X lO'' m'^ approximately, made this artificial dam and created a huge reservoir of depth nearly 
168 meters which took nearly sixty days to fill. It can be quoted from the report. 


“ .... Indeed, the volume of the material was so high that it would be enough to build an engineered 

dam about three times at the same location ...” 

The overall geometric characteristics of Huaccoto dam are shown in Table 5.1 and the dam 
cross-section in Figure 5.1(b). The dam was overtopped during June 6-8, 1974 (Lee and 
Duncan, 1975). The length of the downstream slope length was around 1.069km. The 
overtopping flow during the first two days eroded gradually, a small channel along its 
downstream slope surface. Then after 6-10 hours, there was an increase in the breach channel 
resulting in a trapezoidal .shape approximately 105m in depth and 195-225m in top width with 
side slopes of about 1 ; 1 . 'The breach did not extend up to the bottom of the embankment. As 
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a result, a large lake remained after the breach had subsided. The estimated breach related 
characteristics as mentioned by Vijay P. Singh were : 

Maximum flow depth = 1 5 ~ 20 m 

Time for discharge 0 m'Vsec to 100 mVsec = 16hr 

Time for discharge lOOm’Vscc to peak discharge = lOhr 

hydrograph duration from lOOmVsec rising to 400 m^/sec receding = 32hr 

C. DEED Model 

There are various numerical ftiodels on dam breaching problem as listed in chapter - 2. SEED 
(Breach Erosion of Earth Fill Dams) model is the recent one (1989), encompassing most of 
the features of dam breaching. Since then no further numerical studies on this subject have 
been published and tested on any historical event to the knowledge of author. Seemanpali and 
Singh (1992) have approached this problem using variational calculus but still DEED model 
seems to be more appealing and elegant and widely accepted for application. Hydraulic 
Research Wallingford Limited, Wallingford, U.K. (1995) have carried out large scale model 
study in sand benns, however, the numerical model based on these test data are still in the 
developing stage as reported. So keeping in view this status of knowledge in numerical 
modelling on this topic, BEED model has been chosen for the comparison of the present 
experimental results, 

i. Background of BEED Model 

This model is developed by Singh and Scarlatos (1986, 1987a), Louisiana State University, 
Baton Rouge for Department of the Army of Engineers, Washington. This model has been 
tested with four historical earth dam failures namely, 

1. South Fork Dam, Pennsylvania (1989) 

2. Buffalo Creek Dam, West Virginia (1992) 

3 . Teton Dam, Idaho (1976) 

4 . Huaccoto Natural Dam, Peru (1976) 

Some salient features in this numerical model and the methodology involved in its evolution 
are presented here. At this stage, it is worth to quote the statement of Singh and Scarlatos 
(1989) about this numerical modeling in the report to understand the inherent limitation 
associated with this problem : 
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“ Simulation of the total earth fill dam-breach erosion process is a combination of hydrologic 
events, hydrodynamics, sediment transport mechanics and geotechnical aspects. The real life 
problem is unsteady, non-homogeneous, non-linear, and three dimensional, which is not 
theoretically well understood. Mathematical modelling of the phenomenon requires idealisation of 
the real-1 ifo situation so that the leading physical processes can be described by a set of governing 

For practicality, there is always a trade-off involving complexity, accuracy, and 

efficiency of the model “ 


Keeping the above statement in mind and the necessary idealisation of the system, one should 
look into the five components of the model namely, 

i. reservoir water balance 

ii. breach geometry 

iii. breach hydraulics 

iv. flood routing 

V. sediment routing 

Among the above stated live components, present study is only concerned first three features 
of the SEED model. In reservoir water balance, the discharge through spillway and any other 
structural outlet like sluice gate, during the breach process if any can be accommodated in the 
model. In this model, the longitudinal flow direction is demarcated into two distinct sections 
namely breach section and breach channel. The flow over breach section is similar to flow 
over a broad crested weir. It assumes the validity of Manning’s Equation in breach channel. In 
the breach enlargement, the Einstein-Brown formula is used for the transport of sediment in 
the breach. The sudden breach enlargement is taken care of by slope stability of breach 
channel by contour method. But in the BEED model it is difficult to predict a priori the time 
at which stability will overtake breaching. After instability occurs erosion again predominates 
until another slope stability occurs. 

d. Breach Variables to be Compared 

In the validation of the present experimental results, following parameters are compared with 
observed data of Teton Dam and Huaccoto Dam and simulated results of BEED model. There 
is no observational information on sediment related data like peak sediment discharge 
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intensity and time of occuirence of peak sediment discharge Hence, results on q, 
and t,p are compared only with the simulated result of BEED model. 

Major influencing parameters of breach process to be validated are: 

i. Peak outflow water discharge (Q^yp). 

ii. Time of occurrence of peak outflow water discharge 

iii. Peak sediment discharge intensity (q^p). 

iv. Time of occurrence of peak sediment discharge intensity (t^p). 

V. Final maximum top breach width 

vi. Final depth of erosion {D„). 

vii. Final side slope of the breach cross section (1 V:zH). 

In addition to the prediction of these parameters, the total outflow hydrograph is compared 
with observed and simulated results of BEED model for both Teton Dam and Huaccoto Dam 
in nondimensional form, Qn/Qfyp versus t/t^p. As the observed information of sediment 
discharge with respect to time is not available, it is compared with the simulated result of 
BEED model for the Teton Dam only. 

Before discussing about the validated result of above parameters, it is necessary to 
introduce again at this stage the various nondimensional numbers used in the present analysis 
(which are introduced in Chapter-IV) namely: 

1. Q-peak number (Qpi^) 

2. Q-peak-sediment number (Qpsn) 

3 . Q-peak-cross-sectional number (Q^pn) 

4 . Sediment resisting velocity number 

5. Capacity-inflow number {CAPj^) 

6. Shield’s entrainment number ( t*c) 
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TABLE 5.2a 


VARIOUS NON-DIMENSIONAL NUMBERS FOR THE EXPERRTFvtc 
CONDUCTED ON DAM MODELS (Darn Models with Cohesionkss Sod) 


S, No 

Run Nos. 

Qpn 

Qapn 

Qpsn 

VSR.N 

CAPjj 

(t-c) 

1 

10(}.5 

0.02820 

0.00416 

0.296 

0.0140 

0.00239 

0.0729 

2 

0805 

0.03170 

0.00423 

0.279 

0.0152 

0.00769 

0.0729 

3 

0705 

0.03167 

0.00366 

0.220 

0.0166 

0.00880 

0.0729 

4 

0505 

0.08170 

0.00802 

0.427 

0.0188 

0.00324 

0.0729 

5 

0403 

0.06980 

0.00560 

0.256 

0.0218 

0.00166 

0.0729 

6 

03 05 

0.17040 

r 0.01320 

0.616 

0.0214 

0.00903 

0.0729 

7 

0105 

! 0.13160 

0.00990 

0.455 

0.0218 

0.00921 

0.0729 

8 

3004 

0.06580 

0.00500 

r 0.230 

0.0218 

0.00970 

^ 0.0729 

9 

2904 

0.16950 

0.01378 

0.632 

0.0217 

0.00904 

0.0729 

10 

2804 

0.36700 

0.05530 

2.540 

0.0217 

0.00590 

0.0729 

1 1 

2604 

0.34300 

0.09700 

4.450 

0.0217 

0.00190 

0.0729 

12 

2504 

0.37100 

0.06620 

3.040 

0^0217 

0.01060 

0.0729 


TABLE 5.2b 


VARIOUS NON-DIMENSIONAL NUMBERS FOR THE EXPERIMENTS 
CONDUC TED ON DAM MODELS (Dam Models with Cohesive Soil) 


S. No 

Run No. 

Qpn 

Qapn 

Qpsn 

^SRN 

CAP^, 

(x-c) 

I 

1203 

0.1060 

0.028 

0.640 

0.0440 

0.00295 

0.1612 

2 

2703 

0.0550 

0.016 

0.436 

0.0372 

0.00350 

0.1612 

3 

()4()4 

0.0449 

0.014 

0.404 

0.0342 

0.00330 

0.1612 

4 

1 104 

0.0441 

0.032 

0.446 

0.0315 

0.00410 

0.1612 

5 

2204 

0.1000 

0.111 

1.037 

0.0310 

0.01880 

0.1612 


(T ) + r 

V Oc / cohesion \ oc J friction 

(r ) = 

g^50 

The calculated numerical values of these numbers for Teton Dam and Huaccoto dam are given 
in table 5.2(c). 
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table 5.2c 

VARIOUS NONDIMENSIONAL 
NUMBERS FOR TETON .AND HU.A.CCOTO D.AMS 


Number 

Teton Dam 

Huaccoto Dam 

Qpn 

0.2488 

0.1161 

Qpsn 

9.1325 

0.13259 

Qapn 

0.086 

0.000768 

"'^.SRN 

0.00943 

0.005796 

CAPn 

11.185 

- 

(T-c) 

0.162 

0.0703 

(cni/sec) 

28.48 

23.679 


Q.. 

10421.59 

- 

iohPjr 

a 

6549.87 

- 

'•<« /‘.'V 'ir 

a 

7550.39 

- 


where 




Qpxn 


Qapn 


V 

' SHN 


g., 

■IWd 

Qwp 

SR 

Qwp 

V S^D 


CAP^ 


Capacity * 


T.. 


[x ) 

V ^<^7 combined 

^Pgd,o 


V = 

^ SR 


l~g^5o tan^> + — 

Pw P yv 



Comparison of Dam Breach Results 


172 


5.2 Comparison of Estimated Breach Variables 


a. Duration of Breach 

Duration of breach is taken as the time taken for the 95% of the soil to be eroded by the 
overtopping flood and denoted as The multiple regression analysis is carried out with the 
present expenmental data and field data of Teton Dam for different possible parameters and 
finally it was found that the is the dominant parameter which correlate with magnitude 
of the breach duration well. Since there is an inconsistency in the duration of Huaccoto dam 
breach, as reported in literature, this data is not used in the analysis. From the regression result 
of laboiatoiy data and I eton Dam data it is found that the non-dimensional duration of breach 

time (/«) taken as varies proportional to Then the ftmctional relation 

among them is found out to be with regression coefficient, R = 0.896 as: 


4328650 (5.1) 

Figure. 5.2 shows a good correlation between experimental data and observed data of Teton 
dam in the torm of Equation (5.1). Table 5.3 shows comparison of estimated breach duration 
with observed and BEED simulated result for Teton dam failure. 


TABLE 53 

COMPARISON OF ESTIMATED DURATION OF BREACH (t,;) 
WITH OBSERVED TETON DAM DATA 


Output 

Regression 


Teton Dam 


Variable 

Scenario 






Estimated 

Observed BEED 

%error 

Breach 

Lab Data 

1.46 

1.5 1.35 

2.7 

duration 

alone 

hs 

Lab+Teton 

1.5 

1.5 1.35 

0 

(hours) 

dam data 


The percentage of error listed in Table 5,3 for the estimated value of the variable is with 
respect to its observed value. 
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b. Peak Breach Outflow Discharge (Qfyp) 

Peak outflow discharge of a breach {Q^yp) is one of the important parameters used in the 
analysis of breach process and in estimation of breach dimensions. Regression analysis is 
carried out to find major influencing parameters in developing relation for Q^yp.. Major 
influencing parameters arc Ap, 0„„ Hp, Capacity of the reservoir, Vsp and r.c. 

T*c [ ( g 

where is the tractive shear resistance offered by the sediment bed. The (roc),ohesioniess can 
be found from Shield’s diagram knowing d^o, and But there is no direct approach to get 
the tractive shear stress for cohesive soil. So with this limitation, to get (roc)cohesive> it is taken 
to be 15% higher than shear strength of soil as reported by Raudkivi (1990) in the book 
‘Loose Boundary Hydraulics’, (page-318). Hence (r„.)eohesive =1-15 c. The calculation of r.c is 
shown in appendix , where r is cohesive stress of the soil. Here again, in the case ofHuaccoto 
Dam, inflow information is not available. From multiple regression of experimental result it is 
found that depends on along with height of the dam, capacity of the reservoir, cross- 
sectional area of tlic dam. Shield’s parameter and soil resisting velocity (Kj^j). Since average 
inflow into the reservoir of Huaccoto dam is not available, its data is not used in regression 
analysis. Data of I'eton dam and the laboratory experimental data are taken for regression 
analysis. The functional relation obtained as per regression analysis is indicated in Figure 5.3 
and given in Equation (5.2) as: 


iWP 




2.78 


.Cap VsJ 


0.22 


Xk 


(5.2) 


The values of « and A- along with estimated values of Qwp, and their corresponding observed 
value and simulated value are also listed in Table 5.4. Figure 5.3 contains only the second 
case as indicated in Table 5.4. 


TABLE 5.4 

COMPARISON OF THE PRESENT ESTIMATED VALUE OF Qwp FOR TETON DAM 
BREACH Wnn OBSERVED AND COMPUTED V ALUE BY BEEP MODEL 

I^JresIi^Sa-n'arTo Values of n. k in E qn 5.2 ~ (m /sec) For Teton Dam_bre^ 

n k Estimated Observed BEED % Error 

Present i;xp. data 0.89 5.07 67,028 65,000 68,581 3d2 

.Present t Obsv. 0.89 4.93 66,390 65,000 68,581 2.1 

Telon Darn Data — — 



Qwp/ 


Fig. 5.2: Varia 
(tgs) 


(Q*pn) = CAP,“-^^)°-«%4.93 


00000 Present Experimental Data 
imiiiia Teton Oom Dota 


’4 5 8 7 8 9 


CAPn°-"^ 


Fiq 5 3- Variotion of nondimensional breach peak 
■ discharge (Qwe) with Shield's entrainment 
■number (t.c) pnd reservoir capacity 


numb 
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From the above obsen/ation of the value of.’ and T, it may be reeonmended for further use 
of the Equation (5.2) with n=0.9 and k=5.0. 

c. Time of Occurrence of Peak Outflow Discharge (t„p) 

It is essential to know the time required for the maximum breach outflow to occur in a given 
condition to monitor the flood downstream. Regression of lab data. Teton dam and Huaccoto 
dam data indicates that is a function of , Shield’s entrainment number (r^) and 

peak outflow number (<2;,^.). The relation is given in functional form as indicated in Equation 
(5.3). 


wp 


*95 


0.035 0) 


( 5 . 3 ) 


Wi 


vherc <!>, = y ^ j 

The results are also shown in Figure 5.4. The relation developed for t^.p in Equation (5.3) 
includes the Huaccoto dam data. A comment on the inclusion of Huaccoto dam data in the 
analysis is necessary. Since the ratio is used for regression analysis, any errors in 
and ty, may cancel out to certain extent. This can be observed from Figure 5.4. Error in 
estimation of is shown in Table 5.5. Large error in the estimation of Huaccoto dam is an 
indication of uncertainty of the observation of 


TABLE 5.5 

COMPARISON OF ESTIMATED TIME OF OCCURRENCE OF PEAK OUTFLOW 
DISCHARGE (l„p) FOR TEl ON DAM AND HUACCOTO DAM WITH OBSERVED 
AND SIMULATED VALUE BY BEED MODEL 




Teton Dam 

Huaccoto Dam 

Output Vanalilc 

KcgirruMon S( rriuito 

F.siinuitcd 

Observed BEBD ‘hicrror 

Estimated 

Observed 

BEED 

“terror 

Time to Peak 

Tab i I cton f 







outflow Discharge 

1 luaccoto Dam 

1.27 

1.5 1.35 15.3 

22.8’ 

=17.5 

N.A 

30.3 

(Lrt)hour 









Observed value of 45 is used in the calculation. 
N.A = Not available 


d. Width of Breach (B^,)/ 

Final breach width is one of the major breach morphological features. It represents erosive 
power of the breach flow discharge in embankment fill. It also gives a qualitative picture of 
the breach process. From the regression study of Teton dam data and present experimental 
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data the functional relationship among , Q„ and is shown in Equation (5.4). The 
relation among them is shown in Figure 5.5. 

H J f -r X OPN ] (5.4) 


Here x.c contain both the component of tractive shear stress, due to the effect of cohesion (c) 
and friction(^) of the soil. The functional relation for different situations are presented below. 


K)j_ 

Hn 

//n 


= 1.2 ln(r°,= X f )+ 3.9 (Only lab data) 

0.99 ln(Jr"^^ x 3.56 {Lab and Teton dam data) 

0.74 ln(r 1',. x Q'l^ 3.0 ( Lab, Teton dam and Huac cot o dam data) 

0.87 ln(^rl’,.^ x 3.5 (Re commended) 


(5.5) 


Final estimated breach width for Teton and Huaccoto dam is compared with observed value 
and that of BEED model simulation value in Table 5,6. 


TABLE 5.6 

COMPARISON OF ESTIMATED FINAL TOP BREACH WIDTH 
FOR TEl ON DAM AND HUACCOTO DAM WITH OBSERVED AND 
SIMUIaATED VALUE BY BEED MODEL 




Teton Dam 

Huaccoto Dam 

Output \'ai iablc 

Hrgtcssuni Sccnatiti 

Estimated 

Observed 

BEED %crror 

Estimated 

Obscrv'cd 

BEED %aTor 

Final top Breach 

Only l.ab Data 

206.0 

200 

295 +3 

72 

200-230 

298.5 

66 

Width (B,,), (m) 

eton Data 

202.0 

200 

295 +1 

77 

200-230 

298.5 

64 


Lab""f'"reton ■+ Huaccoto 

182,9 

200 

295 -2.2 

146.8 

200-230 

298.5 

38.6 


Reeom,meiided 

212.5 

200 

295 +6.3 

168 

200-230 

298.5 

22 


From the comparison of the estimated result with observed, the one which gave least 
percentage enror is recommended for use. The recommended equation for the computation of 
final breach width is: 




((Bcr)f/Ho) 
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Y*«X’‘**0.035 
Teton Oom Data 
AA A A A Huoccoto Oom Data 
ooooo Present Data 


; n 


((Vsrh'^O Qph) 

5.4: Variation of nondimensionol time of 

occurrence of breach peak discharge(t*p) 
with Shield's entrainment number, t,c 
and QpN, soil resisting velocity number, 


Y^=0.751n(X)^,0 ,3 

Where X= t*c (Qwp/( 9 ^o j j 


ooooo Present data 
nuuia Teton Dam Dota 
jlaaax Huoccoto Dom Data 


(T.r*(Qwp/(gHo*)“-^)“* 

5.5; Variation of nondimensional final top 
breach crest width, (Bcr)f with Shield s 
entrainment number, t*c Qpn 







Comparison of Dam Breach Results 178 




= 0.87 ln(Z) + 3i 


(5.6) 


where, X'- 




PN 


In the recommended relation slightly higher value of for Teton dam is accepted to get a 
nearer value to Huaccoto Dam data. 


e. Depth of Breach (D^r) 

From the observation of the results obtained from regression analysis of the experimental 
data, it is found that the the maximum lateral erosion is closely related with maximum 
vertical erosion (D,.,.) and the relation among them is shown in Figure 5.6 and the functional 
relation between them may be written as; 

= [0,w X ] * 4.47 (5.7) 


TABLE 5.7 

COMPARISON OF FINAL BREACH DEPTH (D„) 
FOR TETON AND HUACCOTO DAM BREACH 




Teton Dam 

Huaccoto Dam 

Output Variable 

Rcgr'cs.su.;iMfi Scenario 

Estimated 

Observed 

BEED 

Serror 

Estimated 

Observed 

BEED 

%crTor 

Final Depth of 

Lab+'Teton 

90.6 

79 

93 

+14.7 

97.8 

107 

85 

-8.6 


Eqn 5.7 









Breach (D„) (m) 

Recommended 

84.6 

79 

93 

+7.1 

91.1 

107 

85 

-14.8 


Eqn 5.8 










Equation (5.7) is slightly modified as Equation (5.8) and it is recommended for use. 

%^ = [e«,xcrx''.42 (5.8) 

Recommendation is based on minimum error indicated in Teton dam breach, 
f. Side Slope of the Final Breach Cross-section (z) 

During the experiments it is observed that the breach flow cross section remains trapezoidal at 
the crest for most of the breach duration. The side slope of the final breach cross-section is 
expressed as one vertical: z horizontal (TV;zH). Crest being the controlling section of flow 
through the breach, it is essential to predict the final possible side slope of the breach cross 
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section for a given condition. From the regression analysis of the present experimental data 
some field data including Teton dam data, it is found that the side slope of the cross section is 
a function of maximum lateral and vertical breach erosion dimensions as well as height of the 
dam and peak breach discharge (,Qiyp)- Their relation is shown in Figure 5.7 and the relation 
between them be written as: 


where, ij/ ^ = 


z = X 0.107 

D.. 


n 1.49 


HT (b„)T 


X Q 


-0.06 

PN 


(5.9) 


Using Equation (5.9) and the calculated value of and (5 from recommended relations, 
the final side slope of breach section is compared with the observed and the result of the 
BEED model in Table 5.8 as shown below. 


TABLE 5.8 

COMPARISON OF ESTIMATED SIDE SLOPE OF BREACH CROSS- 
SKCrnON FOR TETON DAM AND HUACCOTO DAM FAILURE 



Teton Dam 

Huaccoto Dam 

Output Variable 

Rcgrcssiion Sccnano 

Estimated 

Observed 

BEED 

% error 

Estimated Observed BEED % error 

Final Slope of 







Breach Cross 

Equation 5.9 

0.57 

0.5 

0.675 

10 

= 0.52 1.0 1.67 48 

Section (Z) 








It may be noted here that the percentage error for Teton dam is on higher side. Estimation of 
side slope for Huaccoto dam is relatively very high. 

g. Peak Sediment Discharge Intensity ( 

Sediment erosion in a breach depends on the magnitude of water discharge and soil properties 
of embankment. Eroded sediment volume distribution in a breach is studied by 
nondimensionalising it by its peak value. The magnitude of sediment peak discharge 
intensity {q^^) is computed from the credible sediment volume and its bulk density. 
Magnitude of {q^p) is never measured in the field during breach. It is always computed using 
some numerical models. This value for Teton dam breach is computed from BEED model. 








Fin. 5.7: Variation of side slope (z) 

section with final depth and width of erosion 

at crest 
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The BEED model has numerically simulated the sediment discharge and has found for Teton 
Dam failure that Qsp = 70.53 10^ kg/sec. This Q^pis based on predicted final top breach 
width of 295m. Based on this peak sediment discharge intensity by BEED model (q ) is 
worked out to be 103.95mVsec/m, with average specific weight of soil used for that dam, is 
2300 kg/m^ ■ The peak sediment discharge intensity is found to be the function of height 
of the dam(//^) and soil resisting velocity number (Fsjjy) as shown in Figure 5.8(a). The 
functional relation based on regression analysis of both lab and BEED model’s simulated data 
carried out separately and combinedly, is given as: 

X 10"^ {for only lab data) (5.10a) 

/ / O ' 

^ 3-6S7 X 10"^ {for lab and Teton Dam data) (5.10b) 

\dwp 

In Table 5.9 the sediment peak discharge {q^fj of Teton Dam by the present method of 
estimation is compared with that of simulated BEED result. Percentage of error is quite high. 
The relation is modified as indicated in Equation (5.10c) to reduce the magnitude of error. The 
recommended Equation (5.10c) is as follows: 


= 5,== X 3.67x10-“ 

^WP 


(5.10c) 


TABLE 5.9 

COMPARISON OF ESTIMATED AND SIMULATED BEED 
RESULT AND SEDIMENT PEAK DISCHARGE FOR TETON DAM 



Teton Dam 

Huaccoto Dam 

Out put Variables 

Regression Scenario 

fist j mated 10* 

Simulated BEED Result 

% error 

Estimated 10* Simulated BEED Result % error 

Qsp 

Only lab Data 

79.6 

103.95 

23 

- 

(cmVsec/cm) 

Lab+Teton Data 

74.6 

103.95 

28 

- 


Recommended 

93.8 

103.95 

9.78 

31.9 


♦ for the calculation of q,p , Bcr = 295m is used, simulated result of BEED model 


Analysis is not carried out for Huaccoto dam breach because of nonavialibility of information 
on simulated value of q^p 

The simulated sediment discharge Qsp for Teton dam by BEED model is obtained based on 
computed final top breach crest width 295m. Hence for the validation of present functional 
relation for q,p, the same {8,^)/ value of 295m is chosen instead of the value estimated firom 
the present method. 
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h. Time of occurrence of Sediment Peak Discharge (t^^) 

Like V is also another scaled parameter which is to be estimated for sediment discharge 
prediction with time. From the regression analysis of present data it is noticed that f is 
closely related to Q^p, Ap, and V^p . Two different regression scenarios one with only 
laboratory data and the other with laboratory data and Teton Dam data are undertaken. The 
functional relations obtained from these two scenarios are as follow: 

- — 0.27 ) + 0.68 {only lab data ) (5.11a) 

i_ = 0.1 8 ln((2,.x;^ ) + 0.56 ( lab and Teton Dam data ) (5. 1 lb) 

* wp 

In Figure 5.8(b), lab data and Teton dam data are fitted with the Equation (5.11b). A 
comparative study is presented in the Table 5.10. The estimated value oft in both cases 
predicts less than DEED simulated value. Equation (5.11b) is modified slightly as given 
below and recommended for use. 

j^ = 02\niQpsN) + 0f> (5.11c) 

TABLE 5.10 

COMPARISON OF ES I IMATED TIME OF OCCURRENCE OF SEDIMENT 
PEAK DISCTI AROE FOR TETON DAM AND HUACCOTO DAM BREACH 

t^p in hours 



Teton dam 

Huaccoto dam 


Functional Relation from 


Present 

Present 


Regression Analysis 

BEED 

Estimated 

Estimated 

Regression Scenario 

1 

DC 

1.3.5 

1.22 


Lab and Field Data 

^ wp 





‘sp 

027 In (C|>.W ' 

1,35 

1.56 



^ wp 





Recommended 





^sp 

- 0.2 ]iMQp^^) A 0,6 

1.35 

1.32 

4.45 
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For Huaccoto dam breach the simulated value of 1 ,^ from BEED model is not available for 
comparison with present estimated result. From the laboratory embankment breach analysis it 
is observed that /jp is always less than t^,p. It means that sediment peak discharge reaches its 
peak value before the peak water discharge. This observation is in agreement with the findings 
of Graf and Song (1995) in the study of bed shear stress in non-uniform and unsteady open 
channel flows. 



velocity number, Vsrn 






5.3 
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Muc„ce of Average Inflow Discharge! g.), Capacity of Reservoir 
(Cap), and Cross-Sectional area of Embankment f A t on EilT 
Breach Morphologicai Parameters ^ 


In the comparison of some of the important parameters in dam breach process in last secHon 
it is found that the developed functional relations for Teton dam are in good agreement with 
these parameters. To improve these relations further the effect of capacity of the 
reservoiffcap) and the average inflow( g, ) into the reservoir are included in the analysis of 
final morphological breach parameters like (B„), D,„ A.„ and Q„,. These final breach 
parameters are suitably nondimcnsionlised with An and g, and "plotted against as 

shown in Figures 5.9(a), (b), (c), and (d) respectively. To validate these functional relations 
again Teton dam is taken into consideration. In case of Huaccoto Dam Qm s 0, as reported 
(Singh, 1 995), hence it is not included in the analysis. The Q, for Teton dam is taken as half 
of the peak inflow discharge to the reservoir instead of average Qin. This assumption is used 
because inflow hydrograph was not available for Teton dam. But one should take average Qm 



Volume of inf low water 
Total base time of hydrograph 


In Table 5. 1 1 , a comparison between actual and present estimated values of (5^,.)^, and 

Qnp is made. 


TABLE 5.11 

COMPARISON OF ESTIMATED VARIABLES OF FINAL BREACH 
GEOMETRY FOR TETON DAM BREACH 


Mqiiation Used 

Qm 

a, 

^in 

Qpsn 

Observed 

10,421.6 

6549.87 

7550.89 

9.133 

W.ith Lab Data Regression 

4,688.99 

3834.57 

4138.11 

6.2 

Lab+Teton Data Regression 

6,576.9 

4173.8 

5350.54 

8.8 

Modified 

10,301.2 

5,984.2 

7432.3 

9.06 




OSLSJ9U9 Present dot 

Y = 0.035*(Vsrm)”' 
Wtiere , 


Mg. 5.9(a): 


Variation o 
by average 
nurTiber(Vsf 


((Ao*(gD.r)°'9/(Qi„)ov = 0.0078 VVsrn)' 


Fitted curve (only for present data) 
ooooo Present data 
Teton dam data 


Vs«/(gHo)°-^ 

Fig. 5.9(b): Variation of nondimensional final depth of erosion 
at crest(Der) by average inflow discharge with 
resisting velocity number(VsRNj 
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TABLE 5.12 


FUNCTIONAL RELATIONS FOR FINAL BREACH GEOMETRY 
WITH DIFFERENT REGRESSION SCENAWO 



From Lab Data 

Lab data with Teton Dam 

Modified 


Only 

Data 


^/) ij §(^er) f 

a. 

^SRN ^ 0.01 

X 0.0336 

^SRN 0-035 

^ D V S^cr 

X 0.0035 

F,-^®'x 0.078 

F,-=/ X 0.008 

57” " 

X 0.0076 

X 0.0186 

^SRN ^ 0.02 

C>.,sw = 

CAPfj'^\335 

C/(P/'^^^x2.64 

C4P/^3.4 


The functional relationship evolved from regression analysis for the final morphology are as 
shown in Table 5.12. Three different regression scenarios belonging to: only with laboratory 
data, laboratory data and Teton dam data together, and the modified relationship are 
presented. Modified functional relationship is recommended for use. 

5.4 Breach Hydrograph 

The basic aim of this study is to know the variation of the breach discharge with time in 
addition to the morphological features of the dam breach. Keeping this in mind the breach 
outflow hydrograph is computed as stated earlier for each experimental run both for cohesive 
and cohesionless soil model. In those cases, dam geometry and inflow into the reservoir are 
varying along with the properties of dam fill materials. From the earlier analysis it is quite 
evident that the flow through the breach of the dam is controlled by breach morphology. 
Again the breach morphology is itself a function of dam geometry, inflow condition and the 
properties of the dam fill like soil resisting velocity (Vs^) and Shield’s entrainment number 
(r,.). But, the overall the trend of the flow hydrograph is more or less same; having a rising 
limb, peak and then a falling limb, though the shape of the rising and falling limb of the 
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hydrograph are different, depending on the above three influencing parameters. The rising 
limb reaches its peak very steeply whereas falling limb comes to a equilibrium level 
gradually. To fit a curve to the hydrograph and to generalise the equation to the fitted curv-e 
the following approach has been adopted. 

i. Time versus discharge hydrograph for all the runs are normalised with and Q\up as 
shown in Figure 5.10. 

ii. For the values of Q(/Qiyp, O.S, 0.7, 0.6and 0.4, values are obtained for rising limb and 
falling limb separately for all the runs of experiment conducted including cohesive and 



iii. From the observation it is found that this non-dimensional Qw^Qwp P^ots follows a 
curve of the nature: 

y = (5.12) 


where Y = 


Qw 


and X • 


WP 


wp 


The ‘m’ value is different for rising and falling (m^f) limb for a particular run and 
magnitudes of and mv^ydepend on inflow characteristics, dam geometry and the properties 

of dam fill materials. 

Figures 5.1 1(a) and (b) are developed as an aide to compute directly the magnitude of m,, 
and m^j- by knowing t/t^p and their corresponding Q»/Qfvp ■ These curves are developed firom 
Equation (5.12). These average values of and controls the shape and nature of flo 
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hydrograph. The average values of and are also obtained for two field data, Teton 
Dam and Huaceoto Dam. These average and n,„ values are found to be the fune’tion of 
r., and Fjm from regression analysis. From regression analysis the following relations are 
developed: 

=0.0415(r:f=xK«’ 

a„d m,, =0mS(V„y^‘' I 

The functional relations in Equation (5.13) are shown in Figures 5.12(a) and (b) respectively. 

For the validation of this methodology, two sample hydrograph of run no. 2703 ( dam 
model with cohesive soil) and run no. 1 005 (dam model with cohesionless soil) are tested with 
the observed hydrograph as shown in Figures 5.13(a) and (b). The agreement is encouraging. 
Applying the same approach to the observed outflow hydrograph of Teton dam and Huaceoto 
dam breach shown in Figures 5.14(a) and (b), m,,. and are computed using Equation 
(5.13). These values arc listed below. 

Average Average n^f 
'I'cton Dam 3 .'m 3.03 

Huaceoto Dam 6.00 7.25 

Comparison of normalised outflow hydrograph predicted from the present analysis with 
observed data and BEED model for Teton Dam is carried out in Figure 5.15(a). It may be 
observed that predicted rising limb from the present analysis is not in close agreement with 
observed data. The falling limb of the present analysis predicts fairly close to the observed 
value. Further it may be observed that the BEED model data predicts rising limb in close 
agreement with observed data whereas its falling limb deviates considerably. 

Similar comparison for nomialised outflow hydrograph for Huaceoto dam is shown in 
Figure 5.15(b). It may be observed that the predicted rising limb deviates slightly from both 
observed and BliED model data, whereas the falling limb agrees fairly close with observed 
and BEED model data. In the Huaceoto dam out flow hydrograph analysis, the value of Ap 
used in the predicted curve is taken from the observed value. 

From the validation of Q^’|> with time in the present approach, nonconformity of the 
present result in rising limb may be observed. This is because of the nonavailibility of more 
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Fig. 5.1.5(a): Comparison of observed normalised hydrograph 
with its estimated value for run no. 2204, dam 
model with cohesive soil 



Fig. 5.13(b): Comparison of observed normalised hydrograph 
with its estimated value for run no. 1005, dam 
model with cohesionless soil 
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field data for prediction of exponent appropriately. Considering the non-availability of 
proper data for rising limb of the outflow hydrograph, present result may be accepted with 
some reservation. 



Fig. 5.15(q): C°mparisOT of 

Tor Teton Tom breach by overtopping flow 
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5.5 Breach Sediment Discharge With Time 

The computation of sediment discharge graph due to the dam breach erosion is essential 
along with water discharge for their downstream routing. The variation of sediment discharge 
curve with time have very steep rising limb and steep falling limb in cohesionless soil dam 
model whereas in cohesive soil models, the rising limb and falling limb of sediment discharge 
with time are gradual. Analysis of sediment discharge variation with time is carried out on 
similar lines of breach out flow hydrograph. 

Adopted procedure to compute ; 

a. First, time versus Qs graphs for different experimental data are normalised with t,p and O^p 
respectively. 

b. For the values ol (Js'' Qsp-> 0.6, and 0.4, t/t^p values are obtained for both rising and 

falling limb separately for all the runs conducted in the laboratory. 

c. Knowing !/r,p values for different Qsy^Qsp i-e. at 0.8, 0.7, 0.6 and 0.4 the corresponding ‘m’ 
values i.e., m„. (for rising limb) (falling limb) are obtained using Equation (5.12). Then 
the average and m,., for each set is used for further analysis. 

The /«„. and values controls the shape of sediment discharge graph with time. There is 
no field data information on sediment discharge with time. So only lab data are used to 
develop the functional relation among t»c, Qpn, CAPf^, as shown in Figure 5.16(a).The 
relation among them can be written as: 

= 121517^33 ( 5 - 14 ) 

where x rjf x x CAPy] 

As in most of the experimental runs the rising limb in sediment discharge curve with time is 
steep, the numerical values of w,, are almost same. So instead of developing a functional 
relation among its controlling dam breach parameters, the average value of all for different 
mns are taken. This average value works out to be 2.57. Here with caution, it may be noted 
that this recommended value oi ' based on limited experimental observations. 

To test this approach, the sediment discharge of Teton dam breach as shown in Figure 
5.16(b), simulated by BI-ED model is considered. From Equation 5.14, m,f value for Teton 
Dam is equal to 8.96 and assuming rrisr « 2.57, the non-dimensional plot of t/t,p versus Qs^Qs? 
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is plotted in Figure 5.16(c). The nre-^f^nt U • 

BEED model for Teton dam breach The a^^ simulated result of 



Fig. 5.16(c): Comparison of present estimated normalised 

simulated result 

or dLl-D mod6l for Teton dom 
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S£ Suggested Procedure to Compute the Dam Breach Process 
From the analysis and discussion in the previous chapter and present chapter, it is understood 
that the dam breach process can be seen from three broad angles, i.e. breach morphology, 
breach hydraulics and breach erosion, for a given condition. All of them are implicitiv related. 
To express the breach criterion variables like water discharge and sediment discharge, some 
implicit functional relations are developed from the observed data. This methodolosv can be 
used for any given field situation with fair degree of prediction of the breach characteristics 
due to overtopping flood. This algorithm can suitably be converted to a computer code on 
dam breach analysis. 

a. Methodology for Computation of Parameters Governing the Breach Process 

i. From the given predictor variables like capacity of the reservoir, Q.^ , Hp, Aq, and 

dam fill properties c,(j}, and the following non-dimensional parameters can be found 
out: CAPfj , . 

ii. Peak breach discharge (Qfyp) and breach duration(tp 5 ) can be obtained from Equation (5.2), 
Figure 5.3. 

iii. Predicting Qiyp the following nondimensional numbers: Qw? Qapn^ Qpsn can be 
computed. 

iv. The rj can be computed knowing from Equation (4.20b), Figure 4.38(c). 

V. The reservoir level corresponding to h^p , can be computed knowing from Figure 
4.38(d), Equation (4.20c). 

vi. Reservoir level at any time of breach, (A„,^ ), can be obtained knowing %, Qps^ and 
t* corresponding to front Equation 4.17b, Figure 4.34(b). Reservoir time scale, t* can be 
obtained from Figure 4.34(d), Equation (4.17c). 

vii. Knowing h^p, instant of breach can be found out from 

Figure 4.40(b), Equation (4.20c) 

viii. Knowing and ) at any time of breach, q can be obtained for from 

Figure 4.27(b). 



ix. Knowing and c, can be computed from Figure 4. 
X. Water level at crest can be obtained 
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27(a), Equatioa(4.13a) 
assuming critical state of flow there 








xi. Knowing at any time, drop of water level at crest with respect to n , 

So,S„canbeobtainedfromFigure4.29(a),Eqnation(4.I5a) “ =“cula,ed 

xii. So, Q^y 

xm. Knowing h,„, and , ^can be found 


xiv. Knowing and 

(4.18b) 


out from Figure 4.36, Equation (4.19a). 

^9i, t„ be found out from Figure 4.35(c), Equation 


* z, r„ ,, , ^ 

By can fairly be taken as equal to Bo. 


computed applying continuity 


xvii. From estimated h„ AH. and g,, v«," can be 
equation, (as, K,-AH=h"^). 

xviii. Ej can be computed after the estimation of VW and h \ 

xix. Knowing and AH. i,, at any time (V.J., can be obtained from Figure 4 39 

Equation (4 .21) 

XX. From estimated K, g,’ can be computed from Equation (4.22), Figure 4.40, Qs for given 
Yb can be computed. 

xxi. With estimated h AH V at anv inctcnt r , 

icu rty, mi, K.at any mstant > C = ~ can be computed from Figure 4.41, 

Equation (4.23c) . So can be obtained knowing the corresponding Qjy. 

xxii. The average of above two values of & obtained from previous two steps can be taken 
for further computation . 
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xxiii. With ^95> > time of occurrence of peak breach discharge can be 

computed from Equation (5.3), Figure 5.4. 

xxiv. Final top breach width (B„Vcan be obtained from Figure 5.5. Equation (5.6) for known 

QpNi 

XXV. Final depth of erosion (Z)J at crest can be calculated from Figure 5.6, Equation (5.8) 
after computing (5^^)/- 

xxvi. With estimated value of terminal top width erosion(5,,)^ and depth erosion(D,,), the 
side slope(z) of final breach cross section can be obtained from Figure 5.7, Equation (5.9) 

xxvi. Peak sediment discharge intensity( 9 ,p) can be estimated using Equation (5.10c), Figure 
5.8(a) for a known Qwp, 

xxvii. Time of occurrence of peak sediment discharge intensity {t,p) can be obtained from 
Figure 5.8 (b). Equation (5.1 Ic), once t^p is computed. 

With the above procedure, it is possible to predict the following variables at any instant of 
breach. 

a- 

b. AH, AE 

c. V, 

d. Qw and Qs 

To supplement these results particularly on final breach morphological variables like 

and peak breach discharge another approach has been suggested, taking into consideration 
the average inflow discharge to the reservoir and capacity of the reservoir as follows: 

1. Final top breach width at the crest {B^^fCdca be obtained from Figure 5.9(a) for a given Aq 
and For its better estimation the average of present estimation and the estimated 
value in step xxv, can be taken. 

2. For a given value of the capacity of the reservoir, 2,n j and Qpsn, 1^6 peak outflow can 
be obtained from Figure 5.9(d). The average of the value of Q^yp estimated in step ii and 
present computation can be taken for further use . 
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3. From Figure 5.9(b). Table 5.12, can be computed for a given The average 

of Step xxvi and its present value can be taken. 

4. The final area of erosion at crest can be obtained from Figure 5.9(c), Table 5.12 so r can 

be obtained for(5^,.j^ and D^. 

b. Method to Compute Breach Hydrograph and Breach Sediment Discharge graph 

In the section 5.6a, a detailed methodology to compute outflow discharge and sediment 
discharge with time in the breach knowing dam geometry, and soil properties of embankment 
fill is given. In addition to this, a method based on the shape of hydrograph given in Equation 
(5.12) is suggested in this section. This method is tested for two field data namely Teton dam 
breach and Huaccoto dam breach. 

i. For a given value of can be computed from Figures 5.12(a) and (b) 

respectively. The Equations (5.12) can be referred. 

ii. Using the earlier estimation of Q^yp and t^p in steps ii and xxiii respectively , and 
the breach water discharge with time can be obtained. 

iii. The value of m^j- can be obtained from Figure 5.14(a) , Equation (5.14). For given 
CAPf^ and computed The value of m^^can be assumed to lie within 2.5 to 3.0. With 
the earlier estimation of q^p and t^p in steps xxvii and xxviii, the sediment discharge at any 
time can be computed . 

c. Comparison of Overall Estimated Breach Parameters 

The final breach dimensions are depth of erosion Dg^ width of erosion Bg^ and side slope of 
the bank (z) at the crest of the breach profile. The duration of the breach process is considered 
as the time taken from the start of erosion till 95 percent of sediment is eroded(tp 5 ). In the 
breach process the controlling flow parameters are the magnitude of peak water discharge 
(Qwp), sediment peak discharge intensity {q,p) and their corresponding time of occurrence 
namely t^p and t,p respectively. These parameters are computed for Teton dam and Huaccoto 
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dam data using the procedure suggested in sections 5.6a and 5.6b and are presented in Table 
5.13. 


TABLE 5.13 


COMP.AJRISON OF ESTIMATED RESLTTS FOR 
TETON AND HUACCOTO D.AM BRE.\CH 



This result is presented in a form of bar charts for better overall comparison of estimated 
results with observed and SEED simulated values for Teton dam in Figures 5.17(a) and (b) 
and for Huaccoto dam in Figure 5.17(c). 

d. Flow Chart for Computation of Breach Process 

The procedure suggested in sections 5.6a and 5.6b for the estimation of variables controlling 
the breach process is presented in a logical sequence in a form of flow chart shown in Figure 
5.18. This flow chart suggest a method of computation in a sequence given as reservoir 
level(/ 2 ,,^), crest height of sediment bed level(; 2 <.„), head of water acting over the crest (A„), 
breach width at the crest(5 breach water discharge(2^), drop in water level {AH), energy 
loss (zl£) , average shear velocity(F.) and sediment discharge {Q^. 








Estimated Observed 


/ 1 

Estimated Observed SEED 




8-5.1 7(b). Comparison of estimated final breach dimensions for Teton dam breach with 
Observed and simulated value of SEED model 











Fig. 5.18; Flow chart for computation of breach criterion variables during breach process 
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5.7 Discussion of Results and Method of Computation 

a. Discussion of Results 

Functional relations established through the comparison with two field events can not be 
considered as final. Present analysis does not involve large number of laboratory tests on 
different types of soil and geometry of embankments. Only two field data are incorporated in 
recommending some relations. Still the present study predicts the available field data quite 
closely. It may be emphasized that this problem is governed by flow mechanism which is 
controlled by some major parameters like Vsj^, CAP^ and r^. Though 

there may be slight variation in the relation for prediction of output variables, w'hen different 
combinations of laboratory data and field data are used in regression analysis, the basic trends 
do not vary much. 

It may be admitted here that these experiments are conducted for constant inflow into the 
reservoir (2,„), which may not be tme in case of field. So the experiments with temporal 
variation of C/n are necessary to strengthen these findings. Also the experiments are carried 
out with homogeneous fill embankment dam, which is not necessarily tme for actual dams. 
Keeping these two limitations in view, one may consider the results obtained in the present 
study which can be treated as a fairly good step towards the prediction purpose. 

b. Comments on the Procedure to Compute Dam Breach Process 

As soon as overtopping of flow occurs, erosion of sediment in the embankment sets in. 
Erosion of sediment leads to the fomiation of breach. The breach depth, breach width, side 
slope of breach banks and longitudinal bed slope vary as the breach processes continue. The 
variations of these breach dimensions depend upon the difference in water level upstream and 
downstream of the embankment. Breach dimensions, breach erosion of sediment and outflow 
in breach are related with this difference in water level AH which is function of time. The 
magnitude of variation with time is related with functional parameters involving embankment 
dimensions, inflow to the reservoir and its capacity and erosion resisting properties of 
embankment soil. A methodology to compute this process is described in section 5.6a. All the 
dimensions, quantity of water flow and volume of sediment erosion in a breach can be 
predicted from the initiation of breach process till it ends. 
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Another approach to predict breach outflow discharge and sedin^ent discharge is carried by 
nondimensionahsing hydrograph and sediment discharge graph by their peak values and their 
tune of occurrence in the section 5.6b. These peak values and their corresponding time of 
occurrence can be predicted using the procedure suggested in section 5.6a. The shape of 
nondimensionalised hydrograph and sediment discharge graph are represented in the form of 
Equation (5.12), the magnitude of the exponent ‘m’ is correlated with breach parameters. In 
this method the shapes of the hydrograph are the important controlling parameters i.e. the 
magnitude of and which indirectly are controlled by the breach mechanism. In 

other words, all the mechanism in the breach process is lumped in the values and 

One has to be very careful in predicting these parameters since the observed hydrograph is 
very much sensitive to these parameters. Though in the present study, prediction of the falling 
limb of water and sediment discharge graph with time are done fairly well, more numbers of 
field and experimental observation are required to strengthen the functional relationships of 

m^,f, m^fand . 

c. Limitations of the Present Approach 

i. In the comparison of the result, it may be noticed that the functional relationships 
particularly for terminal breach morphological features are developed using the present lab 
data and the data of Teton dam and Huaccoto dam. The same functional relations are again 
verified on these two field events, which is not a correct approach. These functional 
relations should have been tested in some other field events. But due to the lack of required 
information about other dam breach cases, this has been done specially. However, it may 
be stressed upon that the functional relations only for lab data come very close to the mixed 
data of both lab and field in most of the cases. The developed functional relations among 
the variables of dam breach will be strengthened after testing in some more field events. 

ii. Calculation of velocity at the downstream end of the dam (V^"): 

The width of the breach at downstream end is taken equal to the width of breach at a 
section slightly upstream of it. This needs verification. 

iii. Computation of sediment concentration (C) or sediment load (g * ): 

As shown in Figure 4.41 and Figure 4.40 the correlation is not in good agreement with the 
breach parameters for C or g* respectively. 
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iv. Average inflow discharge into the reservoir: 

As the total inflow hydrograph for Teton dam was not available, Q, has been taken equal to 

half of the maximum inflow discharge. This seems to be arbitrao. although this assumption 
predicted value fairly close to the observed one, it can be suggested to take; 

Q _ Total in flow volume of water 
Time base of hydrograph 

The importance of this is felt, as it is required first to compute peak outflow discharge 

iQwp) and subsequently other related breach parameters. In the present experimental study, 
Qi„ was constant. 

5.8 Concluding Remarks 

From the study on the comparison of results, it is found that the present experimental 
observation of embankment dam breaching gives acceptable relations among the dam breach 
parameters, though it is tested for only two historical cases. This study helps in predicting the 
final values of the breach process of prototype dam failure close to their observed values. It is 
believed that present study will help in the development of numerical models more 
appropriately in simulating the breach process almost near to the field observation. 


PART II 

WASHOUT PROCESS OF FUSE PLUGS 



PART II 


WASHOUT PROCESS OF FUSE PLUGS 


Breach formation due ,o overtopping of flow is a process it. which bed. lateral and 
vertical eroston dominate. In washout of a fcse plug, lateral erosion is restricted because 
of solid boundaty walls and mostly vertical erosion process prevails all over its width 


Experiments am earned out to measure the water level variation and sediment bed level 
variation during the washout process. Method to compute the breach outflow hydrograph, 

volume of sediment eroded and flow parameters associated with washout process are 
presented in Chapter VI. 

Analysis of outflow discharge and sediment discharge over the washout sediment 
beds of fuse plug is presented in Chapter VII. Energy loss associated with sediment 
washout and transportation is analysed using data of mobile bed and its analogous rigid 
bed sediment profiles. A method to predict the washout process of fuse plug is also 


suggested. 



CHAPTER VI 


SCHEME OF EXPERIMENTS ON FUSE PLUG 
WASHOUT 


Fuse plug can be treated similar to an emergency spillway. This comes into action when the 
reservoir level exceeds the assigned maximum water level. Fuse plugs are made up of soil 
embankments. The location of fuse plug depends on the topography along longitudinal axis 
of the dam for disposing safely the washed out soil and excess water downstream without 
endangering the dam body. The bottom level of the fuse plug is placed at full reservoir level 
(FRL) snd its top level is at maximum reservoir level (MRL). Length of the fuse plug depends 
on the design criteria for releasing the excess water to be discharged downstream. Since the 
top level of the fuse plug is well below the embankment level, the excess water spreads 
completely across the length of the fuse plug. This cause the erosion of fuse plug material all 
over the lateral length of the fuse plug. Washout process in the fuse plug can be treated as 
two-dimensional erosion process in which erosion in vertical downward direction and in 
longitudinal direction takes place freely. The erosion in lateral direction is confined as it is not 
allowed to spread to the main dam body. Soil fill in the fuse plug is meant to be washed out 
when flow overtops. Hence, the erosion process in the fuse plug is called as ‘washout 
process’. 

The present study on the washout behaviour of fuse-plug due to overtopping is essentially 
a laboratory scaled model study. All the experiments were carried out in the Hydraulic 
laboratory of Indian Institute of Technology, Kanpur. The major objective of this 
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Experiments were carried out in two phases. In the first phase of experiments the mechanism 
involved in washout process of soil due to overflow on fuse-plug was planned to study. In the 
second phase of experiments sediment bed washout profiles are modeled as rigid bed profiles. 
On these analogous rigid bed profiles flow characteristics corresponding to discharge and 
energy loss are studied. 

In the present chapter, details of flume, fuse-plug model, soil used for model preparation, 
basic data measurement and the procedure to compute basic variables involved in the problem 
are presented. 

6.1 Details of Flume and Fuse Plug Model 

The flume which is used for embankment breach study is also used for the washout study for 
fuse plug. Only the position of test section was varied. The test section for all experiments 
was maintained at 1 5m from the inlet. 

a. Fuse Plug Model Preparation 

The model of earthen dam was made in wood with a rectangular breach section at the centre. 
This breach section was filled with compacted soil and is called as the fuse-plug on the dam 
body. The fuse plug section had the dimensions as crest width (B^) of 0.15m, bottom width 
(Bfb) of 0.875m and height (Hq) of 0.145 m with upstream and downstream face slopes as 1 
vertical to 2.5 horizontal. The central breach cut was a rectangular section of width 20cm in 
the centre. Figure 6.1(a) shows the detailed plan and sections of the fuse plug. Figure 6.1(b) 
shows the photographs of the fuse plug model used in the experiment. 

In all the models of fuse-plug experiments the following geometry of the model were kept 
constant. 

= Height of fuse plug = 0.145m 
Bf= width of fuse-plug = 0.20 m 

Slopes of upstream and downstream faces of the fuse plug = 1(V) ; 2.5(H) 



14.5cn 
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PLAN VIEW UF THE VODDEN MODEL SECTIONAL VIEW 

Fig. 6,li FUSE PLUG MODEL IN PLAN AND SECTIONAL VIEW 
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Fig. 6.1: Photograph of fuse plug model 
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b. Source of Soil used for the Preparation of Fuse Plug 

In the experiments of phase- A, cohesionless sand was used as flise plug fill material. This 
sand was collected from the river Ganga near Kanpur. Details of the sand are already given in 
chapter III in connection with embankment breach study. Since experimental data based on 
the experiments conducted by Rajpal (1995) are used extensively in this study, it is proper to 
describe the soil used by him for his experiments. He conducted the experiments on three 
different soil fill in the breach section of fuse-plug model namely, i. fine cohesionless sand 
(type-A soil), ii. Cohesive soil (type-B soil), collected from L I. T., Kanpur campus, and iii. 
mix of type-A and type-B soil in the proportion, 3:2. The particle size distribution cur\'es for 
these soils are shown in Figure 6.2, and their other properties are presented in Table 6.1. In 
the second phase of experiments the fine sand from river Ganga at Kanpur and coarse sand 
from river Yamuna at Kalpi were used for model preparation. The location of Kalpi is 100 km 
south of Kanpur. The sand from this river is very coarse. The size distributions for these sands 
are shown in Figure 6.2. The particles are well graded, relatively coarse and very angular. 


TABLE 6.1 

TYPES OF SOILS USED FOR THE FILL OF FUSE PLUG 


Types of soil 

dso 

(mm) 

< t ) 

(degree) 

C 

(KPA) 

Yb 

9810 * 

n / m ’ 

w 

% 

Remark 

A 

0.16 

32.1 

0.00 

1.58 

12 

Cohesionless fine sand 
collected from River Ganga 
near Kanpur. India 

B* 

0.28 

33.1 

1.60 

1.98 

14 

Cohesive soil collected from 
LIT, Kanpur campus 

* 

c 

0.18 

29.2 

0.85 

2.03 

15 

Mix of Type-A and Type-B soil 
in the proportion 3:2 

D 

0.60 

36.6 

0.00 

1.62 

12 

Coarse sand from River Yamuna 
at Kalpi 100km south of Kanpur 

E 

0.37 

35.2 

0.00 

1.65 

12 

Mix of Type-A and Type-D soil 
in the proportion 1 :1 


* Soil used by Rajpal (1995) for fuse plug fill 
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c. Preparation of Fuse Plug Model 

The wooden dam model block was fitted inside the flume at a distance of 15 meters from the 
inlet. Sufficient care was taken to prevent the flow of water beneath the dam model bottom 
and the flume bed by making it leak proof. The fuse plug portion of the dam model was open 
rectangular cut. This breach washout section was filled with wet soil with proper compaction, 
layer by layer. The wet soil was prepared by adding required quantity of water to dry soil for 
proper compaction. Each layer was of 5cm thickness. After compacting each layer, the top of 
it was scratched for next layer to adhere properly with the previous layer. The weight of dry 
soil and wet soil before filling the breach section and the left out excess wet soil after filling 
the breach section were measured for each experiment to find out bulk density and w'ater 
content. Soil samples with same bulk density and water content were tested in the Geotech 
laboratory of the institute. Tests were aimed to determine the angle of friction and cohesion of 
the samples. This was earned out because the fuse plug dimensions were small and any 
sampling of soil in it would have weakened the compaction. 

6.2 Experimental Procedure 

a. Aim of the Study 

Experiments were planned to study the erodible characteristics of sediment in the fuse plug in 
the first phase. Erodible characteristics are crest variation, outflow hydrograph, eroded 
sediment discharge graph, erosion rate, relative energy loss and discharge relations. In the 
second phase of experiments, discharge relation developed on eroded crest was compared 
with discharge characteristics of simulated rigid bed profiles. In addition to this, energy loss 
occuiring during the flow over the rigid bed profiles is computed. From these two studies, a 
relation on the magnitude of energy spent in eroding and transporting the sediment is 
postulated in terms of sediment concentration. 

Experiments for the present study were carried out under two phases. 

Phase-I: Erodible Fuse Plug Model 



Scheme of Experiments on Fuse Plug M^ashout 220 

This phase of study refers to the experiments which were conducted to know the mechanism 

of washout of soil in the fuse plug. The major objective of this part of study was the 
following: 

i. From the basic information on variations of water levels in the reservoir, over the fuse plug 
model and the sediment bed variation with time, the involved mechanism in soil-water 
interaction under highly dynamic condition like washout process of fuse plug be analysed. 
Phase-II: Analogous Rigid Bed Profile Study 

The second phase of experiments was carried out with the aim of having following major 
objectives. 

i. To verify outflow discharge relation for a erodible bed profile of constant width. 

ii. Energy loss characteristics during washout of the fuse plug. 

b. Basic Data Collected 

The basic data collected are : water surface variations in reservoir and water surface level and 
sediment bed level variations along the fuse-plug. 

i. Water Surface Variation 

Water surface variation in the reservoir and at 3-sections along the centerline of the fuse plug 
were measured with respect to time. Different persons measured the water surface levels using 
point gauge referring a common digital clock. The typical water surface profiles recorded are 
shown in Figure 6.3. It may be observed that the water surface level gradually increases, 
attains peak and starts decreasing rapidly during the washout of the fuse-plug. Later on, the 
variation in water surface level reaches gradually to a constant value, depending on the inflow 
condition. Erosion of the top surface of the fuse plug starts well before the occurrence of 
maximum water level in the reservoir. The water level variation at other sections indicate 
sudden drop followed by mild rises and then decreasing as the wash out process continues. 
These observations were seen in all the water surface profiles recorded at sections a., b, c and 
d as indicated in Figure 6.3. 

ii. Sediment Bed Profiles 

Sediment bed profiles were marked on both the non-erodible sides of fuse-plug during the 
wash out at different intervals of time. This marking was earned out using coloured pencils. 
The average profile of these two was considered as the profile acting along the center line of 
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This continues for quite sometime till the erosion on the upstream edge of the fuse plug get 
flatten as the crest starts moving upstream, simultaneously decreasing its height. A typical 
variation of the crest is shown m Figure 6.5. It may be noted that the erosion of the fuse plug 

occurs althrough the width of the fuse-plug and level of erosion is almost uniform across fuse- 
plug. 

6.3 Method of Computation of Basic Variables 

Among the computed variables, the water discharge {Q^) and the sediment discharge (^ 5 ) 
with time are the two most important variables which throw enough light on the wash out 
behaviour of a fuse plug. Due to their importance, detailed study was carried out. 

a. Breach Water Discharge {Qfy) 

The water discharge (Q^/) flowing over fuse-plug is computed from the governing equation of 
continuity as follows: 

Qw-Qin=^s^ ( 6 . 1 ) 

Where C,>, = incoming discharge to the reservoir, = Water surface area in the reservoir 

dh 

(constant in the present case as flume was prismatic), — rate of water level variation in 

dt 

the reservoir. In the above Equation ( 6 . 1 ), Q^ , are known quantities and the water level 
variation {dh^. I dt)m the reservoir can be obtained from the measured water surface {h^) 
which is a function of time. Thus the outflow discharge can be computed. A typical 
outflow discharge for a fuse plug is given in Figure 6 . 6 . 

b. Sediment Discharge (Qs) 

Volume of sediment eroded is computed by measuring the area between two consecutive 
sediment bed profiles. This area multiplied by the fuse-plug width gives the volume of 
sediment eroded during that time interval. Computation of sediment eroded at different 
intervals of time for consecutive sediment bed profiles is carried out and is plotted by adding 
in a cumulative way appropriately. Figure 6.7 represents a typical variation of cumulative 
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sediment volume with time. It may be observed that the increase in the cumulative volume 
^Vs is steep in the beginning and later on it remains constant. 



Fig. 6.5: A typical variation of crest of sediment bed profile 
hem (jo'" I'un no 1404) 






5] (Eroded sediment Volume) 


By taking the derivatives of the with time, as A V^/ At = Qs, the sediment discharge (Qs) 
is computed. The typical sediment discharge variation with time is plotted in Figure 6.8. The 
variation is asymmetric having a steep rising curve with gradual receding limb. 



Fig. 6.7: A typical variation of cumulative volume of 
eroded sediment from fuse plug body with 
time for run no 1404 



of sediment discharge with time 


Fig. 6.8; Variation 
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6.4 Details on the Number of Experiments Conducted 

In Table 6.2 number of experiments in different phases of fuse plug are listed. Details of 
expenments on washout process of fuse plug study are given in Table 6.3(a). In Table 6.3(b) 
the experimental information of Rajpal (1995) is presented. Table 6.4 contains the details of 
the experiments conducted mainly with a purpose to evaluate energy loss relations. The 
details of experiments on rigid bed profiles analogous to sediment bed profiles are listed in 
Table 6.5. 


TABLE 6.2 

NUMBER OF EXPERIMENTS CONDUCTED ON FUSE PLUG WASHOUT STUDIES 

IN DIFFERENT PHASES 


Types of 

Soil Used 

Phase A 

Phase B 

No of Experimental Data used in the 

Analysis 

Mobile Sediment bed Profile 

Analogous Rigid bed profile of Run No 1404 

A 

6 + 5* 

2 

6 

B 

5" 

- 

- 

C 

6' 

- 

- 

D 

- 

2 

- 

— ,, 

E 

- 

2 

- 


Soil type-E (Ganga Sand : Yamuna Sand= 1:1) 
Rajpal (1995)’s data 


TABLE 63a 

EXPERIMENTS CONDUCTED ON WASHOUT OF 
FUSE PLUG WITH TYPE - A SOIL FILL 


SNo 

Run No 

Qi. 

cm^/sec 

Soil Type 

Model 

Geometry 

Purpose 

1 

1204 

545.10 

A 

Ho= 0.145 m 


2 

1404 

5770.30 

A 

Bf = 0.2 m 

To study 

3 

j 

1804 

6531.30 I 

A 

Bfl,= 0.875 m 

washout 

•4 

1904 

5210.00 

A 

Bft = 0.15 m 

process of 

5 

2104 

5368.51 

A 

U/S and D/S slope= 1(V):2.5(H) 

fuse-plug 

6 

2204 

4210.95 

A 







table 63b 


EXPERIMENTS CONDUCTED BY RAJP.AL (1995) ON WASHOLT OF 
l^PLUG WITH THREE DIFFERENT SOH. FILLS: .4, B .4ND C 


SNo 

Run No 

Qi. 

cmVsec 

Soil Type 

Model 

Geometrv 

Purpose 

1 

Sl-I 

1373 

A 

Same 

To Study 

2 

Sl-2 

1035 

A 

as 

the Washout 

3 

Sl-3 

1514 

A 

mentioned 

process 

4 

Sl-4 

2247 

A 

in 

of Fuse Plug 

5 

Sl-5 

2320 

A 

Table 6.3 a 





1 

S2-1 

2801 

B 

Same 

To Study 

2 

S2-2 

5926 

B 

as 

the Washout 

3 

S2-3 

2962 

B 

mentioned 

process 

4 

S2-4 

5732 

B 

in 

of Fuse Plug 

5 

S2-5 

4367 

B 

Table 63 a 



1 

S3-1 

6060 

c 

Same 

To Study 

2 

S3-2 

5820 

c 

as 

the Washout 

3 

S3-3 

3800 

c 

mentioned 

process 

4 

S3-4 

7100 

c 

in 

of Fuse Plug 

5 

S3-5 

3850 

c 

Table 6.3a 


6 

S3-6 

5000 

c 




TABLE 6.4 

EXPERIMENTS CONDUCTED IN PHASE-B 


SNo 

Run No 

Qin 

cm^/sec 

Soil Type 

Model 

Geometry 

Purpose 

1 

■n 

1267.80 

A 


To study 



560.85 

A 

Same as 

energy loss 


1802 

942.50 

D 

Phase I 

characteristics 


0303 

904.02 

D 

experiments 

of washout of 


0503 

972.80 

E 


fuse-plug 


0603 

1139.54 

E 















table 6.5 


( IN PHA^-B^ SEDIMENT BED PROFILES 

; ^ -GOGOLS TO MOBILE BED PROFILES OF RLTy NO. 1404 


SI. No 

Run No. 


Qin 

(cm . see) 

hon 

Soil Used 

1 

2501 

0.387 

874.74 

1^4.6 




0.523 

1499.32 





0.681 

2395.29 


Tj-pe - .A 



0.935 

3621.02 




1.304 

7468.19 





2.064 

14535.33 



2 

3001 

0.350 

821.9 

5.0 




0.535 

1860.1 





0.600 

1932.6 





0.614 

1972.8 


Type- A 



1.083 

5580.0 




1.051 

8644.2 





1.319 

117806 



3 

0402 

0.195 

573.63 

7.5 




0.326 

1338.84 





0.485 

2530.16 





0.676 

4570.68 


Type - A 



0.837 

8129.83 




1.110 

10838.3 



4 

0702 

0.12 

385.64 

9.0 




0.20 

765.25 





0.23 

1125.44 





0.35 

1866.62 





0.40 

2851.02 


Type A 



0.42 

3048.53 





0.48 

4032.50 





0.54 

5007.03 





0.57 

5732.40 





0.69 

7333.95 



5 

0802 

0.041 

247.20 

13.5 




0.081 

637.68 


Type - A 



0.116 

1049.94 




0902 

0.047 

331.13 

14.1 




0.085 

484.88 


Type- A 



0.090 

528.13 













CHAPTER VII 


ANALYSIS AND RESULTS OF WASHOUT 
PROCESS OF FUSE PLUG 


The analysis of washout of fuse plug is carried out with the basic variables obtained from the 
experiments for different inflow scenario to the reservoir and variation in the fuse-plug fill 
properties. The analysis has been done under the following subheads: 

7. 1 . Reservoir time scale and length scale 

7.2. Outflow discharge relation 

7.3. Sediment bed profile variation 

7.4. Sediment discharge characteristics 

7.5. Erosion velocity 

7.6. Study on analogous rigid bed profile to instant erodible bed profile 

7.7. Prediction of Qif and with time 

Finally a procedure is suggested for the prediction of the variables controlling the washout 
process of fuse plug. 

7.1 Reservoir Time Scale and Length Scales 

In order to predict the washout characteristics of the fuse plug, the only observation that may 
be possible during the washout is the reservoir water level variation with time as shown in 
Figure 7. 1 (a) as a typical case. Considering this as the major information available to develop 



-y ff UJ i aoc //iWg 




scales for length and time, the gradient of water surface Mt is plotted against time as 
shown in Figure 7.1(b). It may be observed that the gradient has positive peak and negative 
peak crossing over the zero. The time taken between the zero value of the gradient to the 
maximum value of the negative gradient is considered as the time scale(r,). The magnitude 
of the difference in water levels, corresponding to these times is considered as the length 
scale(//*). The velocity scale is defined as the ratio of the length scale to the time scale. 
These scales are used for further analysis. In the analysis of fuse-plug characteristics flow and 
sediment movement variation, the length scale and time scale of the reserv'oir are used 
extensively. So it is necessary to relate these scales with incoming flow , geometry of the 
fuse-plug, and the soil fill properties. Based on these variables, inflow discharge intensity ( 9 ,„) 
is nondimensionalised with V^r and Ha as: . This parameter called as Qi„ 

number in short form, henceforth will be referred as {Qi^, used for analysis of 
nondimensionalised depth scale (//* and nondimensionalised time scale as 

shown in Figures 7.1(c) and (d) respectively. It may be observed that the depth scale increases 
with increase in , reaches a peak and then decreases with further increase in its value. The 
functional relation among them can be written as: 

— = 0.0937 -0.043Z + 0.142 -0.045Z' +0.005Z" -0.0001 (7.1a) 


Where X = 


Qin ^ 


The time scale in a nondimensional form with Vsr mdHa is found to be almost invariant with 
the Qji^i^ parameter and is of the order of 0.02. However very slight variation is observed and 
the relation from regression analysis as: 


^ H, ^ 
T.VsrJ 


= 0.000236 


q- ) 

V. H a J 


+ 0.02 


(7.1b) 


In Table 7.1 the time scale and length scales for different runs of experiments are shown 
along with other major variables in the washout process for the present study and for the 
experiments conducted by Rajpal(1995) on the same fuse plug model with three different fuse 
plug fills. These information are used in further analysis. 
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TABLE 7.1 

VAJUABLES USED FUSE PLUG ANALYSIS 
PRESENT STUDY (TYPE-A SOIL) 


Run No. 

Ho 

Q« 

Q., 

H. 

T. 




X:, 

present 

(cm) 

cm'scc 

cm'/scc 

cm 

see 

see 

cm- sec 

cm sec 


1204 

14,5 

545.10 

12,393.4 

l.S 

34.8 

45.0 

3.00 

0.310 

7.03 

1404 

14.5 

5770.30 

15,528.0 

2.1 

36.0 

40.0 

3.40 

0.232 

S.33 

1804 

14.5 

6531.30 

16,637.0 

-» 

45.0 

32.5 

4.15 

0.318 

9.09 

1904 

14.5 

1 

5210.00 ! 

14,758.0 

2.4 

40.8 

i 

32.5 

4.25 

0.360 

S.Ol 

2104 

14.5 

5368.51 

13,313.0 

2.9 

42.0 1 

31.0 

2.93 

0.260 

7 '*7 

2204 

14.5 

4210.95 j 

13,132.0 1 

3.4 

63.0 

44.0 

4.45 

0.224 

7.31 


RAJPAL’S (TYPE-A SOIL) 


Sl-1 

14.5 

1373 

9,600 

3.45 

56.0 

31.0 

2.50 

0.33 

5.29 

Sl-2 

14.5 

1035 

6,700 

2.75 

52.0 

47.0 

2.25 

0.20 

3.86 

Sl-3 

14.5 

1514 

8,200 

3.15 

65.0 

65.5 

3.65 

0.17 

4.72 

Sl-4 

14.5 

2247 

8,500 i 

3.40 , 

63.0 : 

60.0 

1.95 

0.18 

4.79 

Sl-5 

14.5 

2320 

8,700 

3.00 

55.0 

40.5 

2.20 

0.22 

4.95 


RAJPAL’S ( TYPE-B SOIL) 


S2-1 

14.5 

2801 

3,672 

4.50 

1032 

684 

0.45 

0.015 

0.88 

S2-2 

14.5 

5926 

7145 

6.35 

882 

771 

0.48 

0.020 

1.58 

S2-3 

14.5 

2962 

5140 

8.10 

852 

486 

0.44 

0.026 

1.22 

S2-4 

14.5 

5732 

7822 1 

6.70 

648 

498 

0.60 

0.017 

1.76 

82-5 

14.5 

4367 

5238 

7.00 

1980 

1200 

0.48 

0.007 

1.21 


RAJPAL’S (TYPE-C SOIL) 


S3-1 

14.5 

6060 

8,150.0 

5.2 

450 

261.0 

1.37 

0.0467 

2.85 

S3-2 

14.5 

5820 

10,170.0 

5.1 

213 

96.0 

1.14 

0.0760 

3.63 

S3-3 

14.5 

3800 

6,936.3 

5.5 

264 

184.8 

0.76 

0.0476 

2.55 

S3-4 

14.5 

7100 

9,539.0 

4.6 

333 

279.0 

1.38 

0.0400 

3.24 

S3-5 

14.5 

3850 

8,554.0 

5.9 

306 

171.0 

1.38 

0.0500 

3.15 

S3-6 

14.5 

5000 

7,613,0 

5.8 

456 

187.8 

1.30 1 

0.0600 

2.73 


where, Hd = 

Qin = 
Qwp “* 
tso = 

Vsv = 


Height of Fuse Plug 
Inflow to the Reservoir 
Maximum outflow 

time at which is eroded upto 50% of its height 
Sediment erosion velocity = Qs/A 


X 


22 




^wp 


sediment resisting velocity 
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7.2 Outflow Discharge Relation 

The outflow discharge computed from continuity equation is plotted against the head acting 
on the crest as shown in the Figure 7.2(a). This contains the data of present work and 

Rajpal s (1995) work for fuse plug model with cohesionless type-A soil. It may be observed 
that 2 k- varies almost linearly with instead of power law variation as in the case of weir 

flow. This may be due to rapid erosion of the crest during the wash out process. This analysis 
is further carried out by computing coefficient of discharge C(j as: 



(7.2a) 


The computed value of is plotted against 'y(^ as shown in Figure 7.2(b), where is the 
height of the crest above the bed level. It may be observ-ed that as increases. 


Cj decreases. The relation between them by regression analysis can be written as: 


Cj = 0.63 


u \'°- 

-I 

\Kn.) 


(7.2b) 


combining both the Equation (7.2a) and (7.2b), an equation for discharge can be w'ritten 
as: 


Q,.=0A2L-j2^h','h°J (7.2c) 

The outflow discharge depends on the crest height(/!„„) of erodible sediment bed profile 
along with the head of water over the crest as shown in Equation (7.2c). 



A A A A A P( 

• Pi 

Cc 


A 



I I I I I » > 1 

4.0 


Fig. 7.2(b); Variation of coefficient of discharge (c^) with 
ratio of head of water over erodibie sediment 
crest to the height of the erodibie crest 
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7.3 Erosional Characteristics of Sediment Crest and Bed Profiles 

The crest of the sediment bed, described in last paragraph is denoted as It controls the 
flow passing through the breach. Its location and height vary with the washout of fuse plug. 
Hence, it is necessary to know the crest location and its variation with respect to time. 

a. Crest Location 

The position of the crest of the sediment bed profile recedes from downstream top edge 
towards the upstream as the flow occurs over the crest. This is due to the erosion of the 
downstream face. The position of the crest is measured at various water levels over the crest 
and plotted as receding length versus drop in water level over the crest as shown in Figure 
7.3(a). The receding length, is measured as the distance from the upstream end of the fuse 
plug to the crest horizontally. Receding length decreases as drop increases. The receding 
length is nondimensionalised with base length {L^b) is nondimensionalised with Hq 

and is plotted as shown in Figure 7.3(b). The trend in variation is same as in dimensional case, 
the relation between them is written from regression analysis as: 

i=0.97-225Z+3.28Z--1.92Z' (7.3a) 

^DB 

where X = {prop in water level at crest / 77^) 

b. Sediment Bed Profile Downstream of Crest 

From the eroded sediment profile, it was observed that the upstream surface before the crest, 
varies very gradually or almost remains horizontal from the upstream slope to the crest. But 
the downstream surface varies steeply. The nondimensional downstream sediment surface 
profile in the form / h^^is plotted against {ix-x„)/{x^o - xj) as shown in Figure 7.3(c), 
where = height of the sediment bed at any location X, measured from upstream end of the 
dam, = height of the crest, x„ = the position at which hem occurs and x^q = the distance 

at which he = 05hem ■ Position oih,„, X^, X^o , and are shown in definition sketch. 

Figure 7.3(d). It was observed that he / hem decreases as X increases. The relation between 
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Fig. 7.3(c); Variation of nondimensional erodible sediment 
bed profile 
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them is given by analysis as; 


he 1 



Length scale of erosion surface profile is nondimensionalised with L^g and is plotted 
against t / T^Ai can be observed that (x^q / L^g) decreases with increase in i / as shown in 
Figure 7.j(e). The relationship between them can be written as: 



Fig. 7.3(e): Variation of nondimensional profile length scale (Xjo) with nondimensional 
time (t) 
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c. Crest Level Variation with Time 

In last paragraph, vanation of position of receding length is described. This is found to be a 
function of time. Figure 7.3(f ) shows the variation of nondimensional crest height 
with nondimensional time t/t^Q. Here is the time at which 50 percent of maximum crest 
height occurs. The data plotted in Figure 7.3(f) includes the data of the present study along 
with Rajpal s (1995) experiments. The relation between them may be wxitten as; 





1 + 




(7.3d) 


The variation of time scale is shown in Figure 7.3(g). Time scale is 
nondimensionalised with T* and is plotted against HyHm .The relation between them as 
obtained from regression analysis may be written as: 


= 0.296 
T, 


H, 


\ 0.633 


VHJ 


(7.3e) 


Here both scales T* and H* are the scales obtained from reservoir water level variation. 


d. Rate of Crest Erosion 

Rate of crest erosion / dt), is mainly controlled by the incoming flow velocity (V^- 

FFii ) properties of fill material of the fuse plug. The rate of crest erosion is plotted 

against incoming flow velocity over the crest. It may be observed in the Figure 7.3(h), that the 
rate of crest variation initially increases, attains peak and then decreases very steeply. This 
variation is observed in all the flow runs. Multiple values of rate of crest erosion for same 
value of incoming flow velocity exist near the peak region as can be observed from Figure 
7.3(h). 

Rate of crest erosion is non-dimensionalised with velocity of flow and is plotted 

against t / shown in Figure 7.3(i). It may be observed that the non-dimensional rate of 
erosion is increasing in the beginning of washout process till it reaches maximum value. 
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where, X=t/t5o 
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7.3(q): Variation of nondimensional_ time scale of 
erodible sediment bed profile. (tsp/T.) with 
reservoir length scaie(H.) and height of 
fuse plugCHoj 
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It decreases with further increase in f / The magnitude of maximum erosion coefficient is 
of the order of 2x10 . The functional relation among them can be written as; 

K = -0.0053 + 0.0243,r - 0.0264.V^ + O.OU.V^ - 0.0031.r'’ + 0,00034.V^ - l OT-lo"^ (7.30 

where Y = / dt) I and X = t 1 1 ,, 

This equation is valid upto t / t^a=3.5. 

The relation shown in Figure 7.3(i) does not contain sediment properties of the soil used 
in the fuse plug. To consider this aspect in the analysis of rate of crest erosion, the soil 
resisting velocity ) is incorporated in the nondimensional parameter A' 2 t , which involves 
soil and flow properties. The maximum value of dh^^dt is nondimensionalised with and 
is plotted against A '22 as shown in Figure 7.3(j). The relation between them is obtained by 
regression analysis as: 


M. 


V. 


SR 


where = 


= X 1-9 


(7.3g) 




A/9 C 

S‘^so 

Pw Pw 


It may be noted that the maximum rate of crest erosion varies as gj:” , at^d Hi 
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As the sediment resisting velocity increases magnitude of dh^„/dt decreases and discharge 
intensity of peak flow q^p increases as the magnitude of dh^„/dt increases. These observations 
sounds physically meaningful. 

7.4 Erosion Velocity 

Erosion Velocity is a measure of erosive characteristic of soil of fuse plug fill. Erosion 
velocity (Vsy) is computed as: 

y _ {Volume of se6\xnent eroded between two successive sediment profil^ 

'' ^ AtB^p(Ah^) 

where ~ average height of erosion between succesive sediment profiles. 

At = duration of time in seconds between two successive sediment profiles and Bjp = width of 
the fuse-plug. Erosion velocity V^y is plotted against velocity of flow at the crest = Fgh^ 
as shown in Figure lA(a). It may be observed that ^SV is maximum almost immediately at the 
starting of the erosion and then decreases as increases. There is a sudden decrease in Vsy 
when is of the order of 1 50 cm/sec. This is observed in almost all the runs in present data. 



Fig. 7.4(a): Variation of average sediment erosion velocity (r„) at any time of washout of 
fuse plug with incoming flow velocity (KJ 
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Erosion velocity is nondiraensionalised by its maximum value v,,,, and plotted against 
nondimensional incoming flow veIocity(K. ) with respect to ^ where H, is heigh, of 
fuse plug. 
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Fig. 7.4(b): Variation of nondimensional average sediment erosion velocity (Vsy) at any time 
of washout of fuse plug with nondimensional incoming flow velocity(F„) 


It may be observed that V^y ^'^SVM almost constant for major part of the erosion 

process. During the final stage of the erosion it decreases very sharply which can be seen from 
Figure 7.4(b). 
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Magnitude of erosion velocity scale is nondimensionalised with sediment resisting 

velocity Vs, and plotted against parameter Z,, as sho%vn in Figure 7.4(c). The parameter T,, 
is given as: 




22 



(7.4) 


Magnitude of Vsv\/^sr increases with increase in parameter .^ 22 - The relation among them 
is developed including Rajpal s data for three different types of fuse plug fills along with 
present data. The magjiitude of is of the order of 0.8 and 0.03 for cohesionless soil 

(type-A soil) and cohesive soil(type-C soil) respectively. 



Fig 7.4(c): Variation of nondimensional scale of sediment erosion velocity with 

maximum outflow discharge intensity, reservoir length scale (H. )and V,, 
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7.5 ScdimGnt Discharge Chsrsctcristics 

The rate of cumulative volume of sediment eroded is considered as the sediment discharge 
(a) . A typical variation of & with time is shown m Figure 7.5(a). It may be observed that 
the Qs has steep rising limb attaining a peak and decreases later on. Taking the magnitude of 
the peak value of sediment discharge ..its coiresponding time ^ the time at which 50% of 
Qsp occur on the receding limb noted as t,p/ 2 , the nondimensional plot of (Qs/Qsp) against 
[{t-tspVifsp.Q ~^sp)] is as shown in Figure 7.5(b). The data on the rising limb are very few. There 
is a clear variation of the receding limb as can be seen from Figure 7.5(b). Receding limb data 
is analysed further. The relation between them as shown in Figure 7.5(c) is obtained from 
analysis as: 

Qs _ 1 


a 


SP 




2 . 0 " 

1 + 

^ hp 1 



^ hp/2 ““ ^P ^ 



(7.5a) 


a. Sediment Concentration (C) 

Sediment concentration is obtained by dividing q, by q,,, where q^. is the water discharge 
intensity and q, is the sediment discharge intensity. They are plotted against ihyh,„,) as shown 
in Figure 7.5(d). It may be observed increases, sediment concentration decreases. 

This decreasing rate is related by regression analysis as: 


h 1 

h ] 

2.0 

(h ^ 

3.0 

] 

0.088 0.1134 “ +0.06 


-0.015 

u 

+ 0.001 

u 


U,J 






(7.5b) 


It may be noted here that the sediment concentration as high as 10%, is attained in the 
beginning of the washout period. Later on as washout continues, increases and 

sediment concentration decreases. 

The variation of is studied in Figure 7.5(e) with respect to drop in water level in 

the upstream and downstream of the fuse plug {M), which is nondimensionalised wdth Hp, and 
Qpn 2 - where Q/wi written as: 
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The relation is given using regression analysis as; 



9.9 1- 7.34 X 


Mi 


Hn 


^0 30g 
'pSl 


(7.5d) 



Fig. 7.5e ; Variation of the ratio of head of water on the crest(/zj over the sediment 
bed to its height(/2cm) with nondimensional drop of water AH. 


It is observed that sediment concentration qjq^ is the function of bed shear stress. With 
this in view, concentration (Q is plotted against V./Fsa indicated in Figure.7.5f. 
Concentration (C), varies very steeply in the beginning of the erosion process indicating the 
difficulty in predicting C for given values of This lead to correlate further with other 

parameters. Drop in the water level between downstream and upstream of the fuse plug, (an) 
seems to be a important variable during the overtopping of the fuse-plug. Using regression 
analysis the following functional relation among them as shown m Figure 7.5(g) has been 


established: 


il 


' ah' " 

UJ 


0.53 


X 



- 1.77 


*86.4 


J 


(7.5e) 


and the correlation coefficient « is 0.86. This indicates clearly, that the erosion process of the 
fuse-plug is a function of drop in the water surface level. 
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7.6 Analogous Rigid Bed Profile Study 


a. Discharge Characteristics 

The discharge and the flow depth relation obtained for the crest of erodible fuse plug sediment 
profile was found not to follow the weir type or spillway type relation. The value of exponent 
of the head acting over the erodible crest was of the order of unity instead of 3/2. This 
observation necessitated to verify on the rigid crest having the same sediment profile during 
the erosion. With this as aim, seven models of sediment profiles for the run no 1404, were 
built, one at a time, using cement and alluvial sand collected from river Ganga. This was done 
purposefully to maintain the same roughness characteristics of the mobile bed at a particular 
instant. For each rigid built-in sediment profile, experiments were carried out maintaining 
different conditions. The discharge and head relations are plotted as shown in Figure 

7.6(a). It may be noted that the discharge relationship follows weir relation as: 


0 ^ 

L 


= h 


1.497 


X 20.7 


(7.6a) 


where width of the fuse plug L = 20cm. From the Equation (7.6a), the coefficient of discharge 
works out as 0.707. Considering the contraction on either side of the rigid portion of the 
embankment into the analysis, the effective width is computed as: 


This length is used in analysis and the discharge relationship is written as. 


(7.6b) 


and c„ is found to be 0.71. From this analysis it may be noted that the flow characteristics 
over rigid bed crest follows exactly the weir relation. This leads to a question as to why the 
exponent of h. in case of mobile bed is not equal to 1.5. It may be postulated that the crest 
level goes on changing rapidly with time due to erosion, probably this may be the case for the 

deviation from weir formula. 
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b. Energy Loss Characteristics 

The energy levels at the upstream and the downstream of the rigid model are computed as Ej 
and .This energy ratio E/E, is plotted against as shown in Figure 7.6(c). As 
increases the ratio £ 2 /^/ increases. For higher value of h/h^^ this ratio reaches almost O.S. 
The relations among them for rigid profiles of cohesionless sand from river Ganga and its 
corresponding mobile bed and coarse sand from river Yamuna are shown in Figure 7.6(d), and 
these relations between them may be written in Equations (7.6d), (7.6e), and (7.6f) 
respectively. 

E^ f h '] ( h Y 

-— = 0.14 + 0.86 - -0.31 -~\ (Type- A soil, Rigid ofile) (7.6d) 

-^1 


E ( h ^ 

^ = 0.061 + 0.87 

^KmJ 



[Type- A soil, Mobile Prq/z/e) 


(7.6e) 


£ ( h \ 

-^ = -0.038 + 0.79 



(Type-D and E soil. Mobile Prq/z/e) 


(7.6f) 


Figure 7.6(d) reveals their inter-relationship. This curve relates the energy level at 
downstream due to flow over a rigid bed not involving any erosion or transport of sediment 
over the fuse-plug surface. It is of interest to know how much extra energy loss due to the 


same bed profile occurs if it is made mobile. For this condition energy level at upstream and 


downstream are plotted against h/h^^ for two different soils namely: type-A soil, fine sand 
collected from river Ganga and type-D soil, coarse sand collected from river Yamuna. It may 
be observed that the energy relation curve lies below the rigid bed and for Yamuna sand lies 
below the Ganga sand. This is due to the size variation of the sediment collected from those 
rivers. Median sediment size is 0.16mm for river Ganga sand and for river Yamuna sand it is 
0.6mm. This relation clearly shows the amount of energy spent in eroding and transporting it, 

depending on the size of the sediment as well as . 

For type-A soil, it was of interest to relate the magnitude of energy loss that occurred 
exclusively for sediment erosion and transportation. This is achieved by computing the 
difference of energy loss between the rigid and corresponding mobile bed of sediment profile 
for the same value of £, from Figure 7.6(d). ITus value is plotted to its coiresponding value of 


Qs^Qiy ■ 



Ez/E, 
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ond upstream of fuse plug to the ratio of head 
of water over the rigid crest of analogous 
erodible sediment bed to its crest height 



Fin 7 6fdV Energy loss characteristics in mobile and 
with different hg/hc 


•cm 
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The hyh 

cm is the link to get these values corresponding to each other as the experiments 
conducted in two different phases. The relation between them is plotted in Figure 7.6(e). It 
can be observed from Figure 7.6(e) that the energy loss / £, increases very steeply, the 

concentration increases from 0.01 to 0.05. Later on energy loss gradually increases with 
further increase in concentration. Here (tf£)«,v/ is defined as follows. 

It may be noted that 10% of the upstream energy (Ej) is spent when the sediment 
concentration is 0.14. This is the stage of hyperconcentration flow. This type of analysis could 
not be carried out for coarse sediment of river Yamuna sand because of nonavailability of 
concentration information for a given . 



Fig. 7.6(e): Variation of n 
mobile and its 


ondimensional magnitude of the difference in energy loss between 
analogous rigid sediment bed profile with concentration of 


sediment 
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Hydrograph and Sediment Discharge Graphs 

In fuse-plug washout study, in addition to the information regarding crest movement, crest 
variation and rate of crest variation, one will be interested to know the outflow discharge and 
sediment discharge at any instant of time. The determination of these quantities for any given 
flow situation for a given fuse plug with certain properties of fill materials will enable one to 
understand the mechanism involved in the fuse-plug washout process. Keeping this aspect in 
view both water discharge and sediment discharge at any time are analysed in the following 
paragraphs. 


a. Variation of Water Discharge Intensity with Time 

The discharge intensity over the fuse-plug is plotted against time as a sample case in 
Figure 7.7(a). It may be observed that increases, reaches a peak and then decreases 
gradually, reaches a state of equilibrium discharge intensity as the time increases. This 
relationship is observed in almost all the experimental runs. 

Using the peak value of g^^ and its corresponding time as scales, gjg.,vp is plotted 
against as shown in Figures. 7. 7(b), (c), and (d) for soil types A, B and C respectively. For 
soil type B and C the peak value g^p occurs almost at the end of sediment washout in fuse- 
plug. Hence, Figures 7.7(c) and (d) contain only rising limb of the variation of water 
discharge with time. The relation between gjg,^p and r/typ is written in the form: 




r<Y 




exp 


1 - 


frYl 




(7.7a) 


The value of 'n' is obtained by trial and error method of fitting curve. The values of exponent 


'n' are given below. 


Type of Soil Exponent 

(n) in eqn7.7a 


A 

B 

C 


2.50 

0.50 

0.45 



sec/cm) 
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Fig. 7.7(c): Variation of normalised outflow discharge 
intensity(qv,) with normalised time(t/Upj 
for fuseplug with type— B soil fill 



Fia 7 7(ci): Variation of normalised ouWow discharge 
^ ^ intensity(qw) with normalised _ tinne(t/twp) 
for fuseplug with type— C soil fill 
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It may be noted that the value of exponent 'n‘ for cohesionless soil (type- A) is very high (2,5) 

whereas for cohesive soil or for its mix is of the order of 0.5. This may be due to presence of 
cohesion, which offers more resistance to the erosion and delay the washout process. 

The scale is plotted against inflow discharge in a dimensionless form for all the three 
types of soils used in the present fuse-plug model as shown in Figure 7.7(e). It may be 
observed that q^^,p is related to as obtained from regression analysis: 


^WP 



X 0.0057 


(7.7b) 


To compute the outflow discharge at any time ‘t’, it is essential to know the time scale 
parameter t^^p. This time of occurrence of peak discharge is nondimensionalised with and 

capacity of the reservoir, plotted against nondimensional q^yp of the form Ho) ^ shown 
in Figure 7.7(f). The relation between them may be written as: 


-^ilidf_=o.384 - 0.0907 A + 0.007Z' - 0.0001 
Capacity 


where, X= 






wp 






(7.7c) 


An interesting observation may be marked in Figure 7.7(f). The time to peak discharge closely 
is related to the peak discharge and sediment resisting velocity for a given capacity of the 
reservoir, inflow intensity and height of dam. This value initially decreases, reaches a 
minimum of the order of 0.06 then gradually increases and as the tail variation trend shows it 
may remain constant. This study shows there is a mutual inter-relationship among 

C3.p3.Clty, q^ypy Hd VsR ■ 



q.p/(gHD ) 
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Fig. 7.7(e): Variation of nondimensional outflow discharge 
scale(_qv,p) with average inflow discharge 
intensity and soil resisting velocity, Vsr 



Fig. 


7 7(f)- Variation of nondimensipnalised 
nondimensionalised maximum 
discharge intensity 


time with 
outflow 
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b. Sediment Discharge Intensity with Time 

The sediment discharge per unit width(( 7 j) is plotted against time as can be seen in Figure 
7.7(g). The variation of with time follows similar trend of q.^^, increases initially reaches 
peak and gradually decreases with increasing time. Using the maximum value of sediment 
discharge intensity, {q^y, and its corresponding occurrence of time as scales, q/q^p is 
plotted against t/t^p. All the data in type- A soil model of the present investigation and data of 
Rajpal is incorporated in the Figure 7.7(h). An empirical relation among them may be written 
as : 


h. 




\^zpJ 


-(- 1 “ 


(7.7d) 


The agreement between them may be considered fair. 

The sediment discharge scales {q^^ is plotted with inflow discharge intensity {q,„) in a 
nondimensional form using height of the fuse plug {H^ and sediment resisting velocity 
as indicated in Figure 7.7(i). The relationship between them is obtained from regression 
analysis only for type-A soil. 


I sp 




f \ 

^ K J 


X 0.01 


(7.7e) 


In the computation of 2s at any time, the time scale parameter is essential to be 
estimated. From the regression study it is found that highly depends on t^p and q^.p. The 
ratio of time of occurrence of peak sediment discharge, t^p to the time of occurrence of peak 
water discharge, t^p taken as t^^/t^p is plotted against nondimensionalised peak water discharge 
q,,p as q,,p/{VsR Hq) in Figure 7.7(j). It may be observed that the time scale ratio t,/t,,p 
decreases with increase in intensity of peak water discharge or decrease in soil resistivity 
velocity. The relation between them is obtained as: 


1 



(7.7f) 
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Run no: 1804 (type—A soil) 
A typicoi sediment discharc 
intensity curve 


60 80 
Time (sec) 


Fig. 7.7(g): A typical variation of the sediment discharge 
intensity for run no 1804 with time 
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Fin 7 7(h): Variation of nondimensional segment discharge 
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Fig. 7.7(i): Variation of sediment discharge intensity 
scale(qsp) for type-A soil with peak 
outflow discharge(q„p) and Vsr 



7 7f;v I na in the time of occurrence of peak sediment 
Fig. 7.7(j). L pgQl, to outflow discharge 
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Co Time of Erosion {ig^) 

Sediment erosion in the final stage of the washout process of fuse-plug becomes almost zero 
and only the water will be flowing out. Time corresponding to 95 % of the erosion of sediment 
in the fuse-plug is used as erosion time and is denoted as tpj. In order to find £95 for any given 
situation, it is nondimensionalised with reservoir time scale (F*) and is plotted against 
along with V and Hq as indicated in Figure 7.7(k). The functional relation among them may 
be written as: 

^ = 2.55 -0.72Jr + 0.1^' -0.0013Z' (7.7g) 

J. ^ 



Fig. 7.7(k): Variation of nondimensional duration of washout process with average inOow 
discharge intensity and height of the fuse plug 

It is obsen-ed that the ratio % / T. decreases and reaches minimum of the order of 1.25 and 
increases further with increase in q,.. It may be mentioned here that increase in % may occur 
due to increase in q„ or due to decrease in Vs, for a given height of the dam. This observatton 

can be seen very clearly in the Figure 7.7(k). 
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7.8 Methodology to Predict Washout Process of Fuse Plua 

a. Procedure to predict Washout Variables 

i. From the known value of average inflow discharge intensity (^^), soil resisting velocitv' 
(Vsn), height of fuse-plug (Hg), the reservoir time scale (T. )and drop scale (H. )can be 
computed from Figures 7.1(c) and 7.1(d), Equations (7.1a) and (7.1b). 

ii. Knowing T, and H* and fuse-plug geometry, t^g, X^g can be computed from Figures 7.3(e) 
and 7.3(g), Equations (7.3c) and (7.3e ) respectively. 

iii. The duration of washout (fpj) can be obtained from Figure 7.7(k), Equation (7.7g). 

iv. The crest height of erodible sediment bed (Acm) during washout can be computed by 
knowing t^g from Figure 7.3(f), Equation (7.3d) and hence dhjdt can be obtained. 

V. The sediment bed height {h^ at any location {X) from the crest position can be obtained 
from Figure 7.3(c), Equation (7.3b). 

vi. From Figure 7.5(e), Equation (7.5d) head of water over the crest (A.J at any time of 
washout process can be computed knowing ho- 

vii. From the known values of and discharge coefficient Q can be obtained from 
Figure 7.2(b), Equation (7.2b) and Q;y from Equation (7.2c), Figure 7.2(a). 

viii. Assuming flow to be critical at the crest the water level at crest can be estimated for 
known outflow intensity and subsequently drop of water level at the crest can be 
computed. 

ix. The receding length {XJ of the crest of sediment bed from upstream face of the fuse plug 
can be obtained knowing drop of water level at crest from Figure 7.3(b), Equation (7.3a). 

X. The maximum value of sediment eroding velocity (VsyM) can be obtained from the mean 
curve of Figure 7.4(c) once q^p and are known. 

xi. The maximum value of outflow intensity can be obtained from known values of {q ^, ), 

VsR and Ho from Figure 7.7(e), Equation (7.7b). 

xii. The average velocity of eroded sediment can be obtained from Figure 7.4(b) for known 

K (FsK) washout process. 

xiii. The average eroded height of bed proffle(Ah.)at any time At can be obtained &om the 
predicted bed profile, h, ~Ax). Figure 7.3(c), Equation (7.3 b). 
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XIV. Assuming umform erosion along the wdth of fuse-plug, eroded sediment volume at 
any instant of washout process can be obtained kno^ving (KhJ, 

XV. So the sediment discharge, & = d(lVs )/dt can be estimated. 

XVI. From known value of hjh,„, C= q/q^ can be obtained from Figure 7.5(d), Equation 
(7.5b). Using the predicted value of q^ in step vii, q^ can be estimated. The average value 
of estimated ofq^ and its computed value in step xv can be used. 

xvii.The energy loss (^^{^ can be related to using the Figure 7.7(c). 

xviii.The drop of water level AH can be obtained for a given value of during the 

washout process from Figure 7.5(e), Equation ( 7 . 5 d). 

xix. From the computed values of AH, q/q,„ and known values of and Hp, average bed 
shear stress F. can be computed from Figure 7.5(g), Equation 7 . 5 (e). 

XX. Time of occurrence of peak water discharge {t^p) can be obtained from Figure 7.7(f), 
Equation (7.7c). 

xxi.Time of occurrence of sediment peak discharge {t,p) can be obtained from Figure 7.7(j) 
Equation (7.7f). 

b. Prediction of Water and Sediment Discharge During Washout Process 

Knowing < 7 ,,,^, t,„p, q^p and t^p the normalised hydrograph and sediment discharge graph can 
also be obtained by estimating the exponent ‘m’ in Equation (7.7a) properly. 

In the present method, no functional relation has been suggested to get these exponents 
for a given situation. For normalised hydrograph of type-A soil fill, m=2.S whereas for type - 
B soil fill m = 0.5, as shown in Figures 7.7(b), (c) and (d). 

For normalised sediment discharge graph for type-A soil fill m - 2.0 as shown in Figure. 
7.12 (h). Equation (7.7d). 

c. Flow Chart 

The steps mentioned in the previous section are represented in a more compact form m a flow 
chart. For given geometry of the fuse plug, the properties of its fill material, the reservoir 
capacity and inflow discharge characteristics, the suggested flow chart can be transformed to a 
computer code to know the various washout variables involved m the fuse plug at different 
time of washout. The time step can be chosen for the computation once the duration of 
washout is estimated with the established functional relation in the present study. 
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7.9 Conclusions and Limitations 

Fuse-plug washout and dam breach process by their inherent character look same outwardly. 
But in case of fuse-plug there is restriction on its possible lateral erosion. This causes the flow 
charactenstics to change. Contrary to dam breach, in fose-plug study the importance is laid on 
how quickly it gets washed out by releasing reservoir water without endangering the main 
structures like spillway and dam body. From the analysis of experimental data obtained from 
washout process of fuse-plug model studies with three different soil fills, np^-.4(cohesionIess) 
soil, (ype-5(cohesive) soil, and rype-C(mixture of cohesionless and cohesive soil in 60:40 
ratio)the following conclusion may be drawn. 

a. Conclusions 

i. Introduction of time scale and length scale of the reservoir may be very significant 
parameters to study the washout behaviour of the fuse-plug as it relates implicitly both the 
flow characteristics and properties of the fill material. 

ii. Crest height of the sediment bed profile during the washout and its rate of variation 
can be predicted during the washout of fuse plug which will help to predict the sediment 
bed profile, h^(X). 

iii. Prediction of water discharge passing through the fuse-plug during washout is 
attempted only in the rising part of washout hydrograph. From the study of the 

relation between and %, it may be observed that the washout process is almost in the 
final stages of erosion when water discharge attains its peak value. 

iv. Prediction of some of the key controlling parameters like, peak water discharge intensity 
{(Iwp )? t™® of occurrence of q^p and peak sediment discharge and its time of 
occuirence t^p is possible by the present approach. But the prediction of q^p demonstrated 
only for cohesionless (type-A soil) of a particular value of average grain size {d^o 
=0.1 6mm), which can not be generalised. 

b. Limitations 

The present study has its limitations also as described below. 

i. In the present analysis, geometry of the fhse-plng is maintained constant. Experiments with 
different geometry of the fuse-plug are needed for generalisation of the analysts. 
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ii. The prototype information available in literature could not be used in the present analysis 
because of the limitations on the available information. Hence the validity of present 
procedure needs to be tested with field data. 

Lastly it can be said that being handicapped with these two limitations, the analysis has not 
lost the rationality, rather it has a newness in its own form. 



CHAPTER VIII 


CONCLUSIONS AND RECOMMENDATIONS 


8.1 General 

Bitcicli foiTnstion in fl dsro End WEshout process of e fuse plug due to overtopping flow liEve 
been considered snd snElysed in the preceding chEpters. 

Erosion of sediment, chEnnelisEtion, crest movement, crest widening, sediment trEnsport 
and energy dissipation are the processes associated with breach formation. These processes 
have been experimentally studied with two types of soil embankment models, cohesionless 
and cohesive soils. Water surface levels, sediment bed levels, and width of breach at different 
sections along the breach in the direction of flow, have been continuously measured from the 
instant of start of breach till the stabilisation of the breach. Twelve and five runs have been 
carried out with cohesionless and cohesive soil embankment models, respectively. 

Analysis of data involves development of relations for flow over the breach, widening and 
deepening of the breach, sediment discharge and energy loss associated with the breach 
formation. Analysis of field data along with the experimental data from the present study 
leads to prediction of final breach dimensions and peak discharges associated with them. A 
method to predict the breach process has also been developed and verified with limited 
aviiilable field data. 

A weak granular embankment, which may or may not have a clay cover, is built along the 
axis of the dam body to act as automatic excess flood releasing device. The embankment is 



Conclusions and Recommendations 271 


designed ,o be washed ou, during high noods. Such an auxiliary which augments ,he 
semce spillway by passing large flood flows is called a fuse plug. Experiments axe conducted 
on the fuse plug models built of different types of soils. Water surface levels are measured 
continuously at different locations. Eroded Sediment beds are marked on the fuse plug side 
walls. Analysis of flow over the fuse plug, erosion of sediment and energy loss have been 
carried out. Discharge characteristics of non-erodible fuse plug profiles are measured and 
compared with erodible fuse plug profiles. 

Based on the analysis of experimental data, major conclusions for embankment dam 
breach and fuse plug have been presented under separate headings. 

8.2 Dam Breach Process 

Experimental data is analysed in terms of criterion variables and predictor variables. Predictor 
variables are dam geometry, capacity of reservoir and inflow discharge. The criterion 
variables are variations of breach geometry, outflow discharge, erosion rate of bed and width, 
sediment discharge and energy loss. Some of these functional relations have been modified 
suitably using available field data, so that the prototype breach process can be predicted with 
more confidence. Analysis of the results is presented under the headings namely, breach 
morphology, breach hydraulics, and breach erosion. 

a. Morphological Characteristics 

i. Photographs taken during breach process for cohesionless soil embankment show the 
different states of erosion, channelisation, collapsing of bank, crest formation, and 
channel widening. 

ii. Planfonn of the breach channel in cohesionless embankment and cohesive embankment 
has converging section, throat and diverging section. The position of the throat moves 
upstream as the process continues, this is evident from the photographs taken during the 
breach process. 

iii. The cross-sections of breach during erosion of cohesionless soil embankment is almost 
rectangular in shape, having steep banks, accompanied by bank collapsing. The cross 
sectional shape in cohesive soil embankment appears more like trapezoidal at the 
contraction section (upstream) and more like rectangular at the throat section. 
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IV. Formation of flow cascade and water falls were observed in cohesive soil embankment 
breach. In cohesionless soil embankment breach, these flow patterns were not present. 

b. Discharge Characteristics of Flow 

i. Flow in a cohesionless embankment breach behaves similar to flow over an erodible crest 
spillway. The alignment of spillway is straight in the beginning of the breach but as the 
breach progresses, it becomes curved with concaving in the direction of flow. This is 
clearly observed in photos of breach process. 

Flow in the case of a cohesive soil embankment breach is analogous to that over a broad 
crested weir, with a cascade type of flow on the downstream face in the beginning and 
later it is more like a water fall. 

ii. The discharge coefficient, c^is found to be varying during the breach process. It decreases 
with increase in the ratio of head over crest {h,) to the height of the crest (h^J. This 
variation continues till is equal to 0.4. For higher values of the magnitude of 

remains fairly constant. Rate of crest erosion is higher during hjh^ < 0.4 and 
gradually decrease as hjhc„ increases. This indicates that the coefficient of discharge q 
decreases during crest erosion. 

iii. Variation of reservoir level head acting over the crest of weir (/i„), drop in water 
level between reservoir and down stream of the dam {AH) and height of crest {h^^ ) with 
time from the start of the breach process are analysed combinedly for both types of soil 
embankments. The scales used for normalising these variables are functionally related 
with breach parametric numbers based on peak breach discharge {Qw^^ height of the 
dam(///;), cross-sectional area of the dam (Aq) and soil resisting velocity(F 5 ^). The breach 
parametric numbers are Qapn ^srn-- 

c. Erosional Characteristics of Breach 

i. Variations of crest height h^.^, crest width B^r and crest cross-sectional arez-A^^ are 
functionally related to either time variable or drop in water level at the crest, separately 
for cohesionless and cohesive soil embankments. Variations of crest height with drop 
in water level is carried out for combined data of cohesive and cohesionless soil 


embankment breach- 
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For cohesionl«s soil embaricment, ores, heigh, roeasured from the top of frte dam 
vanes w.ft K power law while the widtt, of erosion a, fte crest is found to vaty with a 
power of 0.57 of the drop in water level at the crest. 

For cohesive soil embankment, width of erosion a, the crest afrains maxinrum value a, 

the beginning of breach erosion and thereafter remains fairly constant at a value equal lo 
1 .2 times the height of dam (//j,). 

li. Rate of erosion of crest heigh, (i„), crest width (B„) and crest cross-sectional area(^„) 
are related to drop of water level, time and velocity of flow at crest. 

For cohesionless soil embankments, rate of crest height erosion is related to drop in 
water level at the crest. Rate of width erosion at crest is related to flow velocity at crest 
and time. 

For cohesive soil embankments, rates of crest height and crest width erosion varies 
with drop of water level between upstream and downstream of dam (AH). 

For combined data of cohesionless soil embankment and cohesive soil embankment, the 
rate of ciest height erosion is related to drop in the water level between reservoir level and 
downstream w'ater level (AH). 

Rate of erosion of crest height, crest width or cross-sectional area varies very steeply 
during initial state of erosion process. As the erosion process proceeds, the magnitude of 
these rates decreases. Contrary to the assertions as in numerical models of dam breach, the 
rates of erosion of the crest dimensions related with any of the flow parameters at the crest 
are found to vary with breach process. 

d. Flow Characteristics 

Mean velocity of the flow, computed from breach flow discharge and the average 
wetted flow cross-sectional area (Aj-) is averaged over the breach length at any instant of 
breach process for further computation of other breach variables. This average velocity, 
average hydraulic radius (RJ, energy gradient (Sj^) and the average velocity of the flow at 

the down stream top edge of the dam is related to, water surface slope, bed shear velocity and 
critical depth of flow for cohesionless embankment breach flow. 
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1. It IS observed from the analysis, that the mean flow velocity varies and 

individually. The applicability of Manning’s equation is tried for breach flow and it is 
found from regression analysis that relation exists considering both and together 

as: V^a . However, when is plotted against Manning’s relation S) , it 

is found to agree closely with regression relation stated above, 
ii. Flow velocity at the downstream top edge of the embankment when nondimensionalised 

with bed shear velocity( V, ) and also with critical velocity at this section -y' g>’^ , are found 

to vaiy with water surface gradient as 5°' and respectively. From these relations 

V . . . _ . 

it may be observed that —?== variation with is comparatively better correlated. 

4sy'c 

e. Sediment Transport Characteristics 

The sediment discharge intensity {q^ is related with flow discharge intensity bed shear 
velocity (F.), drop in water surface levels (AH) along with sediment resisting velocity (Vs^). 

i. For cohesionless sediment embankment breach q^ varies as 

ii. The ratio of cumulative volume of sediment eroded to the cumulative volume of water 
flowed is found to vary inversely with increase in the time of breach process. 

iii. The sediment discharge intensity {q,) is found to vary as square of the bed shear stressCtJ 

iv Sediment concentration defined as — is found to be of very high order, 8 to 10 percent 

< 1 ^ 

in the beginning of breach process, later on its magnitude decreases with increase in 
breach time or decrease in drop level (AH). It may be noted that when sediment 
concentration exceeds 8 percent, the flow can be treated as hyper concentrated flow. 

V. Sediment load g^* is found to increase with increase in square of the bed shear stress for 
combined data of cohesionless and cohesive soil embankment breach, 
vi. The relation between the sediment load (g, ) with ratio of bed shear velocity to sediment 

resisting velocity as f is found to vary with a quadratic polynomial relation obtained 

F CD 
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from regression analysis. This relation has a good agreement for higher values of — 

V 

^ SR 

However for lower values of ~ large scatter is found to exist. 

^SR 

vii. Sediment concentration for combined data is analysed, first carrying out regression 
analysis with sediment and embankment geometry parameters and then their relation is 

graphically represented. Regression relation of the data shows large scatter indicating a 
poor relationship. 


f. Energy Loss Characteristics 

Energy loss is defined as the difference in energy upstream of the embankment to the 
downstream of the embankment. It is related with water surface slope {SJ computed as — 

^DB 


A£ 

in the breach section. Relative energy loss — varies as 


V 


in the cohesionless soil 


embankment breach. For combined data, the quadratic polynomial relation found to agree 
closely in comparison to other relations. 


g. Final Breach Dimensions Based on Laboratory and Field Data 

Final depth and width of breach erosion, duration of the breach and peak flow discharge 
during breach are considered here as terminal breach parameters. These quantities are 
analysed using present laboratory data along with field data available in literature. Analysis is 
carried out like reservoir capacity, height of dam as known parameters. Since, the sediment 
properties and inflow discharge relations are not available for field data, they axe excluded in 
the analysis. 

i. Peak flow discharge is found to control the breach process and breach geometry. It is 
related to capacity of the reservoir {CAP) and height of the dam (iTo). From the analysis, 
it is found that the peak flow discharge increases with K power of the reservoir capacity 
and square root of the height of the dam. 

ii. Breach width is found to vary with peak flow discharge {Qwp) with a power 0.57 

along with reservoir capacity and height of the dam. 

iii. Breach depth of erosion (D J is found to be 0.74 times the height of the dam 
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IV. Duration of the breach tg;, nondimensionalised with depth of erosion {DJ is related to 
capacity of the reservoir, width of erosion (5J, peak flow discharge {Q„p) and height of 
the dam The analysis indicates a poor relationship between these parameters. This is, 
possibly due to lack of inclusion of embankment material properties. 

h. Major Criterion Variables Controlling the Breach Process 

The prediction of breach criterion variables like terminal breach ?eometrv i.e. (B. ),■ D Z 
and the peak outflow water discharge (Q^yp), sediment discharge (Qsp) and their time of 
occurrence, and t^p are very much essential to study the breach erosion process. To arrive at 
a functional relationship for these breach variables in the analysis, Uvo field data, Teton dam 
data and Huaccoto dam data are included along with the present experimental d a ta 

i. The duration of breach (%) varies with soil resisting velocity number as . 

ii. The breach peak discharge (Qiyp) varies with Shields entrainment number as and 
CAPg^ It may be noted that Qfyp is very much sensitive to dso as is a function of d;o 

iii. Time of occurrence of peak breach discharge (t^p) varies as Qp ^ , and 

iv. The terminal breach width at top of the crest (5^.^)/ varies as , Qfpl 

V. The ratio of terminal depth of breach to width of breach is related with 

and numbers and is found to vary as Q^n ^Vsr^ ■ This indicates that the 
geometry of the breach depends more on the soil properties than the breach peak flow. 

vi. The side slope of the breach cross section is found to be dependent on the final depth of 
erosion, width of erosion at crest, and height of the dam. 

The analysis of peak sediment discharge intensity and its time of occurrence(tjp) 

includes the simulated result of BEED model. 

vii. The breach sediment discharge intensity {q,p) is found to vary inversely with as 
F '7 This relation is obtained by non dimensionalising q^ by Qwp^Ho- The peak 
sedLent discharge occurs earlier to the breach peak water discharge in almost ail the 
experimental observations. 

viii. The ratio of to varies in a increasing trend rvith Q™- The time to peak sediment 
discharge is related to time of occurrence of peak flow discharge (t^). 
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ix. In the conclusions stated above from 1 m s in t 

to 8, the peak flow discharge is considered to be an 

important vanable in most of the cases. However, itself is a result of breach process. 

Hence, instead of Q,yp pother discharge parameter was essential to replace Q,yp. For this 

reason, the average inflow into the reservoir is considered to be an appropriate variable in 
further analysis. 

X. The final breach width depth of erosion {D^^) and eroded cross-sectional area j 

at crest in a nondimensional form with and Q. varies with soil resisting velocity as 
^SRN » ^SRN ’ ond respcctlvely. 


i. Methodology to Predict Breach Process 

Breach process involves the variations of breach geometry along with water and sediment 
discharges with time. Prediction of these information is earned out systematically based on 
the analysis of experimental results and also with two field observations. Procedure for 
prediction is presented in an algorithm and in a flow chart form. 

1 . An algorithm to predict reservoir water level, crest height and crest width during the breach 
process is presented. From these predicted values, computation of water discharge and 
sediment discharge graph during breach process are carried out. The distribution of 
normalised hydrograph and sediment discharge graph in the form of a equation 

Y = X" ' ^ where Y represents normalised water and sediment discharges and X 
represents normalised time. The exponent ‘m’ in the above relation is found to have 
different magnitudes for rising limb and falling limb. Magnitudes of these exponents are 
related to parameters derived from dam geometry, soil properties and peak discharges, 
using present laboratory data and one field data of Teton dam. The scales like peak 
discharges and their corresponding time of occurrence are computed from the relations 
developed. The variations of discharge of water and sediment for different time during the 
breach process are compared with observed laboratory data, observed field data for the 
Teton dam and Huaccoto dam and also with predicted BEED model result for these two 

dam breach failure cases. 

2. The breach dimensions like, final width of erosion, depth of erosion, side slope of the 
breach cross section and breach peak flow parameters like peak water discharge and peak 
sediment discharge and their corresponding time of occurrences and also duration of 
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breach are estimated using suggested methodology. These predicted results are compared 
with observed field data for Teton dam Huaccoto dam as shown in Table 8.1. Also these 
predicted results are compared with one numerical model namely BEED model. The 
agreement of the predicted results with observ'ed data is shown in the percentage of error. It 
may be observ'ed that the magmtude of percentage error for Teton dam is far less in 
comparison to Huaccoto dam. Here it should be noted that the inaccuracy in prediction of 
the breach variables for Huaccoto dam probably be due to the reported conflicting peak 
discharge values ranging from 10,000 to 18,000 m^'sec. 

TABLE 8.1 


COMP.4R1SON OF ESTIMATED RESLT.TS FOR TETON .AND HU.ACCOTO D.AMS 



Teton Dam 


Huaccoto 

Dam 


S.No. 

Vanables 

Estimated 

Observed 

BEED 

% Error 

Estimated 

Observed 

BEED 

% Error 

1 

Qwp 

67,028 

65,000 

68,581 

+3.12 






m^'sec 









2 

4sp X 1 0 

93.8 

- 

103.95 

9.75 

31.9 

- 

- 

- 


cmVsec/cm 









3 

^)5 

1.5 

1.5 

1.35 

0 

- 

- 

- 

- 


(hr) 









4 


1.27 

1.5 

1.35 

-15.3 

22.8 

17.8 


30.2 


(hr) 









5 

^sp 

1.32 


1.35 

+3.6 

4.45 





(hr) 









6 

{BJf 

212.5 

200 

295 

+6.3 

168 

200-300 

298.5 

22 


(m) 








-14.8 

7 

Der 

84.6 

79 

93 

+7.1 

91.1 

107 

85 


(m) 







1.67 

AS 

8 

Z 

0.57 

0.5 

0.675 

+10 

0.52 

1.0 


8.3 Washout Process in Fuse Plug 

The results derived from the analysis of water level and sediment bed level dunng washout of 
sediment fill in fuse plug, are summarised under the headings namely outflow discharge 
relation, washout profile characteristics, sediment transport charactenstics, energy loss 
characteristics, computation of outflow of water and sediment discharges. The experiments 
earned out using three types of soil as fill material in fuse plug. Five to six different inflow 
discharges were allowed into the reservoir for each type of soil fill. Experiments on simulated 
rigid washout profiles were carriedout to verify the discharge relatton. 
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a. Outflow Discharge Relation 

i. The discharge intensity (,„) over the erest of the washout sediment profile is found to 
vary linearly with head acting over ctest(ft.). The discharge coefficient ^computed based 
on weir relation is found to decrease with ih, / ft..) having a power of 0.2, where ft,, is 
height of sediment crest at any time of washout process. Discharge relation, considering 
the variation of is found to be related as; 

&= 0.421^2? ft" ft“ 

ii. Discharge characteristics over simulated rigid bed washout profiles were measured. They 
are found to follow weir type of relation having hi' and are independent of This 

leads to conclude that the variation in co-efficient of discharge {cj) is mainly due to 
sediment bed profile erosion during washout process. 

iii. Co-efficient of discharge for erodible fuse plug flows found to be function of the ratio 
/?„ / h ^.,„ . The magnitude of / /i„, is related to drop in the water level betvft,'een upstream 
and downstream of the fuse plug. The magnitude increases steeply during the initial 
stages of erosion from 0.08 to 1.0, later on there is gradual increase from 1 to 6. 

b. Washout Profile Characteristics 

Characteristics of washout profiles includes the variation of crest, sediment profile, rate of 

crest erosion. These are related with time and depth scales derived from the rate of reserv'oir 

depletion during washout. 

i. Reservoir time scale (T*) and depth scale {H*} are measured from rate of depletion of 
reservoir having zero to negative maximum. 

ii. Nondimensionalised maximum crest height(/i£.„///^) is found to decrease with time as 

iii. Nondimensionalised sediment bed profile / h„, ) is found to decrease with distance (X) 
measured from the crest position as ih^ /hcn,)= l/(^+(^/^5o/) 

iv. Scales X^q and t^o are related with reservoir depletion scales T* and H, . 

V. Rate of crest erosion (dhjdt)!^, is found to increase initially with time, reaches a 

peak value of 2 x 10"’ and then starts decreasing with further increase in time. This 
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nondimensionahsed erosion rate can be treated as erosion coefficient for fuse plug. The 
magnitude of this erosion coefficient is found to vary with time. 

Vi. The maximum possible crest erosion rate {dh,Jdt) nondimensionalised with sediment 
resisting velocity is related to peak discharge in the form {q^H, with 

power 3.27. This is carried out mainly to indicate the magnitude of rate of erosion for 
different types of soils. 

vii. Sediment erosion velocity V^y nondimensionalised with its maximum value is 

related to nondimensionalised flow velocity (FJ as magnitude of F,.. 

Vsvs, remains fairly constant as 0.9 when « 0.6. With further increase in 

sediment erosion velocity V 3 V decreases very steeply. This variation is 
seen for all the three types of soils used in the fuse plug. 

viii. The sediment erosion velocity scale V^y;^ nondimensionalised with soil resisting 
velocity is found to be a function of number As the magnitude of A,, increases 
from 1 to 10 the magnitude of Vsy^/ Vsr increases from 0.03 to 0.8. The magnitude of 
V^y xi is higher for cohesionless soil and is of the order of 0.8 times . 

c. Sediment Transportation Characteristics 

i. The magnitude of sediment transported at any time Q 5 nondimensionalised with its 
maximum value is found to increase steeply in initial stages of erosion, reaches peak 
and then decreases gradually with increase in time. More data are available in later stages 
of erosion process. The receding limb of the sediment discharge(g 3 ) is analysed with time 

scale and found to follow a good relation as fe / Qw) = i '' (1 + 

X~(t-t5p)/(t5p/2- tjp). 

ii. Sediment concentration q/q^ is related to It is found that q/q^ decreases with 

increase in . The relation between them is given by polynomial equation obtained 

from regression analysis. 

iii. Sediment concentration was found to be almost independent of bed shear velocity. 
However it is found to be a function of drop in the water level(4E). 
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d. Energy Loss Characteristics 

I . The ratio of energy at downstream end to the upstream of the fuse plug, £/£, is related to 

hjh,„ for simulated rigid bed washout profiles. It is found that as hjh,^ increases £y£. 
also increases. 

II . Effect of erosion of the fuse plug material on the energy ratio E^E, is related to h £ 

^ i I/' QfJl 

along with rigid bed profile data. It was found that considerable energy loss occurs due to 
erosion of sediment at given state of flow. 

iii.The energy loss due to erosion and transportation of sediment is found fi'om the relation 
[(E:) R.g.d-( E2)Mobiie]/Ei against sediment concentration Q^Qw . It is found that the energy- 
loss increases very steeply to a value of 0.9 with increase in sediment concentration to a 
magnitude upto 8 percent. With further increase in energy loss it is found to be nominal. 
Flow with sediment concentration higher than 8 percent is considered as hyperconcentrated 
flow. 

e. Hydrograph and Sediment Discharge 

i. The normalised sediment discharge intensity with its peak value q^p is found to follow 
closely with time ratio t/t^p as {q^ I q,p) = X' exp(l - X‘), where X= t/t^p . 

ii. The peak sediment discharge q^p is found to be a function of q^^yf 

iii. The time scale t^p nondimensionalised with t^p, the time at which maximum water 
discharge occurs is found to decrease with peak water discharge intensity as q'^p 

iv. Duration of washout process of fuse plug designated as % is found to be a function of 
reservoir time scale £♦, inflow discharge, soil resisting velocity and height of the fuse 
plug. The magnitude of % varies firom 1 to 2.5 time T* . 

f. Outflow Hydrograph Characteristics 

i. The value of q^q^yj, plotted against i/i^p is found to increase initially, reaches peak and 
then decreases slowly to a value of mcoming discharge ratio. The occurrence of peak 
discharge depends on the type of soil eroded. For cohesionless soil, the peak occurs 
almost in initial stage of washout process, whereas for cohesive soil fill it occurs at the 
later part of washout process. Because of this nature of vanation, analysis of nsmg limb 
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of hydrograph is earned out in detail. All the rising limbs follow a relation 
(^vv ^ exp(l~Ar ), where The magnitude of exponent ‘m’ for 

cohesionless soil is 2.5 and it is 0.5 for cohesive soil. 

ii. The magnitude of peak water discharge is found to increase with increase in 
magnitude of inflow discharge with a power of 0.45. 

iii. Time scale is the time at which peak water discharge occurs. This scale is found to 
decrease reaching a minimum value and then increases with increase in peak discharge 

. 


8.4 Suggestion for Further Work 

a. Experiments should be conducted for dam breach study as well as for fuse plug with 
temporal variation of the inflow discharge to the reservoir, taking its average 0 :„ , more 
variation of dam geometry and its fill properties. The present developed functional 
relationships can be verified with the breach parameters obtained from these variations. 

b. The automation of data acquisition in terms of time and space should be planned for better 
accuracy of the data. 
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APPENDIX 


Sample Calculation for x*c 

Fine cohesionless Soil (Type-A soil) : 

d 5 o= 0.16 mm 


K. 


Ap 


= 0.27 


gd 


50 


Where = Critieal Shear Velocity 


(t ) +(t ) 

\ / cohesion V / 


friction 


Apgd 


= ( 027 )^“ 


50 


= 0.0729 = r.. 


[For type-A soil, = 0.0] 



Calculation for for soil fill of Teton Dam: 


djo = 3 mm, ^ = 40°, c = 49 kPa 


Ap 


Sd 


50 


= 020976 


d=‘imm 


: (0.20976)' = 0.044 


^Pgd so 

(0-044) X (1.65 X 10^ X 9.8 
= (0.044) X 4856 = 2.1366 kgf/m^ 


( ^ or ) roteive =1-1 5x49 = 56.35 
= 5.74 kgf I 

Combined Shear stress = 2.1366 + 5.74 = 7.88 kgf7m' 




Combined 


7.88 


i^P§dso) 48.56 


= 0.162 




Variation of nondimensional bed shear stress with particle 
size of soil(mnn) obtained fronn Shield's curve 


128596 


128596 
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